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CHAPTER  1 
INTRODUCTION 

Signal  processing  techniques  have  definite  applications  in 
the  area  of  microprocessor  or  computer-based  instrumentation.  A 
new  approach  has  been  evolved  wherein  a  microprocessor-based 
system  can  accomplish  data  acquisition  with  subsequent  processing 
by  a  computer  utilising  suitable  processing  algorithms.  The 
algorithms  used  depend  on  what  exactly  needs  to  be  measured  or 
determined  -  a  specific  example  considered  here  being  the 
measurement  of  phase.  Such  a  system  possesses  the  inherent 
advantages  of  speed,  automation,  flexibility  in  type  of  algorithm 
used  for  the  specific  problem,  and  enormous  computing  and 
processing  capabilities  coupled  with  all  the  advantages  which 
discrete  time  signal  processing  can  offer. 

A  microprocessor  based  phasemeter  has  been  developed  which 
acquires  data  samples  of  the  original  signals  and  sends  them  to  a 
computer  for  subsequent  processing  by  suitable  algorithms.  These 
algorithms  are  based  on  Adaptive  Filtering,  FFT,  Correlation  and 
Numerical  methods.  The  system  was  implemented  on  a  North  Star 
microcomputer  in  conjunction  with  an  HC  6802  interface  and  HP 
Multiprogrammer;  the  whole  forming  an  automated  test  system.  The 
MC  6802  interface  serves  as  a  communication  link  between  the 
North  Star  (NS)  and  HP  Multiprogrammer  (HPM).  The  sampling  of  the 
signals  is  initiated  by  a  command  from  the  NS  software  which 
invokes  an  assembly  routine  in  the  MC  6802  interface  and  uses  the 
High  Speed  Analog  to  Digital  converter  (ADC)  card  on  the  HPM. 
These  samples  are  then  sent  serially  to  the  NS  and  subsequent 


processing  takes  place  by  the  algorithms  stated  above  to 
determine  the  phase  difference.  The  HP  9845B  Desktop  computer  was 
used  to  study  theoretical  simulation  results  and  for  the 
development  of  algorithms. 

Conventional  phasemeters  operate  on  the  principle  of 
measuring  the  time  difference  between  the  zero  crossings  of  the 
two  signals  by  "squaring  up"  the  signals.  Complex  and 
sophisticated  hardware  consisting  of  amplifiers,  precision 
reference  generators,  filters  and  networks  are  normally  used. 
Also,  manual  controls  and  adjustments  are  required  for  making 
measurements.  The  automated  microprocessor  based  system,  however, 
employs  more  software  and  less  sophisticated  hardware.  Though 
this  system  may  not  boast  of  the  same  accuracy  obtainable  with 
the  dedicated  and  sophisticated  phasemeter,  it  still  has  certain 
advantages  and  is  capable  of  further  development  and  improvement. 
The  relative  merits  and  demerits  are  discussed  later. 

Other  applications  can  also  be  found  for  the  system 
discussed  above,  examples  being  spectral  estimation  (PDS), 
frequency  domain  analyses,  filtering,  combating  noise  and 
determination  of  RMS  value  etc.  Thus  this  project  represents  an 
initial  step  in  a  novel  direction  and  the  author  is  optimistic 
that  further  development  and  sophistication  is  possible  in 
future . 


CMEXEB  2 
SOME  PRELIMINARY  CONSIDERATIONS 

The  specific  application  of  signal  processing  techniques  in 
phase  measurement  has  been  selected.  There  are  some  questions  to 
be  answered  (before  we  proceed  further)  such  as  (a)  Why  measure 
phase?  (b)  What  are  the  conventional  methods  of  measurement?  (c) 
How  are  the  proposed  methods  different  and  what  advantages  do 
they  have  over  the  conventional  methods. 

The  measurement  of  phase  is  accomplished  by  a  phasemeter. 
Such  measurements  find  applications  in  impedance  meters,  the 
measurement  of  amplifier  phase  shift,  the  transfer  functions  of 
networks,  the  measurement  of  group  delay,  the  measurement  of 
distance  and  directivity,  and  the  testing  and  adjusting  of  filter 
networks. 

The  principle  of  phase  measurement  used  in  modern  phase 
instrumentation  is  illustrated  in  fig.  2.1.  through  fig.  2.4.  If 
"tn  is  the  time  interval  between  analogous  zero  crossings,  and 
T"  is  the  wave  period,  then  phase  angle  P=360*(t/T)  degrees.  The 
signals  are  then  squared  up  as  shown  in  fig.  2.2.  Most 
phasemeters  use  the  squared  off  waves  to  generate  the  waveshape 
shown  in  fig.  2.3.  The  Dranetz  305  is  one  of  a  series  of 
sophisticated  phasemeters  available.  This  phasemeter,  after 
generating  the  waveform  in  fig.  2.3.  standardizes  the  waveform 
amplitude  and  then  filters  the  waveform  to  yield  the  DC  value. 

The  filtering  in  the  simplest  case  is  accomplished  by  a  low 
pass  filter  ie.  an  RC  network.  The  time  constant  of  this  filter 
must  be  large  enough  to  provide  an  accurate  measurement. 


Fig. 2.1.  Signals  with  a  certain  phase  difference. 


Fig. 2. 2.  Clipning  or  squaring  up  the  signals. 
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Fig. 2. 3.    Triggering  on  and  off  at   the   zero  crossing 
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Fig. 2. 4.  Standardizing  the  vaveform  amplitude. 


Fig. 2. 1. . .2.4.  Principle  of  phase  measurement. 
Vref,   and  all  DC  amplifiers  following  this  circuit  and  other 


circuitry  must  be  stable  and  accurate  with  time  and  temperature. 
The  speed  of  response  is  limited  by  the  time  constant  and  hence 
in  practice  it  is  made  variable  by  automatic  switching  of  filters 
with  different  time  constants. 

The  presence  of  noise,  interference,  high  harmonic  content 
or  various  forms  of  signal  modulation  may  impair  accuracy.  Also 
the  input  signals  may  vary  appreciably  in  magnitude.  This 
variation  is  counteracted  by  switching  attenuators.  Thus  in  the 
Dranetz  system  it  is  found  that  a  great  degree  of  complexity  and 
sophistication  in  circuitry  is  used  to  combat  all  the  above 
mentioned  problems. 

Having  considered  the  conventional  method  of  phase 
measurement  attention  is  now  focussed  on  the  proposed 
techniques.Here  the  analysis  of  signals  is  carried  out  in  the 
discrete  time  (DT)  domain  where  we  work  with  "samples"  of  the 
original  signals  as  shown  in  fig.   2.5. 


Fig. 2. 5.  Samples  of  the  signals  in  the  DT  domain. 
It   is  necessary  according  to  the  Nyquist  criterion  that  if 


all  information  in  the  original  signals  has  to  be  preserved  then 
the  sampling  frequency  "f_"  must  be  greater  than  the  frequency  of 
the  signals  "fi_"  that  is  fs>2*fQ.  Once  the  signals  are 
characterized  by  the  "samples"  then  the  next  step  is  to  develop 
some  algorithm  to  process  these  samples  to  obtain  the  required 
characteristic  (for  example  the  phase)  of  the  original  signals. 
These  samples  may  be  obtained  by  using  the  Analog  to  Digital 
Converter   (ADC). 

The  ADC  produces  a  digital  output  proportional  to  the  analog 
input.  The  parameters  of  importance  here  are  the  conversion  time 
of  the  ADC  (how  fast  it  can  sample)  and  its  resolution  (how  well 
or  accurately  it  represents  the  analog  quantity  in  digital  form). 
Some  ADCs  that  are  available  have  8  or  12  bits  resolution  and 
are  capable  of  moderate  sampling  rates.  The  "Flash"  type  of  ADC 
is  available  with  a  high  sampling  rate  (10  Mhz.)  and  8  bits 
resolution.  A  more  detailed  discussion  on  ADCs  is  found  later. 

It  is  necessary  at  this  stage  to  emphasize  the  fact  that 
once  the  signals  are  "sampled"  or  "digitized"  and  the  resulting 
samples  are  stored  in  memory  then  they  are  absolutely  unaffected 
by  any  time  or  temperature  changes.  Hence  the  only  caution  that 
must  be  exercised  is  prior  to  and  during  digitizing  when  the 
effects  of  noise,  temperature  changes,  offsets  and  all  the 
problems  associated  with  hardware  exist  (as  in  the  conventional 
phase  meter).  Thus  the  single  greatest  advantage  is  due  to  the 
use  of  software  algorithms  which  are  unaffected  by  the  problems 
that  affect  hardware.  A  computer  can  then  obtain  the  samples 


stored  in  memory  and  then  process  them  with  the  pertinent 
algorithm  and  obtain  the  required  result.  All  the  advantages  of 
software  such  as  flexibility,  use  of  statistical  methods, 
computing  capabilities  become  apparent. 

The  techniques  discussed  use  Adaptive  Filtering  which  uses 
digital  filters  (in  contrast  to  the  analog  RC  filters  in  the 
conventional  phasemeter).  Here  the  coefficients  of  the  filters 
can  be  programmed  and  changes  can  be  rapidly  made  without 
changing  circuit  components.  Performance  is  not  affected  by  wear 
and  aging  and  is  more  reliable.  The  response  characteristics  can 
be  easily  changed  by  manipulating  the  coefficients  which  in  turn 
leads  to  more  flexibility.  However,  the  finite  resolution  of  the 
ADC,  finite  precision  or  word  length  of  the  computer,  truncation 
and  roundoff  errors  and  overflow/underflow  problems  are 
associated  with  software  processing.  Nevertheless,  by  proper 
hardware/software  tradeoffs,  it  is  possible  to  obtain  a  suitable 
operating  system. 

Thus  with  these  preliminary  considerations  we  now  proceed  to 
describe  the  proposed  algorithms. 


CHAPTER  2 
THE  ADAPTIVE  FILTERING  ALGORITHM 
3.1   INTRODUCTION 

This  method  uses  a  signal  processing  technique  familiarly 
known  as  the  Widrow's  LMS  (least  mean  square)  algorithm  in  the 
cancellation  mode.  The  schematic  is  shown  in  fig  3.1.  The  two 
signals  y(k)  and  x(k)  form  the  primary  and  reference  inputs  to 
the  adaptive  filter  configuration.  Here  y(k)  and  x(k)  are  derived 
by  sampling  the  CT  (continuous  time)  signals  x(t)  and  y(t)  and 
"k"  is  the  discrete  time  index  which  represents  the  sample  of  the 
original  signals  at  time  t=k.  The  error  e(k)  that  is  generated  by 
comparing  g(k)  and  y(k)  is  used  to  change  the  weights  of  the 
discrete  time  (DT)  digital  filter  W  such  that  the  square  of  the 
mean  error  e(k)  is  a  minimum.  Thus  as  the  processing  continues 
the  error  dynamically  changes  the  weighting  coefficients  of  the 
filter  so  that  the  error  reduces  in  magnitude.  After  a  sufficient 
number  of  processing  steps  the  error  tends  to  a  small  value,  the 
weights  of  the  filter  are  then  used  to  calculate  the  phase 
difference  between  y(k)  and  x(k).  Here  y(k)  is  assumed  to  lead 
x(k). 

3.2   THEORETICAL  ASPECTS 

Consider  the  adaptive  scheme  in  fig.  3.1.  It  is  seen  that 
y(k)  and  x(k)  form  the  primary  and  reference  inputs  and  Wk  is  a 
M+l  weight  filter.  The  error  e(k)  obtained  by  comparing  y(k)  with 
the  filter  output  g(k)  is  used  to  update  the  filter  weights  or 
coefficients.   This  process  is  carried  on  until  the  error   reduces 


to  a  small  value. 
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Fig. 3.1.  Adaptive  filtering  scheme  in  the  cancellation  mode. 
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Hence  we  have  the  following,   in  vector  form: 
W  =  (wQ  Wj  Wj  •••  wm>  considering  Mfl  weights. 
Xk'  =(x(k)   x(k-l)   x(k-2)    ...  x(k-M)) 

where  x(k-l)...  x(k-M)   represent  the  M  past  values  of  x(k). 
x(k)  =  sin(wQkT)  and  y(k)  =  sin (wQkT+P) 
Now  g(k)=  W'Xk  =  w0x(k)+w1x(k-l)+...wMx(k-M) 
e(k)   =  y(k)-g(k) 

When  the  error  e(k)  tends  to  a  small  value  as  a  result  of 
the   adaptive   process   then  g(k)    tends   to   cancel   y(k)    at    the 
summing  device.    Hence  g(k)    "looks  like"  y(k).   At  this  stage 
taking  the  Z  transform  to  obtain  the  filter  transfer  function  we 
have  the  following: 
G(z)   =  X(z)[w0+w1z"1+...+wMz"M] 
The  transfer  function  H(z)   is  then 
H(z)=G(z)A(z)   and  since  z=eJv?'  where  w0=271f0 

f0  being  the  frequency  of  the  signals  and  T  being  the  sampling 
period. 

Hence  the  value  of  transfer  function  at  the  frequency  of  the 
signals  is  given  by 

H(w0)   =  G(w0)A(w.)   =w0+w1e"iw?+...+wMe"3Ml<1, 
The  phase  transfer  function  B(w0)=arctan[   ImH(w0  )/ReH(w0)] 
ImH(w0)  =  [w1sinw0T+w2sin2w(>T+...+wHsinMw0TJ 
ReH(w„)  =  [w0+w1cosw0T+w2cos2w0T+...+wHcosMw0T] 

The  value  of  B  at  wc  gives  the  required  phase  difference  "P" 
between  the  two  signals  x(k)   and  y(k) . 

Thus  from  the  above  discussion  it  is  apparent  that  after  a 
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sufficient  number  of  processing  steps  (or  iterations)  when  the 
error  decreases  the  phase  difference  "P"  can  be  directly 
determined  from  the  coefficients  of  the  fiter  W.  It  must  be  noted 
that  in  the  above  expression  for  B  it  is  required  to  know  the 
values  of  fQ  and  T  for  calculating  the  phase  P. 

3.3  SIMULATION  AND  VERIFICATION 

As  a  first  step,  software  development  for  the  method  was 
done  on  the  North  Star  microcomputer  system.  The  two  signals  x(k) 
and  y(k)  were  generated  as  x(k)=sin(wckT)  and  y{k)=sin(w0kT+P) 
where  w  =2KfQ,  fQ  being  the  frequency  of  the  signals,  T  the 
sampling  period,  k  the  time  index  and  P  the  phase  difference  in 
radians.  The  Widrow  algorithm  was  implemented  in  software.  The 
flow  chart  and  BASIC  program  listing  are  found  in  Appendix  la  and 
lb  respectively. 

The  procedure  followed  is  essentially  iterative  in  nature 
and  the  weights  are  updated  as  follows  for  a  M+l  weight  filter. 
wm(k+l)=wm(k)+ve(k)x(k-m) ;         0<=m<=M+l 

where  v  is  the  convergence  parameter  and  this  parameter  is 
dynamically  updated  as  follows: 

v=«./oi1(k);     tfxJ(k)  =  p^(k-l)+(l-p)x2(k);       <r/(0)=0 

The  parameter  oc  must  be  chosen  with  care  for  proper 
convergence  and  the  phase  is  calculated  as  P  given  by  the  formula 
for  B  in  section  3.2. 

The  most  important  parameters  that  are  involved  in  this 
method  are  the  following: 
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a.  the  frequency  of  the  signals  (fQ) 

b.  the  sampling  frequency  (fs=l/T) 

c.  the  convergence  parameter  (a.) 

d.  the  phase  difference  involved  (P) 

e.  the  number  of  iterations  (K) 

f.  the  ratio  (Vfo' 

Some  of  the  important  results  obtained  after  considerable 
experimentation  with  the  parameters  stated  above  are  as  follows. 

1.  The  theoretical  value  of  sampling  frequency  fs>2fQ,  but  it 
was  found  that  fs>=3fQ  for  proper  working  of  the  method. 

2.  For  a  3  weight  filter  and  with  just  30  iterations  convergence 
was  established  when  the  value  of  sampling  frequency  was  equal  to 
four  times  the  signal  frequency  ie.  fs=4fQ  and  the  calculated 
phase  agreed  very  closely  with  the  actual  phase  (absolute  error 
less  than  0.2  degrees.) 

3.  However  with  a  16  weight  filter  convergence  was  easily 
obtained  for  a  wide  range  of  fg/f0- 

4.  One  of  the  theoretical  aspects  pertaining  to  convergence  is 
that  a  greater  number  of  weights  leads  to  faster  convergence. 
Also  the  convergence  time  is  also  dependent  on  the  convergence 
parameter       as    T=l/*  =1/(1-?). 

5.  Also  the  faster  the  convergence  the  greater  is  the  residual 
error.  Too  fast  a  convergence  can  overdrive  the  filter  W  causing 
large  residual  errors,    (ie.  magnitude  of  e(k)    ). 

6.  Hence  the  number  of  weights  and  the  convergence  parameter 
must  be  chosen  carefully  to  obtain  a  small  residual  error  and 
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hence  a  small  error  in  the  calculated  phase. 

7.  The  convergence  also  depends  on  the  type  of  input. 

8.  The  method  works  for  any  two  zero  mean  signals. 

9.  Since  the  filter  by  itself  cannot  distinguish  a  phase  shift  of 
more  than  90°  suitable  modifications  were  made  in  the  software 
which  fixes  the  quadrant  with  respect  to  the  sign  of  the 
numerator  and  denominator  in  the  equation  for  phase  as  already 
shown.  A  phase  of  360°  could  then  be  easily  distinguished. 

10.  The  above  results  were  established  assuming  y(k)   leads  x(k) . 

Analysing  the  working  of  the  method  in  the  frequency  domain, 
it  can  be  said  that  the  adaptive  filter  merely  adjusts  the 
amplitude  and  phase  of  all  the  components  in  its  input  signal 
x(k)  to  cancel  those  in  y(k).  This  is  a  characteristic  of  the 
adaptive  process. 

3.4  DETERMINATION  OF  PHASE  SIGN 

In  the  above  discussion  it  was  assumed  that  y(t)  leads  x(t) 
by  a  phase  difference  P.  When  this  is  not  known  we  resort  to  the 
following  method  to  determine  which  is  the  leading  or  lagging 
waveform.  If  during  the  period  x(t)  is  positive  there  occurs  a 
negative  to  positive  transition  of  y(t)  then  x(t)  leads  y(t)  and 
if  there  occurs  a  positive  to  negative  transition  of  y(t)  then 
y(t)   leads  x(t)  . 

However  the  sign  is  automatically  reflected  in  the  present 
software.  If  y(k)  and  x(k)  are  interchanged  (assuming  y(k)  leads 
x(k)  )  then  the  calculated  phase  shows  a  value  of  360-P  where  P 
is  the  actual  phase  difference.  Numerical  methods  presented  later 
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can  be  used  to  determine  to  a  moderate  degree  of  accuracy  the 
zero  crossing  from  which  the  sign  could  be  determined. 
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Fig. 3. 2.  Determination  of  phase  sign. 

14 


3.5  PRACTICAL  ASPECTS  AND  SYSTEM  IMPLEMENTATION 

After  having  established  the  validity  of  the  method  in 
theory,  it  was  extended  to  a  practical  situation.  However,  as  is 
always  the  case  ,  the  transition  from  theory  to  practice  is  not 
smooth  and  many  constraints  had  to  be  given  due  importance. 

Before  we  discuss  the  practical  system  we  must  bear  in  mind 
the  following  aspects  of  the  theoretical  simulation  carried  out. 

1.  The  samples  generated  are  correct  to  7  decimals. 

2.  f0  and  T=l/fs  are  quite  accurate  since  they  are  just  values 
given  by  us. 

3.  No  noise  is  included  in  the  signals. 

4.  Phase  shift  is  constant  with  time. 

However  in  a  practical  system,  the  situation  is  a  bit 
different.  The  finite  resolution  of  the  ADC  imposes  a  constraint 
on  the  accuracy  with  which  a  sample  is  represented.  For  example 
with  a  12  bit  ADC  the  accuracy  is  limited  to  around  3  decimal 
places.  The  frequency  of  the  signals  must  be  determined  by  some 
external  means  and  the  sampling  frequency  must  be  accurately 
known.  The  signals  may  be  derived  from  a  function  generator  which 
has  some  noise  in  its  output  and  may  also  have  some  voltage 
offsets.  The  phase  of  the  waveform  generated  may  drift  with  time. 
Hence  all  these  must  be  considered  in  a  practical  system. 

The  basic  configuration  required  for  system  implementation 
of  the  method  is  shown  in  fig. 3. 3. 

The  ADCs  must  be  controlled  to  sample  both  the  signals  x(t.) 
and  y(t)  simultaneously.  These  samples  may  then  be  stored  in 
memory.  Once  stored  in  memory  the  samples  can  be  accessed  by  a 
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computer  or  microcomputer  for  processing  the  samples  with  the 
adaptive  algorithm.  The  CRT  terminal  acts  as  an  Input/Output 
device  to  provide  interaction  between  the  user  and  the  system. 

There  are  different  ways  of  realizing  the  above  system  but 
the  primary  efforts  were  directed/ to  utilise  the  facilities 
available.  A  practical  system  using  the  North  Star  microcomputer 
in  conjunction  with  the  MC  6802  interface  and  HP  Multiprogrammer 
(HPM)  being  used  as  an  automated  test  system.  A  brief  discussion 
of  the  HPM  and  the  interface  follows. (fig. 3. 3. a.) 
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Fig. 3. 3.  System  implementation  of  the  adaptive  method. 
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17 


The  HPM  is  a  master  control  unit  for  bidirectional  systems. 
It  can  be  controlled  using  a  computer.  It  contains  many  types  of 
accessory  cards  which  can  be  used  depending  on  the  requirement. 
For  example  the  V/F  card  in  the  HPM  can  be  programmed  to  make 
voltage  or  meaurements  on  the  device  under  test  (DUT).  The  HPM 
also  contains  a  High  Speed  Analog  to  Digital  Converter  card 
(ADC).  Each  card  can  be  accessed  by  addressing  the  particular 
slot  in  which  it  resides.  The  various  communication  protocols  are 
provided  by  the  following  four  types  of  words  -  control  word, 
output  data  word,  input  address  word  and  return  data  word. 

The  MC  6802  interface  acts  as  a  communication  link  between 
the  North  Star  (NS)  and  the  HPM.  There  exists  an  assembly  routine 
in  the  interface  EPROM  that  can  address  the  ADC  in  the  HPM, 
obtain  the  samples  and  store  them  in  the  memory  available  in  the 
interface  RAM  (IK).  The  procedure  is  to  first  set-up  the  HPM  in 
the  input  mode.  Then  the  ADC  card  is  accessed  by  the  proper  slot 
address.  The  conversions  are  initiated  by  a  gate  signal  and  the 
end  of  conversion  is  indicated  by  a  flag.  The  routine  makes  use 
of  the  above  procedure  to  obtain  the  samples. 

The  ADC  card  is  capable  of  22K  conversions  per  second. 
However  due  to  the  inherent  limitations  posed  by  the  software 
assembly  routine  the  present  sampling  frequency  is  9130  samples 
per  second.  The  number  of  samples  that  can  be  taken  is  fixed  at 
512.  The  ADC  has  three  operating  ranges  (±10V,  ±1V,  ±100mV).  The 
interface  accomplishes  sampling  of  the  signal  and  stores  in 
memory.  It  then  sends  the  data  samples  to  the  NS  using  the  serial 
port.  Thus  the  sampling  is  initiated  by  a  command  in  the  BASIC 
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program  which  invokes  the  assembly  routine  in  the  interface  EPRDM 
and  the  samples  are  then  sent  to  the  NS  memory.  The  schematic  of 
the  interface  and  connections  are  shown  in  Appendix  2. 

A  look  at  fig.  3.3  shows  that  the  two  signals  have  to  be 
sampled  simultaneously  requiring  two  ADCs.  Since  one  ADC  card  is 
available  on  the  HPM  the  following  method  was  adopted  to  enable 
sampling  of  the  two  signals  by  one  ADC  and  yet  make  it  look  like 
they  were  sampled  simultaneously.  The  principle  is  illustrated  in 
fig. 3.4  and  3.5. 
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Fig. 3. 4.  Generation  of  a  trigger  reference. 
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Fig. 3. 5.  Use  of  trigger  reference  to  sample  x(t)  and  y(t) , 
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The  property  of  the  two  signals  that  is  exploited  in  this 
method  is  that  x(t)  and  y(t)  continuous  time  signals  and  are 
periodic.  Hence  the  signals  by  themselves  do  not  change  with  time 
and  with  respect  to  themselves,  a  trigger  reference  may  be  used 
to  initiate  the  sampling  of  the  two  signals. 

In  order  to  obtain  the  trigger  reference  a  zero  crossing 
detector  was  designed  and  built. (fig. 3.6) . 
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Fig. 3. 6.   Circuit  diagram  of  the  zero  crossing  detector. 
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The  output  of  the  zero  crossing  detector  (ZCD)  was  used  to 
initiate  the  sampling  of  x(t)  and  y(t).  In  the  assembly  routine 
just  before  the  sampling  is  actually  initiated,  a  loop  was 
interposed  as  shown  below. 

ADC  assembly  routine 


BACK     LDAA  CRB2 

BPL  BACK 

LDAA  DDRB2 


Sampling  initiated 

We  shall  now  discuss  the  operation  of  the  above  loop.  A 
glance  at  the  schematic  of  the  interface  (Appendix  2)  shows  that 
the  CB2  terminal  (control  register  B,  [pin  18]  )  of  PIA  2  in  the 
interface  can  be  used  as  an  input  terminal  to  which  the  negative 
pulses  of  the  ZCD  are  given.  The  PIA  2  is  initialized  such  that 
in  the  control  register  CRB2  bit  7  is  set  by  negative  transitions 
(♦)  on  the  input  terminal.  Thus  the  output  of  the  ZCD  (the 
transitions  )  triggers  the  sampling  as  the  execution  of  the 
loop  is  terminated  when  bit  7  is  set  by  a  negative  transition. 
The  MPU  read  of  data  register  B  clears  bit  7  and  the  execution  of 
the  sampling  process  is  done. 

However,  this  method  was  later  found  to  have  certain 
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disadvantages  since  it  introduced  some  errors.  The  MPU  clock 
cycle  time  was  found  to  be  1.1  us  since  it  operates  at  920.8  KHz. 
There  is  an  ambiguity  involved  since  the  trigger  pulse  could 
arrive  at  any  instant  with  respect  to  the  two  instructions 
forming  the  loop.  An  analysis  will  reveal  that  there  could  be  a 
maximum  error  of  4  clock  cycles  (4.4  us)  which  at  a  signal 
frequency  of  around  1000  Hz.  gives  rise  to  a  phase  error  of 
around  2.5  degrees. 

Hence  in  order  to  get  around  this  problem  an  interrupt 
service  routine  was  used  in  the  assembly  program  and  the 
properties  of  the  Wait  For  Interrupt  (WAI)  were  utilised.  When 
this  instruction  is  executed  then  the  uP  waits  indefinitely  until 
as  interrupt  occurs  and  then  proceeds  to  the  interrupt  service 
routine  without  any  delay  after  which  the  sampling  is  initiated. 
Thus  the  source  of  ambiguity  and  error  were  removed  totally.  The 
assembly  listing  is  given  in  Appendix  3. 

The  signals  were  derived  using  a  function  generator.  In 
order  to  know  the  frequency  of  the  signals  a  digital  frequency 
counter  was  used.  The  ADC  was  switched  to  operate  in  the  ±10  V 
range  keeping  in  view  the  effects  of  noise.  It  can  be  seen  that 
if  a  sinewave  has  to  be  digitized  properly  then  the  amplitude  of 
the  signal  must  be  adequate  to  cover  the  entire  operating  range 
of  the  ADC.  This  is  discussed  in  a  later  chapter. 

The  signals  x(t)  and  y(t)  were  derived  using  a  simple  RC 
circuit  (fig.3.8).  The  phase  shift  of  this  circuit  is  given  by 
arctan(l/2nf  rq.  This  circuit  also  produces  attenuation  besides 
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the  phase  shift.  Hence  when  y(t)  is  sampled  the  input  is 
increased  to  a  sufficient  value  so  that  the  amplitude  of  y(t) 
covers  the  operating  range  of  the  ADC. 
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Fig. 3. 8.  Use  of  RC  circuit  to  generate  x(t)  and  y(t) . 
Other  considerations  are  that  the  HP3311  (function 
generator)  produces  a  maximum  sine  output  of  +5  V  and  hence  an 
amplifier  using  an  OPAMP  was  used.  Also  since  the  output  of  the 
HP3311  was  noisy  a  low  pass  filter  (LPF)  was  used  at  the  input 
stage.  Offset  adjustment  was  provided  for  the  amplifier,  since 
the  method  requires  that  x(t)  and  y(t)  be  zero  mean  signals.  The 
circuit  diagram  is  shown  in  fig. 3. 9. 
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Fig. 3. 9.  Noise  reduction  and  signal  conditioning  circuit. 
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Accomplishing  all  of  the  above  enabled  a  working  system  to 
be  developed  which  uses  the  adaptive  algorithm.  A  block  diagram 
of  the  system  implementation  is  shown  in  fig.  3.10. 
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Fig. 3. 10.  Block  diagram  of  the  test  system. 
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3.6  WORKING  OF  THE  METHOD  AND  THE  TEST  SYSTEM 

A  schematic  of  the  test  system  is  shown  in  fig.  3.11. 
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Fig. 3. 11.  Configuration  of  the  test  system. 
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The  Silent  700  terminal  acts  as  the  communicating  device 
between  the  operator  and  the  system.  The  entire  operation  is 
controlled  by  the  software  in  BASIC.  A  command  is  sent  via  the 
serial  port  of  the  NS  to  the  MC  6802  Interface  .  This  consists  of 
the  character  "A".  The  main  assembly  routine  in  the  interface 
recognises  this  as  a  command  to  invoke  the  ADC  routine  and 
performs  the  sampling  using  this  routine.  The  samples  are  then 
sent  back  via  the  serial  port  to  the  NS.  The  signals  y(t)  and 
x(t)  are  sampled  successively.  Before  the  processing  starts  the 
frequency  of  the  signals  indicated  by  the  frequency  is  input  to 
the  NS  by  the  use  of  the  INPUT  statement.  The  processing  is  then 
done  and  the  results  are  displayed  on  the  Silent  terminal.  The 
CRD  was  used  to  monitor  the  amplitudes  of  the  signals  to  suit  the 
range  of  the  ADC. 

A  small  modification  that  was  carried  out  at  a  later  stage 
was  that  although  the  ZCD  was  designed  and  used,  the  performance 
of  the  comparator  was  not  up  to  the  mark  (probably  due  to  noise) 
and  hence  a  more  accurate  trigger  reference  was  derived  from  the 
pulse  output  terminals  of  the  HP3311.  (this  being  generated 
internally  in  the  HP3311). 

After  considerable  experimentation,  it  was  found  that  the 
ratio  fg/f0  was  very  significant  in  determining  convergence.  The 
sampling  frequency  being  fixed  at  9130  Hz.,  the  input  signal 
frequency  was  varied  to  study  the  effects  of  the  ratio  over 
convergence  and  the  calculated  value  of  phase.  This  was  done  by 
comparing  the  phase  obtained  with  the  actual  phase  obtained  by 
the    formula    arctan(l/2n  f0  RC)    the   values    of   R   and    C   being 
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determined    to     some    degree    of    accuracy    by    a    digital 
resistance/capacitance  meter. (Digibridge) . 

As  discussed  earlier  in  the  section  3.2,  the  two  parameters 
N  (the  number  of  weights),  and  the  convergence  parameter  played 
a  significant  role  on  the  convergence  (K)  and  the  calculated 
phase  (P).  The  table  below  summarises  the  values  of  N,  oc  ,  and  K 
(the  number  of  iterations  needed  for  convergence)  that  enabled 
proper  working  of  the  method.  These  were  incorporated  in  the 
software     to  cover  the  required  frequency  range. 

I6ELE  SHOWING  Tjffi  VALUES  QE  THE  NUMBER  Q£  FILTER  WEIGHTS 
:  AM!  ITERATIONS  EQE  2JE  COVERED 
FREQUENCY  BflNGE. 

CONVERGENCE  PARAMETER  ITERATIONS 

VI  K 

0.006  200 

0.006  200 

0.006  200 

0.006  200 

0.006  200 

0.004  300 

0.004  300 

0.003  400 

0.003  400 

0.003  400 


29 


FREQUENCY 

WEIGHTS 

F 

N 

200 

8 

300 

8 

400 

8 

500 

8 

600 

8 

700 

4 

800 

4 

900 

4 

1000 

4 

2000 

3 

Thus  the  proper  values  of  N,oC,  and  K  enable  the  filter  to 
converge  smoothly  to  the  actual  phase.  The  range  covered  in  the 
table  is  from  2000  Hz.  to  100  Hz.  and  the  obtained  readings  are 
shown  in  the  Appendix  4. 

Frequencies  less  than  100  Hz.  could  not  be  considered  since 
considerable  attenuation  was  produced  by  the  RC  circuit  and  the 
output  amplitude  was  too  low  to  fit  the  range  of  ADC.  This  thus 
provides  a  frequency  range  fs/f0  from  9130/2000  to  9130/100. 
Hence  if  all  the  frequencies  outside  this  range  are  to  be  covered 
then  the  sampling  frequency  must  be  made  variable. 

Some  of  the  highlights  of  system  performance  are  presented 
below: 

1.  The  sampling  frequency  must  be  at  least  four  times  the 
frequency  of  the  signals. 

2.  The  convergence  parameter  and  number  of  weights  determine  the 
number  of  iterations  for  convergence,  as  a  function  of  the  ratio 

W 

3.  The  working  of  the  method  depends  to  a  great  extent  on  the 

performance  of  the  ADC. 

4.  The  amplitude  of  the  input  to  the  ADC  must  be  properly  scaled 
to  suit  the  operating  range  of  the  ADC  for  proper  digitizing. 

5.  The  trigger  reference  must  be  as  accurate  as  possible. 

6.  The  data  acquisition  time  for  each  signal  is  around  30  sec. 
and  the  processing  time  around  100  sec.  Hence  in  order  to 
decrease  the  time  for  one  phase  measurement  the  two  times  must  be 
minimized. 

7.  The  sampling  frequency  must  be  made  variable  to  accomodate  the 
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entire  frequency  range. 

8.  A  12  bit  ADC  has  been  used;  a  lower  resolution     would  increase 
the  number  of  iterations  for  convergence. 

9.  There  is  always  some  residual  error  due  to  the  fact  that  the 
method  does  not  converge  to  the  Wiener  solution. 

The  performance  of  the  adaptive  method  was  quite 
satisfactory.  The  main  sources  of  error  may  be  due  to  the 
characteristics  of  the  input  signal  (noise,  phase  drift  etc.)  the 
requirement  of  an  accurate  trigger  reference  not  being  adequately 
met,  shifts  in  frequency  of  function  generator  with  changes  in 
amplitude,  accuracy  of  sampling  frequency  and  finally  the 
performance  and  repeatability  of  the  ADC. 

It  may  be  concluded  that  this  method  is  capable  of  being 
implemented  and  also  capable  of  reasonable  accuracy  assuming 
allthe  sources  of  error  are  rooted  out  or  minimized.  A  suitable 
system  for  implementation  of  this  method  is  discussed  in  a  later 
chapter. 
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CHAPTER  4 
THE  £FJI  AfiD  CORRELATION  ALGORITHMS 

4.1  INTRODUCTION 

The  Fast  Fourier  Transform  (FFT)  and  Correlation  algorithms 
represent  other  suitable  methods  for  the  determination  of  phase. 
The  FFT  forms  a  very  powerful  tool  in  the  frequency  domain 
analysis  of  a  sampled  signal.  It  is  possible  to  investigate 
certain  characteristics  of  the  original  signal  by  studying  the 
amplitude  and  phase  spectra  which  the  FFT  makes  possible.  The 
Correlation  technique  is  more  or  less  used  along  with  the  FFT  and 
is  also  discussed.  It  is  shown  that  the  FFT  provides  valuable 
information  about  the  sampled  signal  which  can  be  utilised  to 
describe  system  performance  and  other  characteristics  of  the 
signal. 

4.2  THEORY  OF  THE  FFT 

Consider  the  Continuous  time  signal  x(t)  sampled  at  the 
Nyquist  rate  to  generate  the  sequence  x(m).  The  discrete  Fourier 
transform  (DFT)  is  a  Fourier  representation  of  the  given  sequence 
x(m)  0<=m<=N-l  and  is  defined  as  X(n)=  £  x(m)Wn  ,  0<=n<=N-l, 
where  N  is  finite  and  w=e"327,/N  and  X(n)  is  the  nth  DFT 
coefficient.  It  is  seen  that  the  DFT  coefficients  are 
representative  of  the  coefficients  obtained  by  expanding  the 
given  signal  in  a  Fourier  series.   This  is  shown  in  fig.4.1. 

Thus  the  DFT  coefficients  X(n)   represent  the  sampled  signal 

in  the  frequency  domain  as  follows. 

n-i  __ 

X  (n)  =  ±  x  (m)  W™1         0<=n<=N-l 

„=e-J2W 
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which  implies  that  the  DFT  coefficients  yield  the  Fourier 
transform  of  a  given  sequence  at  a  set  of  points  that  are  w  = 
2K/MT  radians/sec  apart,  or  fn  =1/NT  Hz.  apart. 


m 


SAMPLES     OF     JcC-t) 


time 


Fig. 4.1.  An  N  point  DFT  of  a  periodic  DT  signal. 
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Thus  considering  the  DFT  coefficients  X(n)  we  can  see  that 
"n"  is  the  frequency  index  and  that  X(n)  is  in  general  a  complex 
quantity  having  a  real  and  imaginary  part.  Thus  we  define  the 
amplitude  spectrum  as  pn=  |X(n) I 

and  the  phase  spectrum  as    T-„=arctan[Im(X(n))/Re(X(n))  ] 
Thus   we    see    that    in    the    frequency   domain,     X(n)    shows    the 
frequency  content  of  the  original  signal  x(t),  except  that  this 
information    is   not    available    as   a    continuous    function   of 
frequency  but  only  at  discrete  points  given  by  fn  =1/NT  Hz. 

The  FFT  algorithm  is  merely  an  efficient  and  fast  method 
developed  to  compute  the  DFT.  Thus  a  N  point  FFT  of  a  sequence 
x(m)  is  obtained  by  running  an  FFT  of  N  data  points  obtained  by 
sampling  x(t).  If  the  sampling  frequency  is  made  numerically 
equal  to  the  value  of  N  then  we  have  fn  =1/NT  =fg/N  since  T=l/fs 
ie.fn=l  Hz. 

Thus  for  n=0  we  have   |X(0)l  representing  the  DC  component 
for  n=l  we  have   |X(1)|    representing  the  amplitude  of  the  1  Hz. 
component  and  so  on. 

Similiarly  for  the  phase  spectrum  the  phase  of  the  different 
frequency  components  of  the  signal  is  obtained  at  n=0,l,2..  Hz. 

Now  let  us  consider  two  signals  x(t)  andy(t)  where  y(t)  is 
the  leading  signal.  If  x(t)  and  y(t)  are  sine  waves  of  some 
frequency  f0  then  by  running  a  FFT  of  the  samples  of  x(t)  and 
y(t)  we  have  access  to  their  phase  spectra.  We  determine  the 
phase  of  the  frequency  component  fQ  for  x(t)  and  that  for  y(t) 
and  the  difference  yields  the  phase  difference  between  x(t)    and 
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y(t). 

4.3  SIMULATION  AND  VERIFICATION 

The   simulation  was  carried  out  on  the  HP  9845B  Desktop 
computer.   Two   sine  waves  of   frequency  equal   to  24   Hz.    were 
generated  as  follows  and  sampled  at  128  Hz.  ie.fs=128  Hz.  128 
samples  (N=128)  were  taken  for  running  the  FFT. 
x(k)  =  0.5  sin(w0kT) 
y(k)   =2  sin(w0kT+P) 

P=50  degrees  and  its  equivalent  in  radians. in  the  expression  for 
y(k). 

A  128  point  FFT  was  run  on  the  two  sampled  signals  using  the 
AC  waveform  analysis  package  available  on  the  HP  9845B.  The 
amplitude  and  phase  spectra  were  obtained  on  the  plots  shown. 
Also  a  printout  of  the  FFT  coefficients  was  obtained.  The  plots 
and  the  printout  are  shown  figs. 4. 2  through  4.6. 

Let  us  consider  the  different  spectra  in  turn.  A  look  at  the 
amplitude  spectrum  of  x(k)  shows  that  there  is  a  single  frequency 
component  at  24  Hz.  with  a  magnitude  of  0.5  which  is  justified 
since  x(t)  =  sin(2 n*24*t)(.5) 

A  look  at  the  amplitude  spectrum  of  y(k)  shows  a  peak  at  24 
Hz.  again  with  an  amplitude  of  2  which  is  justified  since  y(t)  = 
2sin(2A*24*t+P). 

By  looking  at  the  printout  in  fig.4.6  we  see  that  the  phase 
of  the  24  Hz.  component  for  x(k)  is  -22.5  degrees  and  that  for 
y(k)  is  27.5  degrees  .  Hence  the  phase  difference  is  given  by 
27.5-(-22.5)  =  50  degrees  which  is  exactly  the  phase  difference 
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between  x(t)  and  y(t)  . 

Thus  the  algorithm  was  verified  by  using  the  FFT  package  on 
the  HP  9845B. 


SMPL1TU0E  OF  X 


FREQ.  HINDOH-6  TO  S 
FREQ.  INTERVAL- 1 

nhx  vslue-s . aeeeeE-a I 


COEFFICIENT 


Fig. 4. 2.  Amplitude  spectrum  of  X. 
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FREQUENCY    DOHA  EN    DflTH   *PHHSE> 

FREQ.     HENDOH-a    TO    S.40eeuE~ 
FREO.     ENJ$RVAt_- 1  .  QQQl 


I    8.BBEJ-BB     I.SBE+BI     2.4BE431    3.2BE*ei    4.B0E4B1    4.aBE*ai    3.5BE4B1    S    IDE- 
FREQUENCY    CHzJ 


Fig. 4. 3.   Phase  spectrum  of  X. 
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AMPLITUDE    OF    Y 


freq.   wtNDOw-e  to  s.*eeeeE-C'! 

FREQ.  INTERVAL- 1 . 9©eeaE*0© 

hhx  vhlue-2 . eeeeeE*ee 


a  18  24  J2  40  4g  S3  M 

, , eOCFFffffENT ! | ! 

a.aeE-nae  i.mmi  2.4ae*ai  a LMRtjt) ».aaE*ai  4.aaEfai  s.saEtai  s.«aE*«u 


Fig. 4. 4.  Amplitude  spectrum  of  Y. 
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FREQUENCY    DOMAIN    DATA     (PHRSE) 

FREO.     HIND(JW«S    TO    S  .  4  3Q33E-0  '. 
FREQ.     INTERVAL-]  .  O0QQi2E-r2G 


] .SCE+C1     2.4SE+01    3.20E*gi    4.0OC+01    -t.SBE+Bl    3.ESE+B1    B.lBE-ni 
FREQUENCY    CHr3 


Fig. 4. 5.   Phase  spectrum  of  Y. 
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FREQUENCY  DOMAIN  DATA 


FREQUENCY  UINDOW«0  TO  6. 40000E+01 [Hz] 
FREQUENCY  INTERVAL*  1 . 00000£*00 (Hz  1 


COEFF. 

FREQUENCYtHl! 

REAL 

IMC 

HACNITUDE 

PHASECDEC1 

DC   TERfl 

0.0OO00E*O0 

-1.06250E-1I 

HAX   FREQ. 

6.400001-01 

6.09375E-12 

1 

1 -OOOOOE-OO 

1  .031221-11 

-1.85450E-U 

2.12193E-11 

-60.92 

2 

2.00000E-00 

1.41345E-11 

4.47449E-12 

1  .482S8E-1I 

17.57 

3 

3.00000E-00 

3.49895E-12 

1 . 09691 E- 11 

1.15136E-U 

72.31 

4 

4.  aooooE-oo 

-1  . 09111E-11 

3.22410E-12 

1 .13774E-11 

163.54 

5 

5.  OOOOOE-OO 

7.38308E-12 

-1 .62487E-11 

1 .79319E-11 

-64 , 98 

6 

6. OOOOOE-OO 

1 . 00808E-11 

2.99170E-12 

1 .03134E-11 

16.53 

7 

7.00000E*00 

9.7232SE-12 

-1.81103E-12 

9.89050E-12 

-10.53 

3 

3.00000E-00 

1.S8410E-11 

6. 411411-12 

1 . 99085E-1 1 

18.83 

9 

9.00000E-10 

1.81I76E-11 

2.07836E-11 

2.75731E-11 

43.92 

10 

1  . OOOOOE-01 

-1  .2174*1-1! 

4.06082E-11 

4.23941E-11 

106.69 

11 

1  . 10000E-01 

-3.30S81E-11 

-1  .96893E-13 

3.30387E-U 

-179.66 

12 

1  .20000E-01 

1.01196E-11 

-4.26504E-12 

1  .09817E-U 

-22.85 

13 

1  . 3OO00E+01 

-9.59942E-12 

1.28S73E-11 

1  .  6043SE-U 

126,73 

M 

1  .40000E-01 

-1.8»773E-1I 

7.21818E-12 

2.03037E-11 

139. 13 

15 

1 .500001-01 

-2. 13964E-11 

3.0I387E-12 

2.14076E-11 

171.98 

It 

1  .  600001-01 

-9.97833E-12 

-1.43867E-1! 

1.75084E-11 

-124.74 

17 

1 . 70000E+01 

-3.93310E-12 

2.12129E-13 

3.93879E-12 

176.93 

18 

1  .80000E-01 

-2.69206E-11 

-2.6241SE-12 

2.70482E-11 

-174.43 

19 

1  .90000E-01 

-1 .66766E-11 

-3.78863E-U 

4.13944E-11 

-113.76 

20 

2.000001-01 

2.25984E-U 

-3.42549E-U 

4.10376E-U 

-36.39 

21 

2.100001-01 

2.06S2SE-11 

1 .01549E-12 

2.07O74E-11 

2.81 

22 

2.20000E-01 

3.40701E-12 

-9.18170E-12 

7.79344E-12 

-69 . 64 

23 

2,300001-01 

2.9617SE-U 

-2.41783E-11 

3.82333E-U 

-39 . 23 

24 

2.40C00E-01 

4.61940E-01 

-1.91342E-01 

5.O0OOOE-O1 

-22.30 

25 

2.50000E-01 

-I.77373E-11 

1 .45107E-11 

2.29323E-U 

140.73 

26 

2.60000E-01 

-6.13491E-12 

a.93933E-12 

1  . 08333E-1! 

124.35 

27 

2.70000E-MI1 

-2.23421E-1I 

1.82666E-11 

2.88589E-11 

140.73 

28 

2.80000E-01 

-3.96842E-1I 

-1.06410E-I1 

4.  10841E-11 

-164.99 

29 

2.90000E»01 

-3.60635E-12 

-3.S7762E-11 

3.97199E-I1 

-102.31 

30 

3. 00O0OE-01 

2.32041E-11 

-1. 73751 B— 11 

3. 96128E-11 

-34.38 

31 

3.10000E»01 

9.07410E-12 

6.263S9E-I2 

1 .10271E-11 

34.62 

32 

3.20000E-01 

-4.21S73E-12 

-1 .26363E-11 

1 .33409E-11 

-108.43 

33 

3.300001-01 

1 .07805E-11 

-1.16101E-U 

1  .58434E-U 

-47.12 

34 

3.40000E-01 

1.90U1E-11 

-1  .05227E-11 

2.17290E-1I 

-29.96 

35 

3.30000E-01 

2.63S12E-11 

9.26184E-12 

2.79313E-U 

19.37 

36 

3.60000E-01 

9.84364E-I3 

1.73739E-11 

1.74038E-11 

86.76 

37 

3.70O0OE-0I 

8.84819E-12 

-2.0*727E-12 

9.09335E-12 

-13.33 

38 

3.80000E-01 

1 . 61 156E-1 : 

1. 60721 E- 11 

2.27601E-11 

44.  92 

39 

3.90000E-01 

-2.11S71E-12 

2.23632E-11 

2.24634E-I1 

93.41 

40 

4.  300001-01 

S.00OOOE-U 

3.75000E-11 

6.2SOO0E-11 

36.87 

41 

4. 1 OOOOE-O 1 

-1  .72640E-12 

4.13274E-12 

4.47884E-12 

112.67 

42 

4.23000E-01 

-9.2**49E-12 

1.24525E-11 

1.57023E-11 

126.32 

43 

4.30000E+01 

-1.44773E-11 

-3.17347E-12 

1 .S5624E-11 

-160.58 

44 

4.40000S-01 

-1.22803E-12 

3.33848E-12 

3.55719E-12 

110.20 

43 

4.53000E-01 

-7.0O51SE-12 

1  ,22428E-t3 

7.0O622E-12 

177.00 

46 

4.6000QE-01 

-3.33S90E-12 

1.34160E-I2 

3.78467E-12 

159.24 

47 

4.70000E-01 

-1.36030E-11 
7.47B3SE-12 

-2.247S1E-12 

1.37874E-11 

-170.62 

48 

4.30000E-01 

-9.49914E-12 

I .22474E-11 

-52.37 

49 

4.90000E-01 

-6.34820E-12 

9.04331E-12 

1  .11633E-U 

125.91 

50 

s.oooooe-o: 

-9.90145E-12 

-3.83301E-12 

1 .32638E-11 

-138.26 

51 

5.10000E-01 

1.7B326E-I3 

-9.34734E-12 

9.34907E-12 

-88.91 

52 

3.20000E-01 

1  .07659E-11 

-6.07360E-12 

1  .  2361 OE- 11 

-29.43 

53 

5.33000E-01 

t.H360E-U 

1 .16109E-11 

1    60880E-11 

46.20 

54 

5.40000E-01 

-4.07135E-13 

1.28334E-11 

1 .28398E-11 

91.82 

55 

5.30000E*01 

-1.81793E-U 

1 .80911E-11 

2.56471E-11 

135.14 

56 

5.60000E-01 

-5.11848E-U 

-1 .29436E-11 

5.27960E-11 

-165.81 

57     - 

3.7O3OOE*01 

-1 .4S437E-11 

-2.07O93E-U 

2.53072E-11 

-125.08 

58 

3.80000E-01 

-3.31462E-12 

-3.33420E-11 

3 . 38649E-I 1 

-99 .  89 

59 

5.90000E*01 

2.48705E-I1 

-1 .98252E-1! 

3. I8053E-11 

-38.56 

60 

6.00000E*01 

8.60474E-12 

3.589UE-12 

1  .02606E-11 

33.01 

61 

6.10000E-01 

1  . 19537E-11 

-5.27144E-12 

1.30690E-U 

-23.79 

62 

6.20000E-01 

I.12363E-11 

-3. 16946E-13 

1  .  12409E-H 

-1  .62 

63 

6.30000E*01 

1.55675E-12 

-7.03644E-12 

7.20639E-12 

-77 .  32 

Fig. 4. 6.  Printout  of  FFT  data  for  X. 


40 


FREQUENCY  DOMAIN  DATA 


FREQUENCY  UINDOU.O  TO  6. 40000E-01 IHi 1 
FREQUENCY  INTERVAL-  1 . OOOOOE-OOtHll 


COEFF . 

FREQUENCTIHI) 

REAL 

IMAC 

MACNITUDE 

PHASE  I  DEC  1 

DC  TERM 

0. OOOOOE-00 

-1.80469E-11 

MAX  FREQ. 

6.40000E-01 

2.91406E-U 

I 

1 . OOOOOE-00 

3.68841E-13 

-2.47474E-11 

2.47704E-U 

-89.13 

2 

2. OOOOOE-00 

2.24841E-11 

4.97397E-12 

2.30296E-11 

12.47 

3 

3. OOOOOE-00 

-1.50707E-11 

3.07509E-11 

3.42453E-U 

116.11 

4 

4.00000E*00 

-4.37S83E-U 

-1.18837E-11 

4.72763E-11 

-145. 44 

5 

3. 30000E-00 

2.94318E-11 

-3.23S06E-11 

4.11319E-11 

-61.22 

6 

6. OOOOOE-00 

7.92958E-12 

-1 .21028E-11 

1  .44491E-11 

-34 .  77 

7 

7. OOOOOE-00 

4.00o59E-l! 

-5.43444E-11 

4.91538E-11 

-54.39 

3 

3. OOOOOE-00 

1 .04142E-10 

3.59250E-13 

1 .04142E-10 

.20 

9 

9, OOOOOE-00 

7.4947SE-I1 

4.47302E-U 

8.72807E-U 

30.33 

11 

1 . OOOOOE-01 

4.23474E-11 

1.207O4E-I0 

1 .33947E-10 

62.61 

11 

1 . 1OO0OE-01 

-4.44133E-U 

9.20721E-11 

1.13S23E-10 

123.90 

12 

1 .20000E-01 

-4.87474E-12 

2.27892E-U 

2.3803SE-11 

104.79 

13 

1  . 30000E-01 

-7.84951E-U 

7.32328E-1I 

1 . 07332E-10 

t34.99 

M 

1 .40000E-01 

-8.49932E-11 

3.97237E-12 

B.32Q47E-11 

173.98 

11 

1 . 50000E-01 

-4.S3150E-U 

8.2S440E-12 

4.40404E-11 

149.48 

It 

1  60000E-01 

-3.44826E-11 

-3.54039E-U 

4.93430E-11 

-134.08 

17 

1 .70000E-01 

-2.99080E-U 

9.6O078E-12 

3. 14U2E-11 

162.20 

19 

1  30000E-01 

-1 .24332E-10 

-4.93340E-11 

1 .33937E-10 

-158.23 

1? 

: .90000E+01 

8.22449E-12 

-1 .5423BE-10 

1.344S4E-10 

-84.99 

20 

2. 0O0OOE-01 

1 .10939E-10 

-4.29784E-11 

1  .18973E-10 

-21. 18 

21 

2.10000E-01 

1.49743E-U 

3.70855E-11 

3.90173E-11 

73.30 

22 

2.20000E-01 

-3.94353E-11 

-4.07295E-11 

3.66923E-11 

-134. 08 

23 

2.30000E-01 

9.42521E-11 

-8.44410E-11 

1 .28024E-10 

-42.39 

24 

2.40000E-01 

1.77402E«00 

9.23497E-01 

2. OOOOOE-00 

27.30 

25 

2.50000E-01 

-2.22918E-11 

-3.72059E-U 

4.33729E-11 

-120,93 

26 

2.6OOO0E-O1 

-1  .62640 E-  12 

7.01844E-I2 

7.20443E-12 

103.03 

27 

2.70000E-01 

-8.92625E-I1 

2.97424E-11 

9.40933E-11 

141.54 

29 

2.30000E-01 

-8.51247E-U 

-1.24549E-10 

1 .50860E-10 

-124.33 

2? 

2.90000E-01 

9.52150E-11 

-1 .30260E-10 

1.41349E-10 

-33 . 83 

30 

3. OOOOOE-01 

1.014S2E-10 

9.83270E-12 

1 .01937E-10 

3.33 

31 

3. iooooE-01 

4.601S4E-I1 

4.54929E-12 

4.62397E-1T 

5.43 

32 

3.20000E-01 

3.67188E-U 

-1.84375E-U 

4. 10B73E-11 

-26.44 

33 

3.30000E-01 

7.54804E-11 

2.31134E-11 

7.9S493E-U 

18.40 

34 

3.40O00E-01 

3.43407E-I1 

5.40154E-U 

4.31023E-11 

34.  07 

35 

3.30000E-01 

1.24373E-11 

4.89966E-11 

5.0S303E-11 

75.74 

36 

3.40000E-O1 

-1.IHM-II 

5.19379E-11 

4.O50O7E-11 

120.84 

37 

3.70000E-01 

-2.04997E-12 

-t  .08416E-U 

1.10337E-11 

-100.71 

38 

3.80000E-01 

3.43730E-11 

4.25438E-U 

5.41031E-U 

49.32 

39 

3.90000E-01 

-4.68920E-11 

5.23645E-U 

7.02915E-11 

131.84 

40 

4, OOOOOE-01 

2.10938E-10 

-1  33123E-10 

2.406S7E-IO 

-33 .  98 

41 

4. 10QO0E-01 

3.S5653E-11 

-6.04535E-13 

3. 3371  Of -11 

-.97 

42 

4.2OOO0E-01 

-2.76S09E-11 

4.72S34E-11 

S.47731E-11 

120.32 

43 

4.3000QE-Q1 

-1.37932E-U 

-5.S7581E-11 

4.03434E-U 

-103.04 

44 

4.40OO0E-O1 

4.4181BE-11 

3.67949E-11 

3.74969E-11 

39.79 

45 

4.500O0E-01 

-4.83474E-U 

2.17144E-I1 

5.22008E-11 

133.91 

46 

4.40000E-O1 

-4.44132E-12 

-1  .42780E-U 

1 .30194E-11 

-108. 08 

47 

4.70000E-01 

-1  .65805E-U 

3.18693E-12 

1 .48840E-11 

169.12 

48 

4.3O0O0E-01 

3.46701E-11 

-3.72991E-12 

3.63392E-11 

-s.ai 

49 

4.9O000E-01 

-3.34344E-11 

6.39539E-U 

7.22617E-11 

117.74 

50 

5. OOOOOE-01 

-4.44793E-11 

-8.30738E-12 

4.323S4E-U 

-169. 17 

SI 

S.1000OE-O1 

-1.712S3E-U 

-8.58489E-12 

1  .91604E-U 

-133.37 

52 

5.30000E-01 

2.S7403E-11 

-9.45244E-12 

2.74211E-U 

-20 .  14 

53 

S.300O0E-01 

8.41300E-12 

4.17708E-11 

4.24096E-11 

78.61 

54 

5.40O00E-01 

-3.49198E-11 

2.35754E-11 

4.31214E-11 

140.71 

55 

S.SO000E-01 

-5.68347E-11 

3.41459E-11 

4.43134E-U 

14B.99 

56 

5.60O00E-OI 

-1  .24454E-10 

-5.27437E-11 

1 .33178E-10 

-137. 02 

57 

5.70000E+01 

-I .30833E-10 

-9.04259E-U 

1.39041E-10 

-145.35 

58 

5.80000E-01 

1 .B4022E-11 

-1  .88539E-10 

1 . 39433E-10 

-84.43 

59 

5.90000E+01 

1 .12401E-10 

-2.53973E-11 

1.13234E-10 

-12.73 
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4.4   PRACTICAL  ASPECTS  AND  WORKING  OF  THE  METHOD 

In  order  to  demonstrate  the  working  of  this  method  in  a 
practical  situation,  use  was  made  of  the  system  consisting  of  the 
North  Star  microcomputer,  the  MC  6802  interface  and  the  HP 
Multiprogrammer  discussed  in  chapter  3.  The  same  system  was  used 
to  obtain  samples  of  x(t)  and  y(t)  generated  using  the  function 
generator  and  the  RC  circuit.  These  samples  were  then  keyed  into 
the  HP  9845B  computer  to  form  the  data  on  which  the  FFT  was  to  be 
run  using  the  AC  waveform  analysis  package. 

The  package  is  a  standard  one  available  on  the  HP  9845B.  It 
has  two  modes  of  data  input  -  manual  and  from  disc  storage.  The 
simulation  experiment  was  carried  out  by  storing  128  samples  of 
data  of  x(t)  and  y(t)  on  disc  and  then  made  use  of  by  the  FFT 
program  in  the  package.  However,  when  the  actual  data  samples 
were  taken  using  the  NS  system,  the  manual  mode  was  used  and  data 
was  input  to  the  HP  9845B  sample  by  sample  in  this  fashion  and 
finally  stored  on  disc. 

Certain  aspects  had  to  be  given  due  importance  and 
consideration.  The  sampling  frequency  is  9130  Hz.  and  just  128 
samples  were  used  which  means  that  the  frequency  index  "n"  in  the 
FFT  components  represents  a '  frequency  separation  of 
fn= 1/(128* [1/9130])=  71.33  Hz.  Thus  if  the  input  signal  frequency 
was  to  figure  in  the  plots  then  some  integral  multiple  of  fn  had 
to  be  chosen  as  the  value  for  the  signal  frequency.  Hence  utmost 
care  was  taken  to  accurately  adjust  the  signal  frequency  to  this 
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value  before  sampling  was  done. 

The  results  are  shown  in  Appendix  5.  The  phase  difference  by 
the  FFT  method  was  compared  with  the  actual  phase  difference 
obtained  by  the  formula  for  the  phase  shift  of  the  RC  circuit  at 
that  frequency. 

These  plots  can  be  compared  with  those  obtained  by 
simulation.  It  can  be  seen  in  the  amplitude  plots  for  the 
simulation  case  there  is  a  distinct  peak  at  the  signal  frequency 
and  the  amplitude  is  zero  at  all  points  but  this  is  not  the  case 
in  the  plots  obtained  by  using  actual  data  samples.  These  plots 
in  fact  show  that  the  sampled  signal  contains  a  small  amount  of 
other  frequencies  as  well.  These  effects  are  considered  later  in 
this  chapter.  It  can  also  be  seen  that  the  error  in  the  case  of 
simulation  was  zero  but  in  this  case  there  is  a  finite  phase 
error  (  obtained  as  the  difference  between  the  obtained  phase  and 
the  actual  phase  of  RC  circuit).  This  may  be  attributed  to  the 
following  factors. 

1.  Finite  precision  of  ADC  (12  bits) 

2.  Just  128  samples  were  used  for  the  FFT.  Since  the  sampling 
frequency  is  9130  Hz.  at  least  as  many  samples  would  have  to  be 
used  for  a  very  good  frequency  resolution  which  would  minimise 
phase  errors. 

3.  Presence  of  noise  in  the  signals. 

4.  Inaccuracies  in  the  trigger  reference  and  the  characteristics 
of  the  function  generator. 

However   in  spite  of  the  above  drawbacks  the  error  was  only 
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0.99  degrees  as  can  be  seen  from  the  plots  and  improvements  in 
the  above  stated  factors  would  minimise  this  error  further. 

The  FFT  is  a  very  powerful  tool  for  analysis  in  the 
frequency  domain.  It  is  capable  of  a  high  resolution  in  frequency 
provided  enough  samples  are  used.  As  seen  from  the  plots  one  can 
clearly  see  the  noise  that  exists,  the  presence  of  harmonics  (due 
to  harmonic  distortion  of  signal  from  generator  etc.)  and  also 
any  dc  offsets  that  may  be  present.  Thus  apart  from  finding  the 
phase  the  amplitude  plots  provide  extra  information  as  to  the 
performance  of  the  ADC,  noise  characteristics  etc. 

In  the  earlier  discussion,  it  was  mentioned  that  for  a 
resolution  of  1  Hz.  the  sampling  frequency  must  be  equal  to  the 
number  of  data  points  or  samples.  Thus  in  a  practical  system,  the 
sampling  frequency  must  be  made  variable  to  cover  the  entire 
frequency  range.  For  example,  the  sampling  frequency  could  be  256 
Hz.  to  cover  frequencies  from  1  Hz.  to  128  Hz.  and  a  128  point 
FFT  is  required.  Similarly  the  value  of  sampling  frequency  could 
be  stepped  up  to  say  4096  Hz.  and  4096  points  are  necessary.  This 
means  a  memory  of  4K  would  be  required  in  hardware  and  the 
maximum  signal  frequency  that  can  be  analysed  using  the  FFT  would 
be  2048  Hz.  If  a  memory  of  16K  could  be  afforded  then  fs  would  be 
approximately  16  Khz.  and  the  maximum  signal  frequency  would  then 
be  8  Khz.  If  such  huge  memory  requirements  cannot  be  met  and 
assuming  a  4K  memory  is  available  then  the  frequency  range  is 
limited  to  2048  Hz.(assuming  a  resolution  of  1  Hz.  is  necessary). 
In  order  to  effectively  increase  the  range  the  following  approach 
is  recommended. (fig.4.7). 
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In   this   schematic,    the   signals  are  heterodyned  with  a 
oscillator  to  provide  an  intermediate  frequency  of  <=2048  Hz.  and 
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Fig. 4. 7.  Scheme  for  increasing  effective  frequency  range. 
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the  resulting  output  is  fed  to  a  low  pass  filter.  Once  this 
frequency  translation  takes  place  the  output  of  the  LPF  still 
contains  all  the  phase  information  of  either  x(t)  or  y(t)  and 
hence  can  be  sampled  and  processed  by  the  FFT  algorithm.  The 
mixer  may  be  realised  using  an  analog  multiplier. 

With  the  existing  system  the  sampling  frequency  is  fixed  at 
9130  Hz.  and  yet  can  be  used  to  obtain  the  phase  response  of  a  RC 
filter  from  0  Hz.  to  say  9130/2  =4500  Hz.  The  only  restraint  is 
that  the  signal  frequencies  must  be  raised  in  steps  from  0  Hz.  as 
integral  multiples  of  fn=l/NT  Hz.  Here  fR  =  71.33  Hz.  and  the 
phase  response  of  the     filter  can  be  determined. 

It  may  be  noted  in  passing  that  FFT  chips  are  available  to 
accomplish  a  very  fast  Fourier  transform  of  the  data  samples. 
Also  this  may  be  achieved  by  an  assembly  routine.  These  may  also 
be  considered  as  possible  methods  in  a  practical  system. 

Thus  the  FFT  method  is  another  practical  method  that  can  be 
used  in  the  determination  of  phase.  As  mentioned  earlier  a  lot  of 
extra  information  is  also  provided  by  this  method  and  may  be  used 
to  advantage. 

4.5  THE  CORRELATION  TECHNIQUE 

Efforts  were  also  directed  to  obtain  a  suitable  method  for 
phase  measurement  using  Correlation  techniques.  We  first  consider 
an  attempt  to  use  cross  correlation  in  attaining  our  goal.  By 
using    the  FFT  along  with  this  method  a  suitable  method  for  phase 
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determination  is  described. 

Theory  of  Correlation: 

Consider  two  sequences  x(k)=Vxsin(w  kT)  and  y(k)=Vysin(w  kT+P). 

The  cross  correlation  sequence  is  given  by  the  following: 

N-l 
C(mT)=[l/N]  £  x(k)*y(k+m)   0<=m<=N-l; 

and  with  m=0  we  have 

N-l 
C(0)         =       (V,V,/N)     £    sin  (wokT)*sin(w0  kT+P) 

N-l  *° 

=   (VxVy/2N)  £  [cosP  -  cos (2w0 kT+P)] 

N-l 
It  can  be  shown  that   £.  cos(2w0kT+P)  will  go  to  zero  if  N  is 

chosen  to  cover  one  full  period  or   integral  multiple  of  the 

period  of  x(k)  or  y(k) . 

Hence  C(0)   =(V<Vy/2)     cosP     or  cosP  =  2C(0)/  V»Vy 

If  V  and  V  are  known  accurately  along  with  C(0)   then  the  phase  P 

can  be  determined  using  the  above  relation. 

Attempts  were  made  to  utilise  this  relation  and  simulation 
was  done  on  the  HP  9845B.  After  the  generation  of  x(k)  and  y(k) 
the  two  sequences  were  multiplied  to  obtain  the  cross  correlation 
sequence  C  .  This  was  then  averaged  over  a  period  of  x(k)  and 
C(0)  was  computed.  Some  unreliable  results  were  obtained  (for  eg. 
value  of  cosP>l)  and  the  results  were  inconclusive  and  this 
approach  was  abandoned. 

A  slightly  different  approach  produced  encouraging  results. 
Consider  the  signals  x(t)=Vxsin(w0t)  and  y(t)=V  sin(w0t+P) . 
The  cross  correlation  between  x(t)  and  y(t)   is  x(t)*y(t) 
C,/t)=  V„Vy  sin(w0t)sin(w0t+P) 
=  [V„V2]    <cosP  "  cos(2w0t+P)) 
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=  [V»V),/2]cosP    -   [V,V;/2]cos(2w0t+P) 

In  the  frequency  domain  the  first  term  represents  a  dc  term 
independent  of  frequency  and  the  other  term  is  a  harmonic(2w0) 
term.  Now  if  an  FFT  is  run  on  the  cross  correlation  Cy(t)  then 
the  amplitude  spectrum  will  show  a  dc  term  which  is  equal  to  that 
given  by[V,  V,cosP]/2.  An  FFT  of  x(k)  and  y(k)  in  turn  would 
enable  Vx,  Vy,  to  be  determined  from  the  amplitude  plots. 
Thus  the  dc  term  ie.  C(0),Vx,Vy,  enable  cosP  and  hence  P  to  be 
determined. 

4.6  SIMULATION  AND  VERIFICATION 

Using  the  same  example  as  was  used  for  the  verification  of 
FFT  by  simulation  the  following  were  generated. 
x{k)=0.5  sin(2A*24*k*T) 
y(k)=2     sin(2**24*k*T+P) 

128  samples  of  x(k)  and  y(k)  were  used  to  generate  the  cross 
correlated  sequence.  An  FFT  of  this  sequence  was  run  and  the 
amplitude  plot  is  shown  in  fig.4.8.  From  the  printout  (Appendix 
6)  the  DC  term  is  found  (3.21394E-1)  which  gives  C(0).  The  values 
of  V  and  V  were  obtained  from  the  amplitude  plots  (0.5  and  2). 
Thus  cosP=2*C(0)/[0.5*2]  and  from  this  P  was  calculated  and  found 
to  be  50  degrees  and  this  was  exactly  the  phase  difference 
between  the  two  signals  x(t)  and  y(t).  Thus  the  algorithm  was 
verified. 

4.7  PRACTICAL  ASPECTS 


All  the  practical  aspects  discussed  for  the     FFT  method 
apply  here.  Other  considerations  are  the  presence  of  harmonics  in 
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Fig. 4. 8.  Amplitude  spectrum  of  the  correlated  sequence. 
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the  signals  x(t)  and  y(t)  due  to  the  harmonic  distortion  present 
in  the  function  generator  output.  The  presence  of  harmonics  will 
cause  errors  in  the  value  of  the  DC  term  C(0).  Also  additional 
frequencies  will  be  introduced  as  a  result  of  the  mixing  action 
produced  by  the  cross  correlation  (x(t)*y(t)).  In  order  to 
overcome  this  problem  filtering  has  to  be  done.  Use  of  analog 
filters  is  not  practical  and  hence  the  use  of  digital  filtering 
is  recommended. 

A  simple  low  pass  digital  filter  can  be  easily  implemented. 
The  digital  filters  have  many  advantages  such  as  flexibility, 
reliability  and  ease  with  which  response  can  be  changed.  Thus  a 
LP  digital  filter  can  get  rid  of  the  harmonic  content  in  x(k) 
and  y(k)  after  which  the  correlation  sequence  may  be  formed  and 
analysed  using  the  FFT. 
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0&ETEE  5_ 
NUMERICAL  METHODS 

5.1  INTRODUCTION 

Numerical  and  statistical  methods  may  also  be  used  in  the 
determination  of  phase.  This  method  relies  to  a  great  extent  on 
an  accurate  determination  of  zero  crossings  of  the  signals.  This 
method  outlines  how  the  determination  of  zero  crossings  of  the 
signals  can  be  used  to  calculate  the  phase  difference  .  Numerical 
methods  like  curve  fitting  and  Newton's  method  for  finding  the 
root  of  a  polynomial  are  described. 

5.2  THEORY  OF  THE  METHOD 

The  method  uses  the  fact  that  the  accurate  determination  of 
the  zero  crossings  of  the  two  signals  enables  the  determination 
of  phase  difference.  Consider  the  situation  shown  in  fig.  5.1. 
Let  us  focus  our  attention  to  the  zero  crossing  zone;  we  shall 
consider  the  positive  to  negative  zero  crossing  for  the  sake  of 
uniformity  for  both  signals.  Consider  the  samples  being  taken  at 
times  tp  t2»  t«i  ti+l,"'-n  sPaceo'  by  an  amount  T  where  T  is  the 
sampling  period  and  T=l/fg  where  fs  is  the  sampling  frequency. 

Consider  a  straight  line  joining  v(t;)  and  vtti^j.)  making 
intercepts  x  and  y  on  the  reference  axis  as  shown.   It  is  required 
to  obtain  the  zero  crossing  time  t  .  Here  x+y=T,  and  using  the 
property  of  similiar  triangles  it  can  be  shown  that  : 
tx  =  t|  +  T[l/(l+y/x)]  =  tj  +  T[Vd+v(trt,)/v(ti)] 

If  v(t;+,  )  and  v(t,-)  are  known  then  tx  can  be  evaluated. 
Similiarly  ty  the  time  of  zero  crossing  for  the  second  signal 
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Fig. 5.1.  Use  of  zero  crossings  of  signals  to  find  phase. 
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y(t)    can   be   determined.    Once   tx  and   tv   are   known   then   the 
difference  of  time  between  zero  crossings  can  be  determined  ie. 


At  =  t. 


Once  At  is  known  the  phase  difference  can  be 


calculated  as  follows: 


Fig. 5. 2.  Phase  determination  for  a  periodic  sine  wave. 
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Since  a  sine  wave  is  periodic  the  starting  point  is 
relatively  unimportant  as  can  be  seen  from  fig.5.2.  If  the  ratio 
of  fs/f0  is  say  6  then  there  would  be  6  samples  per  cycle  of  the 
waveform.  Hence  6T  would  represent  1  full  cycle  or  a  phase  of  360 
degrees.  Hence  by  computing  tx-t„=4t,  the  phase  would  be  given  by 
P=4t(fQ/fs)*360  in  degrees.  Thus  the  phase  difference  could  be 
evaluated. 

It  can  be  seen  from  fig.5.1  that  this  method  represents  a 
straight  line  fit  for  the  sine  wave  at  the  zero  crossing.  Thus 
this  is  an  approximate  fit  for  the  sinewave  at  the  zero  crossing. 
Simulation  was  done  to  study  the  working  of  this  method.  It  was 
found  that  the  method  produced  large  phase  errors  at  some  points 
possibly  due  to  certain  discontinuities.  This  may  be  due  to  the 
fact  that  the  straight  line  fit  is  an  approximate  fit  for  the 
sine  wave.  Also  further  investigation  of  the  expression  for  tx 
reveals  that  if  the  ratio  v(t  t+1  )/v(t[ )  tends  to  -1  then  the 
denominator  becomes  zero  and  this  causes  a  discontinuity.  Hence 
this  method  was  not  considered  further. 

A  slightly  modified  approach  was  then  tried  out.  Since  the 
straight  line  fit  is  approximate  a  higher  degree  polynomial  is  a 
more  accurate  representation  of  a  sine  wave.  Hence  a  third  degree 
polynomial  was  used  to  fit  a  curve  to  the  sampled  data  at  the 
zero  crossing.  The  following  procedure  outlines  the  method. 

1.  Obtain  samples  of  x(k)  and  y(k) . 

2.  Find  the  values  of  k  when  there  is  a  negative  to  positive 
transition   (zero  crossing)   of  x(k)    (this  can  be  determined  from  a 
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sign  change  of  the  successive  samples  from  -ve  to  +ve).  Repeat 
this  for  y(k) . 

3.  Fit  a  third  degree  polynomial  to  the  obtained  samples  at  the 
zero  crossings.  This  requires  two  data  points  before  and  two  data 
points  after  the  zero  crossing  point  under  consideration.  Fig.5.3 
illustrates  this. 
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p3(x)=f(ai);  0<=i<=3,  a.^  are  coefficients  of  the  polynomial. 
Fig.5.3.  Fitting  a  polynomial  at  the  zero  crossings. 
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The  third  degree  polynomial  is  fitted  at  each  zero  crossing  using 
two  data  points  to  the  left  and  to  the  right  of  the  crossing. 
Repeat  for    y(k) . 

The  polynomial  fit  is  done  using  the  divided  difference 
table  as  shown  in  fig. 5. 4. 


TJJE  DIVIDED  DIFFERENCE  X&BLE 


k 

x(k) 

1 

x(l) 

1 

x'(l)=x(2)-x(l) 

2 

x(2) 

x'(2)=x(3)-x(2) 

x"(l)=x'(2)-x'(l)/2 

x"(l)= 
x"(2)-x"(l) 

3 

x(3) 

x'(3)=x(4)-x(3) 

x"(2)=x'(3)-x'(2)/2 

4 

x(4) 

1 

PjW-xUJ+x'U)  (x-l)+x"(l)  (x-1)  (x-2)+x"(l)  (x-1)  (x-2)  (x-3) 
which  reduces  to  a0+a1x+a2x2+a3x3  where  the  coefficients  are 
given  by  a^   0<=i<=3. 


Fig. 5. 4.  The  divided  difference  table. 
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The  coefficients  of  the  polynomial  are  determined  from  the  table 
and  the  polynomial  is  given  by  the  expression  shown.  The  same 
procedure  is  adopted  for  the  signal  y(k)  to  obtain  the  polynomial 
at  the  zero  crossings. 

Once  this  is  obtained  then  the  Newton's  method  is  used  to 
obtain  the  root  of  this  polynomial,  in  other  words  the  zero 
crossing  time.  An  iterative  scheme  is  used  to  determine  the  root 
using  the  Newton's  Method. 

Thus  the  zero  crossing  of  x(k)  and  y(k)  are  obtained  and  the 
value  of  at  is  used  to  calculate  the  phase. 

5.3   SIMULATION 

Simulation  was  done  for  this  method  and  it  was  found  to  work 
reasonably  well.  The  software  used  is  shown  in  Appendix  7a  and 
the  results  in  Appendix  7b. 

Certain  very  important  results  were  deduced  from  this.  The 
ratio  fs/f0  had  to  be  fixed  at  6  for  the  best  results  (least 
phase  error).  Any  decrease  or  increase  in  this  ratio  increases 
the  phase  error.  Also  when  random  noise  is  injected  into  the 
samples  the  error  increases  as  shown  in  the  results.  A  phase  of 
360  degrees  could  be  distinguished  by  suitable  modifications  in 
the  software.  This  method  is  very  sensitive  to  noise  and  the 
ratio  of  sampling  frequency  to  the  frequency  of  the  signals  is  to 
be  fixed  at  6,  since  the  polynomial  represents  the  best  fit  when 
this  condition  is  satisfied.  However  an  averaging  effect  is 
produced  by  determining  the  phase  at  several  crossings  and  then 


57 


obtaining  the  statistical  average. 

5.4  WCRKIN3  OF  THE  METHOD 

When  actual  data  samples  were  used  (obtained  using  the 
system  discussed  in  chapter  3)  and  processed  using  this  method 
(  the  ratio  fg/f0  was  fixed  at  6)  the  error  obtained  was  around  4 
degrees.  However  by  running  a  predictor  of  the  signals  (discussed 
in  chapter  6)  which  acts  as  a  signal  enhancer  the  error  reduced 
to  around  2  degrees. 

As  a  conclusion  this  method  has  a  constraint  on  the 
frequency  ratio  being  fixed  at  6  and  also  the  method  is  sensitive 
to  noise  and  needs  a  more  accurate  representation  of  the  signals. 
Hence  the  method  may  not  be  very  practical.  It  may  be  possible  to 
use  higher  degree  polynomials  but  this  might  require  an  increase 
in  the  frequency  ratio  (>6)  thus  requiring  a  higher  sampling 
frequency  and  also  limits  the  useful  frequency  range  of  operation 
of  this  method. 
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CHAPTER  £ 
H£E  QE  PREDICTOR  MP  NOISE  CONSIDERATIONS 

6.1  INTRODUCTION 

In  this  chapter  the  effects  of  noise  on  the  signals  that  are 
used  for  sampling  are  studied.  A  suitable  method  for  combating 
noise  by  the  use  of  predictor  is  presented  and  also  other  effects 
like  finite  resolution  of  ADC,  finite  word  length  of  computer, 
truncation  etc.  are  also  considered  and  studied. 

6.2  USE  OF  A  PREDICTOR 

Sufficient  care  must  be  taken  to  ensure  that  there  is  no 
noise  in  the  original  signals.  This  can  be  done  by  the  use  of 
analog  filters  like  low  pass  filters  etc.  In  spite  of  this  there 
still  tends  to  be  some  noise  which  also  gets  "digitized"  along 
with  the  signal.  This  noise  is  assumed  have  a  zero  mean  and  to  be 
random.  It  is  now  shown  that  the  adaptive  predictor  can  be  used 
as  a  signal  enhancer  and  gets  rid  of  the  noise. 

The  schematic  of  the  predictor  is  shown  in  fig.6.1.  This 
scheme  is  similiar  to  the  adaptive  scheme  in  the  cancellation 
mode  (mentioned  in  chapter  3)  except  that  the  reference  input  to 
the  filter  W  is  a  delayed  version  of  the  input  x(k).  In  other 
words  =1  denotes  that  the  input  to  the  filter  at  time  k  is  x(k- 
1)  ie.  the  past  value.  In  vector  form  for  a  M  weight  filter  ; 
Xk'  =  [x(k-l 
g(k)  =  W'Xk 
e(k)  =  x(k)+n(k)-g(k) 

Minimizing  the  mean  square  error  e(k)  causes  g(k)  to  look 
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Fig. 6.1.  Schematic  of  the  adaptive  predictor. 
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like  x(k)  or  in  other  words  the  error  e(k)  tends  to  look  like  the 
noise  n(k).  Thus  the  filter  W  removes  as  much  of  the  correlated 
energy  from  x(k)  and  the  enhanced  output  is  available  as  g(k). 
Since  the  filter  takes  a  few  iterations  to  adapt  itself  the 
output  g(k)  must  be  taken  after  a  few  iterations  for  subsequent 
processing. 

The  predictor  can  also  be  implemented  as  an  iterative 
scheme.  Thus  the  predictor  accomplishes  enhacement  of  the  signal 
corrupted  by  noise.  An  adaptive  predictor  with  a  10  weight  filter 
was  used  to  clean  up  signals  and  the  software  is  shown  in 
Appendix  8.  The  predictor  helped  to  enhance  the  signals  used 
especially  in  the  Numerical  method  and  was  instrumental  in 
reducing  the  phase  error. 

6.3  OTHER  CONSIDERATIONS 

Certain  considerations  must  be  given  due  importance  when 
dealing  with  implementation  of  DT  (discrete  time  )  digital 
systems.  The  following  are  the  points  to  be  noted. 

1.  Input  quantization  error  is  introduced  due  to  a  finite  number 
of  bits  used  to  encode  or  digitize  the  signal  x(t)   to  form  x(k) . 

2.  Coefficient  inaccuracy  error  that  results  from  using  a  finite 
number  of  bits  to  represent  each  coefficient  wk. 

3.  Round  off  errors  that  occur  each  time  a  product  W|4*x(k-i)  is 
rounded  to  match  a  prescribed  word  length.  (Truncation  errors  are 
worse  than  roundoff  errors.) 

4.  Certain  types  of  undesired  oscillations  that  may  occur  due  to 
limited  accuracy  and  adder  overflows. 
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These  are  the  general  considerations  and  we  must  take  into 
account  only  the  pertinent  ones.  Basically  there  are  two  types  of 
arithmetic:-  fixed  point  and  floating  point.  It  can  be  shown  that 
if  a  satisfactory  representation  is  required,  the  number  of  bits 
M  required  for  an  accuracy  of  D  decimal  places  is  given  by 
M>=3.3D. 

Input  quantization  error:  This  error  occurs  when  an  ADC  with  a 
finite  resolution  is  used  to  sample  the  signal  x(t).  If  the 
encoded  output  is  of  B  bits  and  the  ADC  has  a  step  resolution  of 
q=Vref/(2B)  then  there  is  an  ambiguity  of  q/2  referred  to  as  the 
quantization  error.  Hence  we  represent  the  quantized  value 
xQ(n)=x(n)  +  e(n)  where  e(n)  represents  the  error  and  is  a  random 
quantity  which  is  uniformly  distributed  over  (-q/2,  q/2).  The 
overall  effect  of  e(n)  is  assessed  in  terms  of  its  variance  or 
power  which  is  denoted  by  d.  It  canbe  shown  that  the  value  of 
this  variance  is  equal  to  q2/12  or  k*2_2B/12.  Hence  as  B  tends 
to  0  the  value  of  the  variance  or  power  of  this  noise  tends  to 
zero.  When  xQ(n)  forms  an  input  to  a  DT  filter  then  the  output 
noise  is  modified  depending  on  the  transfer  function  H(z). 

It  can  be  shown  that  the  effects  of  product  roundoff  errors, 
coefficient  inaccuracies  and  adder  overflows  all  contribute  some 
amount  of  noise  power  to  the  overall  signal  content.  These  are 
given  due  consideration  especially  in  special  purpose  hardware 
implementation  of  DT  syatems.  However  in  a  software  implemented 
system  only  the  pertinent  errors  need  be  considered. 


62 


CHAPTER  2 
DESIGN  CONSIDERATIONS 

After  having  discussed  the  different  algorithms  it  is 
essential  to  consider  the  design  aspects  of  a  suitable 
microprocessor  based  system  which  will  accomplish  data 
acquisition  in  the  most  efficient  manner  and  also  process  the 
data  with  speed  and  accuracy.  A  discussion  of  the  existing  system 
is  followed  by  the  design  considerations  for  a  more  efficient 
system. 

In  the  existing  system  the  configuration  being  indicated  in 
fig. 7.1  it  can  be  seen  that  the  microprocessor  MC  6802 
essentially  acts  as  a  communication  link  between  the  data 
acquisition  device  HP  Multiprogrammer  (HPM)  and  the  data 
processing  device  North  Star  (NS).  This  important  role  of  the 
microprocessor  is  elaborated  later  in  this  chapter.  An  estimation 
of  the  data  acquisition  time  and  the  data  processing  time  will 
enable  us  to  make  further  improvements  in  the  existing  system.  It 
can  be  seen  that  the  data  acquisition  is  initiated  by  a  command 
from  the  NS  via  the  serial  port  to  the  ACIA  in  the  MC6802 
interface.  The  main  assembly  routine  recognises  this  command 
("A")  and  branches  to  the  ADC  subroutine.  The  procedure  is  as 
follows.  There  are  two  PIAs  (Peripheral  Interface  Adapter)  in  the 
interface  each  with  two  8  bit  ports  available. (refer  to  the 
schematic  in  Appendix  9).  The  PIA1  is  used  to  handle  all  the 
addressing  to  the  HPM  (16  bits)  and  the  PIA2  is  used  to  handle 
the  data  from  the  HPM  (12  bits)  by  utilising  the  two  ports  of  8 
bits  each.  When  the  data  is  from  the  HPM   (12  bits)   the  remaining 
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4  bits  are  masked  off. 

The  ADC  subroutine  causes  a  control  word  and  the  slot 
address  to  be  input  to  the  HPM  using  PIA1  (address  2000)  setting 
up  the  ADC  card  in  the  HPM  for  input.  A  gate  signal  is  then  given 
to  the  ADC  to  initiate  conversion  and  at  the  end  of  the 
conversion  signal  indicated  by  the  flag,  the  data  in  PIA2  data 
registers  is  stored  in  a  convenient  memory  location.  This  process 
is  continued  until  all  the  samples  (512)  are  taken  and  stored  in 
memory.  These  samples  are  then  sent  to  the  NS  via  the  serial  port 
using  the  ACIA  in  the  interface.  It  was  experimentally  found  that 
the  time  taken  for  the  whole  process  is  around  30  sees.  Hence  for 
sampling  both  signals  the  time  taken  is  60  sees.  It  has  been 
found  that  the  data  processing  time  is  around  100  sees,  and  hence 
the  total  time  for  a  phase  measurement  is  around  160  sees.  A 
method  to  decrease  the  data  acquisition  time  is  now  discussed. 

It  can  be  seen  that  the  ADC  samples  at  9130  Hz.  and  hence 
the  time  for  taking  512  samples  is  given  by  512/9130  =  0.5  sec. 
Hence  for  sampling  the  two  signals  it  is  around  1  sec.  Thus  it  is 
apparent  that  the  use  of  serial  transmission  is  the  main 
disadvantage  since  data  is  sent  in  series.  Hence  if  the  data 
acquisition  time  is  to  be  minimised  the  parallel  transfer  of  data 
is  to  be  done.  This  may  be  accomplished  by  using  the  parallel 
port  of  the  NS  as  shown  in  fig. 7. 2. 

The  data  processing  time  however  cannot  be  reduced  unless  a 
faster  processing  machine  is  used.  Attention  is  now  focussed  on 
the  data  acquisition  part  and  the  errors  generated  and  the 
improvements  suggested  for  the  same.    It  is  shown  that  if  an  ADC 

65 


rir    MUUTlPHQaRAMMCRl 


APC 
CARD     -• 


MC  6802 

KTE  RFACeL 


NORTH 
STAR 
IS    373   LATCH 


Fig. 7. 2.  Use  of  the  parallel  port  of  the  North  Star. 
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has  an  operating  range  of  ±1V  then  the  signal  amplitude  must  also 
be  of  this  range  to  enable  proper  digitizing  of  the  signal.  This 
is  shown  in  fig. 7. 3. 
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Fig. 7. 3.  Effect  of  signal  amplitude  on  "digitizing". 
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Hence  in  order  to  enable  the  "dynamic"  adjusting  of  the 
signal  amplitude  for  sampling  by  the  ADC  of  a  given  operating 
range  the  following  is  suggested. (fig. 7. 4) 
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Fig. 7. 4.  Dynamic  adjusting  of  signal  amplitude  to  suit  ADC  range. 
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The  operation  is  as  follows.  The  peak  detector  obtains  a 
measure  of  the  amplitude  of  the  incoming  signal.  This  is  then 
compared  with  a  standard  reference  Vrej  and  the  resulting  error 
voltage  is  used  to  control  the  gain  of  a  programmable  amplifier 
so  that  the  output  of  the  amplifier  is  now  suitable  for  the 
operating  range  of  the  ADC.  If  instead  of  using  a  separate 
reference  voltage  the  Vref  point  is  connected  to  VQ  and  derived 
from  the  output  then  we  have  a  feedback  control  system.  The 
programmable  amplifier  may  be  realized  by  using  a  digital 
attenuator  or  variable  resistance  as  one  of  the  resistances  in  a 
closed  loop  OPAMP  circuit  as  shown.  Thus  irrespective  of  the 
signal  amplitudes  the  output  VQ  is  always  scaled  to  suit  the  ADC 
range. 

Errors  due  to  delays  in  triggering  are  now  considered.  Since 
the  output  of  the  zero  crossing  detector  or  some  suitable  trigger 
reference  is  used  it  is  of  primary  importance  that  these  be  very 
stable  and  accurate  with  respect  to  time.  Fig.  7.5  illustrates 
the  errors  involved. 

Let  tl  be  the  time  after  the  zero  crossing  of  the  signal 
x(k)  when  the  trigger  pulse  initiates  the  sampling  process  for 
x(k)  and  that  for  y(k)  be  t2.  If  t2  is  different  from  tl  then 
this  leads  to  phase  error  given  by  p=(t2-tl)*2*n*f  where  f 
represents  the  frequency  of  the  signals.  For  a  difference  of  4us 
between  tl  and  t2  and  a  signal  frequency  of  around  1  KHz.  the 
phase  error  can  be  as  high  as  2.5  degrees.  This  difference  may 
arise  due  to  the  comparator  characteristics  in  the  ZCD  or  more 
possibly  due  to  the  noise  at  the  zero  crossings.  Hence  in  order 
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Fig. 7. 5.  Phase  errors  due  to  inaccuracies  in  trigger  reference. 
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to  overcome  this  problem  synchronous  sampling  is  suggested  using 
two  ADCs  unless  of  course  a  very  stable  and  accurate  trigger 
reference  is  provided. (fig.  7.6) 
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Fig. 7. 6.  Synchronous  sampling  using  two  ADCs. 
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Also  the  use  of  an  ADC  with  lower  resolution  (say  8  bits) 
would  increase  the  convergence  time  and  the  phase  error  to  a 
certain  extent.  In  fig.7.6  the  two  signals  x(k)  and  y(k)  are 
simultaneously  sampled  and  the  samples  are  stored  in  memory  1  and 
2  respectively  which  are  then  accessed  by  the  microprocessor.  The 
synchronous  sampling  would  reduce  the  phase  error  to  great 
extent. 

Having  considered  all  of  the  factors  discussed  a 
description  of  a  suitable  system  for  data  acquisition  and 
processing  as  follows.  The  schematic  is  shown  in  fig. 7. 7. 

The  two  signals  are  initially  conditioned.  This  includes 
reduction  of  noise  and  automatic  magnitude  scaling  of  signals  to 
suit  the  ADC  range.  The  control  circuits  and  logic  control  the 
simultaneous  sampling  of  the  ADCs  and  storing  of  these  samples 
in  memory.  The  microprocessor  then  obtains  these  samples  and 
sends  them  to  the  computer  for  processing.  Interaction  is 
provided  by  the  Input/Output  device.  The  most  important  part  of 
this  system  is  the  portion  dealing  with  the  acquisition  of 
samples  and  storing  them  in  memory.  This  portion  is  now 
discussed. 

There  are  basically  two  methods  of  controlling  the  data 
acquisition  : 

1.  By  using  the  microprocessor 

2.  Direct  Memory  Accessing  (DMA) 

The  choice  of  the  method  depends  on  the  sampling  frequency 
required.  At  nominal  values  of  sampling  frequency  the  first 
method  is  possible  but  at  higher  sampling  rates  the  second  method 
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Fig. 7. 7.  Suitable  system  for  implementing  the  proposed  methods. 
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is  desirable.   The  first  method  has  certain  limitations  which  are 
now  elaborated. 

If  it  is  assumed  that  the  microprocessor  operates  with  a 
clock  of  1  MHz.  (for  the  MC6802  )  then  the  MPU  cycle  time  is  1 
us.  This  is  the  time  taken  for  one  cycle  and  each  assembly 
instruction  takes  a  certain  number  of  MPU  cycles.  The  minimum 
software  required  to  obtain  the  data  samples  is: 
MAIN         LDAA  SAMPLE  4   CYCLES 

STAA  SOMEWHERE     4   CYCLES 

DK  4   CYCLES 

CPX     #       FIN  ADDR.      3  CYCLES 

BNE  AGAIN  4   CYCLES 

TOTAL     NUMBER  OF  CYCLES  19  CYCLES 

We  shall  assume  that  the  total  number  of  cycles  is 
approximately  20  which  means  that  it  takes  a  minimum  time  of  20us 
for  executing  these  instructions.  If  the  control  signals  (gate 
and  flag)  for  initiating  ADC  conversion  and  checking  end  of 
conversion  are  also  to  be  incorporated  in  software  then  a  minimum 
of  another  20  us  may  be  necessary  (20  cycles)  for  execution  and 
hence  the  total  execution  time  is  40  us.  To  this  must  be  added 
the  ADC  conversion  time  which  from  the  manual  is  given  to  be  50 
us.  Hence  the  overall  time  is  100  us.  Thus  the  sampling 
frequency  is  the  reciprocal  of  this  time  ie.  l/100us  =  10  KHz. 
The  value  of  sampling  frequency  experimentally  determined  is  9130 
Hz.  which  is  approximately  the  value  obtained  in  theory.  Thus 
higher  rates  of  sampling  are  not  possible  with  this  method. 

The  DMA  method  which  uses  direct  memory  accessing  is  capable 
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of  a  much  higher  sampling  frequency.  Fig. 7. 8  shows  the  schematic. 
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Fig. 7. 8.  Direct  Memory  Accessing  (DMA) , 
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This  method  is  especially  suitable  when  using  ADCs  which  can 
handle  high  sampling  rates  eg.  the  Flash  ADC.  Its  sampling  rate 
is  around  10  MHz.  However  it  is  necessary  at  the  same  time  to 
have  a  memory  with  a  short  enough  access  time  so  that  it  can  take 
in  the  data  samples  given  by  the  ADC  at  such  a  high  rate. 

In  the  schematic  shown  two  3-state  buffers  one  for  address 
and  one  for  data  are  used  to  access  the  memory  by  the  external 
devices.  Here  a  high  speed  counter  is  used  to  provide  the  address 
to  the  memory.  The  counter  is  to  be  synchronized  properly  with 
the  ADC  conversion  signals  and  the  other  buffer  (DB)  for  data 
must  also  be  simultaneously  used  to  access  the  data  bus.  Thus  the 
microprocessor  must  essentially  disconnect  itself  from  the 
address  and  data  bus.  This  may  be  done  by  HALTing  the  uP  or  by 
other  means  depending  on  the  UP  used,  (example  cycle  stealing,  TSC 
control  etc.)  .  The  DHAC  (direct  memory  access  controller)  could 
also  be  used  for  achieving  DMA. 

In  order  to  make  the  sampling  frequency  variable  a 
programmable  frequency  generator  (VCO)  may  be  used.  A  counter  may 
be  used  to  determine  the  signal  frequency  which  in  turn  can  be 
used  to  program  the  VCO  to  obtain  the  sampling  frequency.  It  is 
essential  for  proper  working  of  this  system  all  the  different 
parts  must  be  properly  synchronised  and  this  requires  suitable 
control  logic. 

However  at  high  sampling  rates  the  access  time  of  the  memory 
available  is  a  very  important  aspect  to  be  considered  and  poses 
the  greatest  problem  in  design  and  implementation.  Also  the  rise 
and  fall  times,  the  propagation  delays  and  the  logic  switching 
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delays  must  also  be  taken  into  account  for  the  system  to  operate 
well.  It  is  essential  to  study  different  types  of  memories 
available  CMOS,  Schottky  TTL,  ECL  etc.  and  compare  these  with 
respect  to  speed,  cost,  amount  of  memory  available,  power 
dissipation  etc. 

In  the  system  discussed  the  microprocessor  acts  as  a 
control/communicating  device.  The  UP  is  well  suited  for 
manipulation  at  the  bit  level  and  hence  can  do  some  scaling,  bit 
manipulation,  sign  insertion  etc.  The  data  can  then  be 
transferred  to  a  computer  for  processing  and  probably  storage  on 
disk  for  future  processing. 

As  far  as  the  processing  portion  of  the  system  is  considered 
there  are  again  two  methods  1.  using  assembly  language  2.  high 
level  language.  While  the  first  method  can  handle  only  highly 
specific  and  dedicated  algorithms  the  later  has  a  lot  of 
flexibility.  The  high  level  language  algorithms  in  the  Computer 
can  be  easily  written  and  changed  to  suit  the  type  of  data 
involved  and  the  processing  algorithm  itself  may  be  changed  for 
adifferent  problem. 

Thus  by  taking  into  account  all  these  factors  discussed  the 
final  design  may  be  implemented. 
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CHAPTER  £ 
CONCLUSIONS 

Before  we  attempt  to  conclude,  we  must  make  a  comparative 
study  of  the  four  methods  developed.  The  overall  performance  is 
not  only  dependent  on  the  algorithm  itself  but  also  on  the 
performance  of  the  system  used  to  implement  it.  However  with  the 
existing  system  the  Adaptive  method  and  the  FFT  method  seem  to 
perform  well  as  compared  to  the  other  two  methods.  A  comparison 
of  all  the  methods  with  respect  to  a  few  features  is  now  dealt 
with. 

As  far  as  simulation  is  concerned  the  Adaptive,  FFT  and 
Correlation  methods  performed  very  well.  The  numerical  method  had 
a  larger  phase  error  compared  to  the  other  three.  When  actual 
data  samples  were  used  the  adaptive  and  FFT  methods  worked  the 
best.  The  effects  of  noise  induced  large  errors  in  the  numerical 
method  and  had  to  be  combated  by  the  use  of  predictor.  Harmonic 
distortion  produced  errors  in  the  Correlation  method  but  did  not 
seem  to  affect  the  Adaptive  and  FFT  methods.  The  FFT  method  has 
an  advantage  that  it  provides  other  valuable  information 
regarding  the  signal.  The  adaptive  method  has  an  advantage  that 
it  automatically  cancels  the  noise  since  it  uses  the  cancellation 
mode.  All  the  methods  require  that  the  sampling  frequency  be 
variable  in  order  to  cover  the  entire  frequency  range  and  be 
known  accurately  and  also  require  that  the  frequency  of  the 
signals  be  known  accurately.  Having  considered  all  of  the  above 
features  it  seems  that  the  adaptive  and  FFT  methods  have  the 
greatest  advantages  and  accuracy.   Hence  the  Adaptive  and  FFT 
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methods  are  recommended  for  an  implementation  for  the  measurement 
of  phase. 
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CHAPTER  a 
OTHER  APPLICATIONS 

There  are  several  other  applications  of  signal  processing 
techniques  in  instrumentation.  The  use  of  a  suitable  system  for 
data  acquisition  enables  a  variety  of  algorithms  to  be  used  for 
processing.  The  use  of  microprocessors  in  such  applications  is 
increasing  and  some  new  dedicated  chips  called  "signal  processor" 
chips  are  now  making  their  appearance.  A  brief  discussion  of 
other  possible  applications  follows. 

Besides  the  phase  and  the  FFT,  the  Power  Density  Spectrum 
(PDS)  of  the  signal  can  also  be  obtained.  Also  digital  filters 
can  be  used  to  shape  the  signal,  alter  its  characteristics  etc. 
The  predictor  can  be  used  for  signal  enhancement  and  the 
combating  of  noise.  The  adaptive  algorithm  is  a  very  powerful 
processing  tool.  A  method  for  obtaining  the  amplitude  of  a 
sinusoidal   signal   is  suggested  as  shown  in  fig.9.1. 

It  is  required  to  find  the  amplitude  of  x(k).  Let  r(k)  be 
the  reference  input  generated  within  the  computer  at  the  same 
frequency  as  that  of  x(k).  The  reference  input  r(k)  could  also  be 
derived  by  sampling  a  1  volt  reference  signal  of  same  frequency 
as  that  of  x(k). 

As  before  when  the  error  e(k)  decreases  in  magnitude  to  a 
sufficiently  small  value  then  r(k)  is  transformed  in  amplitude 
and  phase  by  W  to  cancel  x(k).  We  obtain  the  amplitude  transfer 
function  as: 

H (w0) =w0+w1e-3w°T+ . . . +wMe_3woMr 
The  magnitude  of  the  transfer  function  is  equal  to  the  square 
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Fig.9.1J\  method  to  find  the  RMS  value  of  a  sine  wave. 
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root  of  Re(H(w0))2  +  Im(H(w0))2. 

Since  g(k)  cancels  x(k)  at  the  summing  device  the  magnitude  of 

x(k)   is  given  by  the  magnitude  of  the  transfer  function. 

This  actually  gives  the  maximum  value  of  the  signal.  So  the 
maximum  value  divided  by  a  factor  of  root  2  would  yield  the  RMS 
value  of  x(k).  So  this  is  one  possible  method  to  obtain  the  RMS 
value  of  a  sinusoidal  signal     given  its  samples. 

Other  applications  would  be  in  other  measurements  and  in 
transducer  applications.  The  output  of  a  transducer  can  be 
sampled  and  analysed  in  the  frequency  domain  to  obtain  certain 
important  characteristics  like  frequency  content  etc.  Also  after 
processing  the  sampled  signal,  the  digital  output  can  be 
reconverted  to  analog  form  by  the  use  of  EAC  and  used.  Real  time 
signal  processing  is  now  in  the  limelight. 
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OiA£IEB  10. 
Hfll  PACKAGE 
A  plot  package  was  developed  for  the  HP  9845B  which  could  be 
useful  in  plotting  data.  This  data  may  be  available  on  a  floppy 
disk  and  the  program  can  be  modified  to  access  this  data  from  the 
disk  for  subsequent  plotting  of  the  data.  The  program  may  also  be 
included  in  the  main  program  as  a  subroutine  and  accessed 
whenever  a  plot  is  required.  The  program  listing  is  given  in 
Appendix  9.  The  program  automatically  finds  the  maximum  and 
minimum  values  in  the  data  block  and  plots  the  data  and  labels 
the  axes  etc.  It  can  be  used  for  plotting  random  data  also.  It 
was  used  extensively  used  to  plot  data  obtained  from  simulation 
analysis  of  the  different  algorithms.  It  is  capable  of  further 
modifications  and  improvements  to  achieve  greater  sophistication. 
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APPENDIX  i  a. 
Flow  chart  for  BASIC  orogram.  (Pcta^t  iva  alro^ithm^ 
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APPENDIX  1  b. 
BRSIC  arograrn  listino.  (Simulation) 


10  REM  »••••••«••••••••*«••••••••••«••••«*•••*••«••*••«•*«*»•»*•••• 

20  REN  THIS  IS  A  SIMULATION  PROGRAM  UHlCH  GENERATES  TUO  SINE 

30  REM  UAVES  UITH  A  XNOUN  PHASE  DIFFERENCE.  THE  FREQUENCY  OF 

40  REM  THE  SIGNALS  AND  THE  SAMPLING  FREQUENCY  CAN  BE  INPUT  TO 

SO  REM  THE  PROGRAM  AND  IT  USES  THE  ADAPTIVE  FILTERING  SCHEME 

60  REM  (UIDROU'S  LMS  ALGORITHM)  TO  COMPUTE  THE  PHASE  DIFFERENCE. 

78  REN  THIS  VALUE  CAN  BE  COMPARED  UITH  THE  ACTUAL  PHASE  INPUT  TO 

80  REM  THE  PROGRAM. THE  AMPLITUDE  OF  THE  SIGNALS  CAN  BE  CHANGED  IF 

90  REM  NECESSARY.  ALSO  THE  AMPLITUDE  TRANSFER  FUNCTION  AT  THE 

100  REM  SIGNAL  FREQUENCY  IS  ALSO  DETERMINED. 

1 1 0  REN  *•••»•••**««»»»••*«*»•***••«»»««»»«»»»»»»«»••*•»«»•»»•»*»*•«» 

120  REM  FORM  THE  ARRAYS  FOR  THE  FILTER  INPUT  X<»>  AND  THE  FILTER 

130  REM  WEIGHTS  U<») 

1  «0  REM  •*»•••>**«*•*«•••■■•••.<••**>•••••■••*«  ••««••»•••••.,., .  HNHHHUHHHM 

130  DIN  XU6>,U(16) 

1 60  REM  *•••*•*•«»***•••••••***«•««*••«*•••**••****»•••*•«•••••«•»*•* 

170  REM    INPUT    THE   VALUES    OF    THE    SIGNAL    FREQUENCY(F) .THE    SAMPLING 

180  REM   FREQUENCY(Fl)    AND   THE   REQUIRED    PHASE   DIFFERENCE(P)    BETWEEN 

190  REM    THE   GENERATED   SIGNALS 

200  REM    — It— W—  HWHMHtOIMHW  l»mi>  mfWHWHMI  MHHH>— IHHHRMHHI 

210  INPUT  F,F1,P 

220  P0«2»PXttF 

230  T-l/Fl 

240  P0-P*PI/180 

230  REM  •■««•«■»■■■  llM»»MMWWmtMMlWttM*«tf»«WM*M»MW« 

260  REM  INITIALISE  COEFFS.  OF  THE  FILTER 

270  REM  •••■*■•••«••*•••*••*•*■*■■•••>•••*••>».«••■••«••• #*«•«•«•»»««*. 

290  REM  ••••••••••••**•«•••*.•••••*•.••••**»••...••.#.•••**,*,••••«.•• 

290  REM  FIX  THE  NUMBER  OF  UEIGHTS(COEFFICIENTS  OF  THE  FILTER) [N] 

300  REM  •»•••••••«»#*••••*•»•«•••»•»»#»•#•••••  »•••»»••»#•»»«•••  •**«•«•« 

310  N-io 

320  FOR  I-l  TO  N 

330  U(I)-0 

340  NEXT  I 

330  REM  ••••*••••*•*•«•••«••«•#*••••••••••«•«••••••«••*••••••*••••**«•• 

360  REM  FIX  THE  NUMBER  OF  ITERATIONS 

370  REM     W»tXHHHHf  *tt««*»«t««W««— tW  «f  «i»««i«0»*MM 

380  Nl-30 

390  REM     IHMHHIW Mt«*M«tW»««tt»i»t«t»i»»«t«»»f  ■ 

400  REN  FIX  THE  VALUE  OF  CONVERGENCE  PARAMETER  (VI) 

41 0  REM  #•*•*«•*•***••»•»**••«■»••««»* »•••••••*••**»*««•••« »••••••■••••• 

420  V1-.002 

430  REM     MIII»ltWI«M«WHtM»WMMIMHtHWtMttW»WWtWiW««MtMt 

440  REM  START  PROCESSING 

430  REM  M«imtit»»*i«imM««Mit»>mnn>«wM»tM*tMHw*(Mn<wttMM 

460  FOR  K-i  TO  Nl 

*80  REM  G£N  Y  THE  PRIMARY  INPUT  TO  THE  FILTER 

490  REM    WMWiHI ■Mti»<m>tM»w«n»«t«MW<mn»wt 

300  T-SIN<F0«<»T*P») 

310  REM  <MI»— WW»>WMt««tWWWWM«»MMMl*»Wlt» 

320  REM  GEN  XI  THE  REFERENCE  INPUT  TO  THE  FILTER 

330  REM  •«»»»«««•«»«»««««»»—«—»»«««««»««««■»■«»«»«««« 

540  X1-SIN(F0»K«T) 

350  REM  •---------.-----.--....--...-.--.............». 

560  REM  ESTIMATE  THE  INPUT  SIGNAL  ENERGY 

570  REM  TO  UPDATE     THE  CONVERGENCE  PARAMETER  V 

580  REN  •*•••**•••••*»••«•••»••*»••*•••••« 

590  U2»V2»<t-Vl)+Vl»Xl«Xl 

600  IF   V3<.001    THEN  630 


as 


410  V-U1/U2 

621    COTO  650 

630  U-0 

6-40  REM  *•««••••«•«••*••»«•«•«•«••*••••***••««»««•«•»•«««««*«»•««••«««««« 

610    REM  FIND  THE  N  PAST  VALUES  OF  XI  TO  FORM  X  THE  INFUT  TO  THE  FILTER 

660  REM  TO  FIND  THE  FILTER  OUTPUT  C 

670  REH  *»•••#»»«»•»•»»»••»«••••»••«••»••«•#«»■»#•»»»»••»"«»#*»»*»*»•»»• 

680  FOR  J-l  TO  N 

690  M-J-l 

700  X<I>-SIN<FOAT«<K-H>) 

710  NEXT  J 

730  G-0 

730  REH  •  •»*••■•••**»»»•»••••••••*••»-«•••••••»••«»«•■•»*•*»«•••*•«»•••« 

7«0  REH  FIND  THE  FILTER  OUTPUT  G 

760  FOR  L-l  TON 

770  C1-W<L)»X<L) 

780  S-C-Cl 

790  NEXT  L 

800  REH  ••«•«»•«»*••»••■•««•••«»«•««••••«••»•»«••»•*•••«««••*•••••«*»«•«• 

81 0  REH  COHP  ERROR  BETWEEN  THE  PRIMARY  OUTPUT  AND  THE  FILTER  OUTPUT 

820  REM  ................................................................. 

830  E-Y-C 

840  PRINT  E, 

830  REH  •«•«»•<••««•«••»•»••»•-*»•■»•••««•••.•*•••...••»••...*.».•••.•*«.. 

860  REH  UPDATE  COEFFICIENTS  OF  THE  FILTER  USIMC  THE  ERROR  E 

870  REH  •••••»••••••••«»««•»*««.*»••»•■>»**•■*«•*••-••••••»•••..•«•*««*••..». 

880  FOR  1-1  TO  N 

890  U(I)-U(I)*U»E«X<!) 

900  NEXT  I 

9t0  NEXT  K  ._ 

920~'  PRINT  ■  •  

930  PRINT  -COEFFICIENTS' 

941  FOR  Ia|  TO  N 

930  REM    PRINT    COEFFS. 

960  PRINT  U(X>, ■ 

970  NEXT    I 

980  REH     »W4— tWWMiMWMt WWWMwmw^BM,, 

990  REN  FIND  PHASE  FROH  COEFFS.  OF  THE  FILTER.  S  IS  THE  REAL  PART  OF 

1000  REM  THE  TRANSFER  FUNCTION  AND  Z  IS  THE  IMAGINARY  PART. 

1010  REM  ••""••••»•.»•»•••...-..•.•...........•.............,,.„..„.. 

1020  S-0 

1030  FOR  1-1  TO  N 

1040  S-S*U(I)«COS<F0»T«<I-1>) 

1030  NEXT  I 

1 060  Z-0 

1070  FOR  J-l  TO  N 

1080  2"Z*U<J)«SIN(F0»T»(J-D) 

1190  NEXT  J 

1100  Z— Z 

1110  R-Z/S 

1120  PRINT  -RESULT  -,R 

1131  C'ATN(R) 

1140  C-CP180/PI 

1130  REH  •"**•«*••»••»••«*•••*•••»«•««*»•*•••••»«■•••»»»..»,..».„......, 

iiH  !E  JfiUSI  I™  UALUE  0r    THE  °»T«NED  PHASE  TO  COVER  THE  PHASE  RANGE 

1170  REH  FROH  0  THRO  360  USING  THE  INFORMATION  OF  THE  SIGN  OF  S  AND  Z 

1 180  REM  ••••••••••••••••••••••••••••••#»»»##»#»#»##»#M##„#_fc«#M„-__ 

1190  IF  <3<0>  AND  (Z>0)  THEN  C1-180*C 

1200  IF  (S<0)  AND  <Z<0)  THEN  C1-180-C 

1210  IF  <S>0)  AND  (Z<0)  THEH  C1-340*C 

1220  IF  <S>0>  AND  (Z>0)  THEN  C1-C 

1240  REH  PRIHT  THE  UALUE  OF  PHASE  ON  THE  TERMINAL 

IH!  nmrr"-^'."' " * 

1270  REM    •••••••••<«»»«..«Mi...«....«.....„„.....„.....„.„. 

SI  .'£  22SJSSZ&*  *CTWL  'H"SE  T°  K  C0"""q  hith  the 

1301  REM  ••••••••»••»••••«•••«•«••«.••.........«.„„....„., 

1311  PRINT  -THE  ACTUAL  PHASE  DIFFERENCE  IS-.P  ••••••« 

1320  REM  »•••»•••*•••••••••«••#••••«##»«•**»»«««,„.„.„„„__ 

1330  REM  FIND  THE  AMPLITUDE  TRANSFER  FUNCTION 

13*0  REM  ••••••••••••«•••••»•••••••••••••»»#«»»«#,»,,„--__«„„„ 

1331  A«ol-SORCZ»Z*S»S)  ■■■■■■■■■ 
1360  PRINT  -AMPLITUDE  TFNw-  Anol 

1371  END 
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APPENDIX  1  e. 
BASIC  program  listing  (Actual  data  samples) 


£■2  BSrt  THIS  PROGRAM  DETERMINES  THE  PHASE  DIFFERENCE  BETWEEN  TWO 

33  RE."!  SIGNALS  BY  USING  THE  ADAPTIVE  FILTER  ALGORITHM. 

43  REM  ACTUAL  SAMPLES  Or  THE  TWO  SIGNALS  ARE  TAKEN  USING  THE 

53  R=H    ADC  CARD  I.N  THE  HP  MULT I  PROGRAMMER  AND  SENT  TO  THE 

63  REM  NORTH  STAR  FOR  PROCESSING  BY  THE  ALGORITHM 

73  REM  THE  SAMPLES  OF  THE  FIfWT  SIGNAL  ARE  TAKEN  INDICATED 

S3  REM  BY  :hE  FLICKERING  OF  LEDS  ON  THE  HPM  PANEL 

93  REM  Tn£  SAMPLES  ARE  THEN  SENT  TO  THE  NS  IN  APPROX  30  SEC5 

190  Ra.M  AFTER  WHICH  THE  SAMPLES  OF  THE  SECOND  SIGNAL  ARE  TAKEN 

lie  rem  indicated  BY  THE  lEDS  On  THE  HPM  PANEL. 

123  REM  AFTER  PROCESSING  BY  THE  ALGORITHM  TnE  RESULTS  ARE  DISPLAYED 

133  REM  Di\j  THE  SILENT  TERMINAL. 

143  REM  IT  IS  ASSURED  THE  FIRST  SIGNAL  LEADS  THE  SECOND  SIGNAL 

153  REM  IN  PHASE...  TmE  RESUlTS  ShOw  A  VALUE  OF  363-P  WHERE 

1S3  REM  P  IS  THE  PHASE  DIFFERENCE  IF  THE  SIGNALS  ARE  INTERCHANGED 

160  REM  FORM  THE  ARRAYS  FOR  STORING  THE  SAMPLES  OF  THE  TWO  SIGNALS 

£30  REM  Y  AND  Xl  ARE  USED  FOR  THIS 

£10  REM  X  USED  INITIALLY  FOR  STORING  SAMPLES  AND  THEN  THIS  IS 

£20  REM  S'CRsD  IN  X:,  X  IS  THEN  USED  FOR  CURRENT  PROCESSING 

£30  DIM  X<512),  Y(5l2),W!l£),Xl  (512) 

£40  REm***********************************************************^ 

£53  REM  C_ZAR  THE  ACIA  IN  INTERFACE  CMC  6832) 

£S0  REM****************************************************,*****.*.,. 

£70  FOR  1=1  TO  10 

280  iNPJTfcl.D* 

293NEX7  I 

330   REM* mutim »»»»»» » mm>m«t»»ii>imni «» i tum ,|,>t 

310  REM  SEND  COMMAND  TO  INVOKE  THE  ADC  SUBROUTINE  IN  INTERFACE 

323  P= 

330  PRINTsl, "A" 

343  REM* 

353  REM  GET  T~£.   312  SAMPLES  AND  STORE  IN  Y(») 

3S3  REM**************************************************,.,..,..,.,.,.*.,.* 

370  FOR  I»l  TO  512 

3S3  INPUT*!,  D4 

330Y(I)=VAL(D«) 

433  NEXT  I 

413  REM**** 

423  REM  CLEAR  THE  ACIA  AGAIN 

430  REM** 

443   FOR  1-1  TO  13 

453  :nput#i,ds 


************* 


********** 


*************** 


37 


4S3  NEXT  I 

47a  repi************************************************************ 

4S0  REM  ISSUE  COMMAND  TO  SAMPLE  THE  SECOND  SIGNAL 

A  90  R=tt************************************************************ 

'500  PRINT*!,  "A" 

£ 1 a  R£m.*»********************************************************** 

523  REM  SET  512  SAMPLES  OF  SECOND  SIGNAL 


R£;»l************************************************************ 
FOR  :■   1  TO  512 
INPUTtl,M 
X(I)»VAL<DS> 
NEXT  I 


533 

540 
55a 

56a 

57lo 

583 

590  REM  SCflLI.MG  THE  SAMPLES  TO  OBTAIN  THE  ACTUAL  DECIMAL  VALUES 

60a  REM  FOR  THE  FIRST  SIGNAL 

620FOR  L-l  TO  5i2 

630  IF  Y<L) <2048  THEN  GOTO  65a 

640  Y(L>-Y(L>-4096 

650  Y(L)«Y<L>*.  S05 

660  NEXT  L 

£70  REM********************************************************** 

680  REM  PRINT  ABOUT  25  SAMPLES  OF  FIRST  SIGNAL 

690  REM«»»il»»»»«»»«i'»>»»«»'"'»">'"'l»»'"''»»'l»"»"»»""»llltl>lt«»l"t'""'»»''l> 

700  FOR  L«l  TO  25 

710  !  Y(L>. 

720  NEXT  i_ 

730  PRINT"/" 

750  REM   SCAs-INS  FOR  THE  SECOND  SIGnAi. 

760  REM********************************************************** 

770  FOR   I»i  TO  512 

760  IF  X(I>  <2348  THEN  GOTO  S00 

790  X(I)=X(I)-40SS 

S00  X(I)=X(I)*.005 

610  NEXT  I 

620  RErc*********************************************************** 

333  REM  PRINT  25  SAMPLES  OF  THE  SECOND  SIGNAL 

840  REM********************************************************** 

353  FOR  I»l  TO  25 

860  !  X(I), 

673  NEXT  I 

860  Ra*;.***-********************************************************* 

690  REMSCALE  THE  AMPLITUDES  OF  THE  SIGNALS  TO  IV 

530  RE?"************************************************************ 

Sia  FOR  »■  1  TO  512 
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9£2Y(I)=.  i*YC) 

538  XI  <I)«X(D*.  i 

S40NEXT  I 

95a  I "  " 

9S0  SE>**************************«*************+******«-»****»****** 

973  SEM  SET  THE  FREQUENCY  OF  The  SIGNALS 

988  SS!^*+**********-)Hf***»***********#*****»************************ 

99«  !"FREQ", 

isaa input  f 

i  a  i  a  rb^i************************************************************ 

102a  REtf  START  PROCESSING 

1030  fiE^**+-»Ht»******+****#*##****#***#**#»##**+***********»*+******* 

1348  F0=E*(££/7)*F 

iasa  Fi»5i3a 

:3Sa  T=l/9. 1323 

10Sa  REM  FIX  THE  VALUES  FOR  THE  NUMBER  OF  WEIGHTS   FOR  THE 

10912  REM  DIGITAL  FILTER 

1100  REM  FIND  THE  RATIO  OF  SAMPLING  FREQUENCY  TO  THE  SISnAl  FREQUENCY 

ma  rem****************+**»******************************»********* 

1123  IF  Fl/F  (»15  THEN    GOSUB  1160 
U30  IFCFl/F  <»il  AND  Fl/F>15)  THEN  GOSUB  l£00 
1142  IF  (Fl/F<=9  AND  Fi/F>:i)  THEN  GOSUB  1£40 
115a  IF  <Fl/F(»5  AND  Fl/F>9)  THEN  GOSUB  1280 

H60  N=a 
Ii7«  VI".  MB 

iiaa  M=aaa 

11 3«  RETURN 
12S0  N=4 

1210  vi=. ea* 
122a  M»3aa 

123a  RETURN 
1240  N'=4 

i£sa  vi=.  aa3 

12S3  N 1*400 
l£7a  RETURN 
1283  N=3 
1290  Yl». 003 

:3aa  M=4aa 

13 10  RETURN 

13£^RZ^************************************************************ 

1330  REM  START  PROCESSING  USING  THE  DETERMINED  VAi-UES  OF 

1340  REM  THE  WEIGHTS  (IV),  CONVERGENCE  PARAMETER <V1 >,  AND 

1350  REM  THE  NUMBER  OF  ITERATIONS (Nl > 

t3S&FtZ*.******************+***+***+******************************-**+ 

1370  INITIALIZE  THE  WEIGHTS (COEFFICIENTS)  OF  THE  FILTER 
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i3&iiRz™*+**#*****+***************+*************+** 

i3Si2  "OR   I"   :   TO   16 

It j£   w  i I J    =_  B 

14:8   NEXT    l' 

142S   FOR    K=iv   TO   Nl 

l  »3v>R£w»»*»«*»»«  »  minimum II  »» nnt  »<mmhmh>  »«»»»■  »  »imn n  H  »«»»»»» II  II IMI  || »» 

1448    RE,»:    FIND    INPUT    SISMOi.    ENESSY    TO    UPDATE    ThE 

1438   REM  C0NVER5ENCE   PftRfiMETER 

l460RSi»i<MMHHMnnnt  »»«mcinm«  **********************  «  mummmiimiinimimnnm 

1474  VS»V£*(1~V1)+V1*X1 (K)*Xi (K> 

1463    1?    V2(.0Oi    TriE;\    1310 

1498   V-V1/V2 

:588   3u~0    1530 

iSlio   V=v'. 

: 5£ia$E,'lHHHMHHHHMHMHHm '«■»  » » ********  ««»»»»<  *  «  mum  »  mm  m  »»■»»»<»»»«  m  «  unit 

153aREM  BEN.  X 

1S48SSM«»«  mm  »««♦»» »»»»inm  »«*»«»»««» »>  imimimimimimmiiniiimmmniimiMim 

1558      FOR   J=i   TO   N 

:5£a      :Y'-j'-J. 

:37a  x(j)=x: <k-«> 

1=55      NEXT   J 

:59a    3=a 

:6agR£:'1»«»«*<H>»»»<M>»»»»»<MM>#»»«it»i>iiii»iiiiiiiiii»ii»»Kiiiiilllli»»i(»»ii»ii»»iiii»ii» 
ISia   REM   FIND   THE   OUTPUT   OF   TnE   FILTER   3 

lSa0REin»**«*»»»»*«HHH>»»».|MMM|.»»»»»»«.»»»»»»IMMi.<M»»»<|||<|  |m  «  »» »»,|m  »»»«*« 
■  CT-7i    ~Tw  ■  —   i   —  ~  Hi 

:3"Ci   3l=,v<i_)*X  •;_> 
1&S8   3=3+Gi 

.358   NEXT  u 

l^QfP.i^**********^***********************^**^^^^^^^^^^^^^^^^^^^^^^^ 
iSaa  RE^  CO:»iP.  ERROR  BETWEEN  PRIMARY  INPUT  &  FILTER  O'JTPiJT 
:s3»RE^******************************#«**»»*****»****#*#*»***#^^+^+ 

:7ae  e=y(k^-3 

17.'?   FOR  1  =  1   TO  N 

:  7Z2?=******-************+*****************++*******+*+++++„++++^++ 

173C-  Rtii  LPBSTE  CGEFFS. 

^wRE:';***********************^*********^^ 

17s?   w<I>"W(I>+V»S#X<l) 

;7£ei  ^ex^  : 

1778   NEXT  K 

: 7sC   !  " ,' " 

1798   :  "COEFFICIENTS" 

1838  FOR  I  =  ;   TO  N 

i3l8    REM    PRINT    C0EFF3. 

:SiC    i    w(I). 

l&Zii   ;\EXT    I 
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**■*■*  ■*•*■*■#■*-*  •*■■*•***■*  -*-.f  *■**■*"*  -R-  -***-*•*  *  *  * 


i.353    SEil   FIND    PHASE   FRCi*.   CQEFFS. 

1  SSflRSfli******  ft  *****  »*»»»*#*»#**  «»  It********************  It  H  a  IHHM'HIISrfHO 

l373  REM  S  IS  REAL  POST  OF  THE  TRANSFER  FUNCTION 
lBS®REft****************#****#************"****iMt*HHHHMHHH*******#**** 

LS9<a  s-a 

.9*3  FOR  I  "J   TO  N 

1913  S=S+w( I ) »C0S<F8#T*<i-l) ) 

i'3£3  i\EXT  I 

L93itREm************+********************************»**+-f  •»*****■**** 
:34a  REM  Z  IS  THE  IMAGIimARY  PART  DF  THE  T8AM9FER  FUNCTION 

1 95fl3£M<HHMHMHMHHMMHt»  »**»»*  »#»*  It  It  *  H  »» II  »»»  » It  It********* »  »«»*»»«  *»»««* 
196*     Z=0 

1370   FOR    J=i      TQ    :\ 

1338    Z=Z-t-*J<J )  *SI.'i  (F0*T*<  J- i)  ) 

L993  NEXT  J 

=333    Z=-Z 

;3l3   R=Z/S 

;3£2!    i     "  RESULT",     R 

29318   OATN<R3 

2SJ40   C=C*133/  (£2/7) 

2353R£i**********************************+***^*+*******-**^ 

£3&'2  RE.*!  ADJUST  Tr;E  VfiuuE  CF  TriE  OBTAINED  P.-ASE  TD 

2373  SEW  COVER  "THE  ENTIRE  RANGE  (3  THRO  363)  USING  THE 

i'3a3  REM  INFORMATION  FROM  S  A.mD  Z  5I*jCE  THE  FILTER  CAN 

i3S3  REM  CAN  DISTINGUISH  A  RAN3E  (3  THRO  93  ONLY) 

i::33  RE.'i  PRINT  THE  OBTAINED  PhASE  AFTER  ADJUSTING  D.m  ~-E  SILENT 

ill  13  REM  TERMINAL. 


J123RE.11***** 

1:133  IF  S<3  A;vD  Z>3  THEN 

1:143  IF  S<3  AND  Z  <3  THEN 

il53  IF  S>3  AND  Z  <3  THEN 

11 63  IF  S>  3  AND  Z>3  THEM 


THE  PHASE  DIFFERENCE  IS",  lfl«*C 
THE  PHASE  DIF-ERE.vCE  15". C+:e3 
7n=    PHASE  TJZfFZRS.^CZ    IS",368+C 
TmE  ?~SaE  DIFFEREivCS  _'S".C 
Sl73fiEftt»***»  »  *  »  it*************************************************** 
c:l33  RE."  PRINT  THE  ACTUAu  VALUE  GF  PHASE  FRO*  RC  CXT. 
£193  REM  HERE  THE  VALUES  GF  R  A.*:D  C  ftRE  FIXED 

£S38R£;1»»»'»*  »»»***ihh>**4hmhhh*ihhmii»»ii  »»«»»** »*#*»int**it* »**»*»****** 
££13  ! ATN(1/  <F3*:3l.  l5E3*4.  394E-9;  )*133/  (££/7) 
2££3REM**  ***+***+***+*********+*****+**********■***  >****^ 
£233  RE.1  THE  ACTUAL  PiiASE  CAN  BE  COMPARED  wITH  THE  OBTAINED  P-ASS 
£S43S£fl*******in»  »»n  »»***»***■****♦»***»«*»  «■»  tmnMHHhr***»)Hi»int**innHHi 
££53  END 
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APPENDIX  2 

Schematic  diagram  of  the  M6802  Based  Interface 


IfeH  fgjTI  ggsggssgei 


..Tftfrme 


FFFPFFPFPrFJ 
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APPENDIX   3 
6802  Monitor  Program  Listing. 
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00001 

NAM 

MICROl 

00002 

0050 

STAKPT 

EQU 

HO 

0C003 

AOOO 

PTMCi 

EQU 

1A000 

• 

0000* 

AOO* 

PTM8 

EQU 

1A004 

00005 

A001 

PTMC2 

EQU 

SA001 

OC006 

0000 

INOL 

EOU 

too 

00007 

00  60 

INOH 

EO.U 

160 

00008 

0011 

MEM 

ECU 

SU 

C0009 

0012 

MEM1 

EOU 

112 

00010 

2001 

CRA1 

EOU 

$2001 

00011 

2003 

CRB1 

EOU 

$2003 

00012 

2000 

0ORA1 

EOU 

S2000 

00013 

2002 

0ORB1 

EOU 

12002 

0  001* 

*001 

CRA2 

EOU 

$4001 

00015 

4003 

CR32 

EOU 

1*003 

00016 

4000 

00RA2 

EOU 

$4000 

00017 

4002 

DOR  82 

EQU 

$400  2 

00018 

9000 

AC  IAS 

EOU 

S8000 

00019 

3000 

AC1AC 

EOU 

ACIAS 

00020 

8001 

ACIARO 

EOU 

S8001 

00021 

8001 

AC  I  AT 

EQU 

ACIARO 

00022 

0010 

CP 

EQU 

110 

00023 

OOOE 

VALL 

EQU 

$0E 

0002* 

OOOF 

VALH 

EQU 

JOF 

0C02S 

OOOC 

TMPL 

EQU 

SOC 

00026 

0000 

TMPH 

EOU 

SOD 

0  0027 

0008 

KNT 

EQU 

SOB 

00028 

OOOA 

MARK 

EQU 

$0A 

00029 

0006 

LS03 

EQU 

S06 

00030 

0007 

LS02 

EQU 

107 

00031 

0008 

LS01 

EQU 

$08 

00032 

0009 

LSO 

EOU 

109 

00033 

0001 

SUSHI 

EOU 

$01 

0003* 

0000 

SUSLO 

EQU 

$00 

00035 

0003 

MINHI 

EQU 

$03 

00036 

0002 

NINLO 

EQU 

$02 

00037 

000* 

MEG 

EQU 

10* 

0C038 

0005 

CARRY 

EQU 

$05 

0C039 

001* 

TEMP 

EQU 

«1* 

000*0 

0016 

TYPE 

EQU 

$16 

000*1 

0017 

TVALL 

EQU 

117 

000*2 

0018 

TVALH 

ECU 

$18 

000*3 

0019 

FFLAG 

EQU 

119 

000** 

EOOO 

ORG 

SEOOO 

000*5 

EOOC 
E001 
E002 
£003 
E00< 
EOOi 

00 

5F 
00 
17 
00 
05 

Fca 

100.15F.iOO 

$17,100,105,100,102 
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00046 

00047 

00048 

00049 

OOOSO 

00051 

00052 

0  00  53 

00054 

00055 

00056 

OC0S7 

00058 

00059 

00060 

00061 

00062 

00063 

00064 

00065 

00066 

00067 

OC068 

00069 

00070 

00071 

0C072 

00073 

00074 

00075 

00076 

00077 

00078 

00079 

00080 

00081 

00082 

O0O83 

00084 

O0085 

00086 

00087 

00088 

00089 

00090 

00091 

0C092 

00093 

00094 

00095 

OC096 

O0097 

0C098 


E006   00 
E007  02 


PIAiACJA,   AND   PTH    INITIALIZATION 


E008 

E008 

EOOO 

E010 

E013 

E015 

E018 

E01B 

E010 

E020 

E022 

E025 

E027 

E02A 

E020 

E030 

E033 

E035 

E03S 

E03A 

E030 

E03F 

E042 

E044 

E047 

E04A 

E04C 

E040 

E04E 

E04F 

E051 

50  53 

E055 

E057 

6059 

E05B 

E050 

E05F 

E062 

E064 

E067 

E06A 

E06C 


SE    0050 
86    OC 
B7    2001 
97    2003 

86  FF 

87  2000 
37  2002 

86  04 

87  2003 

86  3C 
B7  2001 
36  00 
B7  4001 

87  4003 
87  4000 
87  4002 

86  04 

87  4003 
86    3C 
B7    4001 

86  03 

87  8000 
36    01 
87    3000 
86    4002 
84    30 
44 


38  00 
97  12 
36  EO 
89  00 
97  11 
OE  11 
EE  00 
36  93 
87  AOOl 
36  00 
B7  AOOO 
FF    A0O4 

86  83 

87  AOOl 


LOS 
LOA  A 
STA  A 
STA  A 
LOA  A 
STA  A 
STA  A 
LOA  A 
STA  A 
LOA  A 
STA  A 
LOA  A 
STA  A 
STA  A 
STA  A 
STA  A 
LOA  A 
STA  A 
LOA  A 
STA  A 
LOA  A 
STA  A 
LOA  A 
STA  A 
LOA  A 
AND  A 
LSX  A 
LSR  A 
LSR  A 
AOO  A 
STA  A 
LOA  A 
AOC  A 
STA  A 
LOX 
LOX 
LOA  A 
STA  A 
LOA  A 
STA  A 
STX 
LOA  A 
STA  A 


ISTAKPT 

•  SOO 

CRA1 

CRB1 

•8FF 

O0RA1 

00RB1 

1104 

CRB1 

»13C 

CRA1 

•too 

CRA2 

CRB2 

00RA2 

D0RB2 

•804 

CR82 

«»3C 

CRA2 

8803 

AC  I  AC 

•  801 

AC  IAS 

O0RB2 

8S30 


•  INDL 
MEMl 

•  INOH 

•  SOO 
MEM 
MEM 
0,X 
1593 
PTMC2 

•too 

PTMC1 
PTMB 

•  833 

PTMC2 


POINT  CRA1  AT  00RA1 
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00099 
OOIOO 
00101 
00102 
00103 
00104 
00105 
00106 
00107 
001C8 
00109 
00110 
00111 
00112 
00113 
00114 
00113 
00116 
00117 
00118 
00119 
00120 
00121 
00122 
00123 
00124 
00125 
00126 
00127 
0C12S 
00129 
0C130 
00131 
00132 
00133 
0C134 
OC135 
00136 
00137 
0C138 
00139 
0C140 
00141 
00142 
00143 
00144 
00145 
00146 
00147 
00148 
00149 
00150 
00151 
00152 
00153 


E06F  80 
E072  81 
E074  26 
E076  BO 
E079  20 
E073  81 
E070  26 
E07F  80 
E082  20 
E084  81 
E086  26 
E088  80 
E088  20 
E080  81 
E08F  26 
E091  80 
E094  2  0 
E096  81 
E098  26 
E09A  30 
E090  20 
E09F  81 
E0A1  26 
E0A3  80 
E0A6  20 


MAIN  ROUTINE  TO  CHECK  FOR  VARIOUS  COMMANDS  C, ?.«.*.  .'A* 


»*****»•******«**•*•-*•»«»»»»••»«-•*.»•■»*** 


E18F 

26 

05 

E100 

F4 

3F 

05 

E200 

E8 

23 

05 

E400 

E2 

25 

OS 

E50  0 

09 

21 

05 

E60  0 

00 

41 

cc 

E19  8 
CT 


LOOKCO  JSR 

CMP  A 
Ml 
JSR 
BRA 

CKQN    CMP  A 
BNE 
JSR 
BRA 


CKLB 


SETONE 
i*26 
CKOM 
OOAT 
LOOKCO 
Ktf 
CKL8 
RONSNO 
LOOKCO 
CMP    A      IS23 
BNE  CKREL 


JSR 
BRA 


TIMER 
LOOKCO 


CMP    A      1(25 

BNE  CKVCO 


JSR 

uu 

CKVCO  CMP  A 
8NE 
JSR 
BRA 
CMP  A 
BNE 
JSR 
BRA 


CAOC 


RELAY 

LOOKCO 

•121 

CAOC 

VCO 

LOOKCO 

»S4l 

LOOKCO 

AOC 

LOOKCO 


GET   THE    NEXT   CHAR   FROM   NS    WHEN   SENT 


IF   NOT    7  THEN  CHECK    IF    IT    IS    < 

REAO      HP        ANO    SEND      RESULT   TO    N.S7AR 

CHECX    IF   CHAR   ECUALS    I 

IF    NOT    »•«    THEN    CHECK    IF    RELAY    COM    "X" 

INITIALIZE  V/F    CARD    £    TIME    IT 

CHECK    IF    IT    IS    A    RELAY    COMMAND 

IF   NOT      RELAY   THEN   CHECK    IF    IT    IS   VCO   COM 


CHECK    IF    VCO    COMMANO    IE    "    « 

IF   NOT    VCO   COMMANO   THEN   CHECK    IF    AOC   COM 


CHECK    IF    A/D   CONVERTER    COMMANO.    IE    "A" 


«M*uatMiM**nrn*Mi**n«ttwayMmiataimtlslaM1Munl)M 

E   N    0        OF        MAIN      ROUTINE 


»»•••>*»*»»•»»•***»«»•»»*«*»•«»»«»»»••»»». 


OOAT    SUBROUTINE      -   TO    ACCEPT    CHARACTERS    FROM    N.STAR    UNTIL    J 
CONVERT   TO    BINARY.    ANO    SEND   RESULT    TO    HP 


E100 
E100  BO 
El  03  96 
E105  27 
El  07  BO 
E10A    39 


E10B    COAT 

10 

03 

E16B 


ORG 
JSR 
LOA    A 
BEQ 
JSR 


XODAT      RTS 


IE100 

NSREAD 

CP 

XOOAT 

STOHP 


GET    CHARS   FROM    N.STAR    UNTIL    "J" 
IF    NOTHING   SENT    THEN    EXIT 

SEND  VALUE    IN   VALL    C   VALH   TO    HP 
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0015* 

00159 

00136 

00157 

00158 

00159 

00160 

00161 

00162 

00163 

00164 

00165 

00166 

00167 

0C168 

00169 

00170 

00171 

00172 

00173 

0017* 

00175 

00176 

00177 

00178 

00179 

00180 

00181 

0C182 

00183 

00184 

00185 

00186 

0C187 

00188 

00189 

00190 

00191 

00192 

00193 

00194 

00195 

00196 

00197 

00198 

00199 

00200 

00201 

00202 

00203 

00204 

00205 

00206 

00207 

0C208 


•*****••»»»»••*•»«•••»•»•»»•••»»» 


NSREAO   -   SUBROUTINE   TO   REAO   N.STAR    S    CONVERT    TO   BINARY 

^^rJ?UT^e,MUL   RKE,VE   ASC"    D,SlTS   "ON   THE    ™ST«   SnTIL   A- 
CONVERT      TO   16    BIT    BINARY    I    STORE   RESULT    IN    VALL    C   VALH. 


E10B 

El  00 

El  OF 

£111 

El  13 

E115 

6117 

E119 

E11A 

E11C 

EUF 

E121 

E123 

E125 

E126 

E129 

E12B 

E120 

E12F 

E130 

E132 

El  34 

E136 

E139 

E13S 

E130 

E13F 

E141 

E142 

E143 

E144 

El  45 

E147 

E149 

E14A 

E14B 

E14C 

E140 

E14E 

E150 

El  52 

E155 

E157 

E159 

E15B 

E150 

E15F 

E162 


86  00 
97  10 
97  OE 
97  OF 
97  OC 
97  00 
97  OB 
4C 
97    OA 

80  E18F 

81  24 
27  C8 
80  30 
36 

7C    0010 
20    Fl 

96  10 
27  3B 
32 

06  OB 

26  00 

97  OE 
7C    OOOB 

96  10 
91    OB 

27  28 
20    EE 
5F 
OC 
48 
59 

97  OC 

07  00 


NSREAO  LOA  A  1100 
C? 
VALL 
VALH 
TMPL 
THPH 
KNT 


STA  A 
STA   A 

STA    A 

STA   A 

STA   A 

'      STA   A 

fNC   A 
STA   A 

L00KC2    JSR 

CMP   A 
8EQ 
SUB    A 
PSH   A 
INC 
BRA 

CHECK      LCA   A 
BEQ 

NEXT    PUL  A 
LOA  B 
BNE 
STA  A 
INC 
LOA  A 
CMP  A 
BEQ 
BRA 

KNTPOS  CLR  B 


XI 0 


59 
48 
59 
OC 

9B  OC 
09  00 
7A  OOOB 
26  EB 
98  OE 
09  OF 
97  06 
07  OF 
7C  OOOA 
96  OA 


CLC 
ASL  A 
ROL  B 

STA  A 
STA  8 
ASL  A 
ROL  B 

ASL  A 
ROL  B 

CLC 
AOO  A 
AOC  B 
OEC 
BNE 
AOO  A 
AOC  B 
STA  A 
STA  B 
INC 
LOA  4 


MARK 

SETONE 

8824 

CHECK 

8*30 

CP 

L00KC2 
CP 
XNSRO 

KNT 

KNTPOS 

VALL 

KNT 

CP 

KNT 

XNSRO 

NEXT 


TMPL 

TNPH 


SET  THE  NEXT  CHAR  FROM  NS  WHEN  SENT 
IF  "J«  THEN  GET  OUT 


TMPL 

THPH 

KNT 

X10 

VALL 

VALH 

VALL 

VALH 

MARK 

MARK 
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PACE 


00209 

00210 

00211 

00212 

00213 

0021* 

00215 

00216 

00217 

00218 

00219 

00220 

00221 

00222 

00223 

00224 

0  0223 

00226 

00227 

00228 

00229 

00230 

00231 

00232 

00233 

00234 

0C235 

00236 

00237 

00238 

00239 

0C240 

00241 

00242 

00243 

00244 

00245 

00246 

00247 

00248 

00249 

00250 

00251 

00252 

00253 

0  0254 

00255 

00256 

00257 

00258 

00259 

00260 

00261 

00262 

00263 


E164  97  OS 
E166  91  10 
E168  26  C5 
E16A  39 


XNSRO 


THI 
THE 


E168 

El  60 

E170 

El  72 

E175 

E176 

El  77 

E17A 

E17C 

E17F 

El  82 

£184 

E186 

El  89 

E188 

E18E 


96    OE 
B7   2000 
96    OF 
87    2002 
02 
02 

B6   4001 
84   F7 
87   4001 
86    4002 
84   40 
27   F9    ' 

86  4001 
SA    08 

87  4001 
39 


STOHP  LDA  A 
STA  A 
LDA  A 
STA  A 
NOP 
MOP 
LDA  A 
ANO  A 
STA  A 

L00KC3  LOA  A 
ANO  A 
8EQ 
LOA  A 
ORA  A 
STA  A 

XSTOHP  RTS 


STA  A  KNT 

CMP  A  CP 

BNE  NEXT 
RTS 

*" "  ■«———.. mini 

STOHP   -   SEND    TO    HP    SUBROUTINE  ••«••«...«....„., 

.f"!??"!""5    WILL    SEND   z    BVTES   STORED    IN    VALL    £    VALH    Tn 
i*JiL»lfWT  OF  HP  -   WUI   T0   HIGHcYmTc    VAU^lWr. 


VALL 
D0RA1 
VALH 
00RB1 


SETONE 


CRA2 
#*F7 

CRA2 
00RB2 

"«C3        CHE"    "*   FL*S    ™QM    """-TI PROGRAMMER 

CRA2 

»»08 

CRA2  MHEN   FOUND,    RESET   THE    SATE 

—***—**— r ,- 


SUBROUTINE  TO  GET  I  CHAR  „N  REG  Ai  FROM  MS.  -HEN  SENT 


E13F  36  8000 
£192  47 
E193  24  FA 
E19S  86  8001 
E198  39 


GETCNE  LOA  A   ACIAS 
A$R  A 

BCC     GETONE 
LOA  A   ACIARD 

XGETON   RTS 


«tio    r^tri    CHRACTER    IS    SENT 
REAO    IT    INTO   REG    A   WHEN    SENT 


SE0NO-tSeT?STtJ:I   lifSIK    '^ll'lSr    !    S"'L«    ™ 

BE    IN   SLOT   ioroS'TSrMUL'T^ROGRA^MER00   M°    "    A"U"E°    T0 


E19B 
E19B  CE 
E19E  OF 
E1A0  BO 
E1A3  CE 
E1A6  OF 
E1A8  CE 
E1AB  BD 
ElAE  B6 
E181  A7 
E1B3  08 
E1B4  B6 
6187    A7 


'"' '"" ■ 


SOFO 

OE 

E16B 

0050 

OE 

600  0 

E16B 

4000 

00 

4002 
00 


AOC 


ANEXT 


ORG 
LDX 

STX 
JSR 
LDX 
STX 

LDX 
JSR 
LOA 
STA 
I  NX 
LOA 
STA 


JE198 

MBOFO 

VALL 

STOHP 

M0050 

VALL 

IJ6000 

STOHP 

0CRA2 

O.X 

D0R82 
0,X 


SEND   CONTROL   WORO      -11U   0000    1011    0000- 

AOC    CARO    IS    IN    SLOT   60S 
*_    MODIFICATIONS     START     HCIE. 
"e?RTHS|    liS^I    "1TH    ""    SIGNAL 
STORE    IT    AT    LOCATIONS   6000    ONWAROS 
SET    THE    UPPER    BYTE   OF   THE   SAMPLE 
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00264 

00265 

00266 

O0267 

00268 

00269 

00270 

00271 

00272 

00273 

00274 

00275 

00276 

00277 

0  0278 

00279 

00280 

002S1 

00282 

00283 

00284 

00285 

00286 

00287 

00288 

00289 

00290 

00291 

00292 

00293 

00294 

0C295 

00296 

0C29T 

0  0298 

00299 

00300 

00301 

00302 

00303 

00304 

00305 

00306 

00307 

00308 

00309 

00310 

00311 

00312 

00313 

00314 

00315 

00316 

0  0317 

00318 


6189   08  INX 

EISA    8C    6400  CPX  886400 

6180   26    EC  •      BNE  ANEXT 


E1SF   CE    6000  LOX  816000 

E1C2  A6    00        AASAIN   LOA  A      O.X 

E1C4  43  COM  A 

E1C5   97    00  STA  A      00 

E1C7   08  INX 

E1C8   A6    00  LOA  A      O.X 

EICA    43  COM  A 

E1CB    84    OF  AMO  A      HOF 

E1C0    97    01  STA  A      01 

EICF   08  INX 

ElOO   DF    13  STX  $15 

E102  CE    003F  (.OX  (800 3F 

E105    80    647E  JSR  TIMOUT 

E108   80    E224  JSR  WSENO 

El  OB    BO    6209  JSR  L0OKC1 

6106  06    15  LOX  815 

61 60    SC    6400  CPX  1*6400 

6163    26    DO  8NE  AAGAIN 

61E5   39  XAOC         RTS 

£200  ORG  1E200 


•*•***•*•*•**»*»•* 


COLLECT    512    SAMPLES    IN  ALL.      IE    IX   BYTES 


G6T    READY   TO    SEND   THE    SAMPL6S    TO   NORTH    STAR 

GET   THE    SAMPLE    STORED    IN   MEMORY 

TAKE    l'S   COMPLEMENT    OF    THE    LOWER    BYTE 

GET    THE    HIGHER    BYTE   OF    THE    SAMPLE 
MASK    OFF    THE    HIGHER    4   BITS 

SAVE    THE    SAMPLE    POINTER 

WAIT    FOR    ACIA   TO   GET    REAOY 

SENO   TH6    SAMPLE 

SEND    A    CARRIAG6    RETURN 

REPEAT    UNTIL    ALL    512    SAMPLES   ARE    SENT 

,*****  *•***••***•  *•*•*••••••••*«•**.•..»,.««».„ 


RONSNO    -   REAO   AND   SENO    TO   N.STAR    SUBROUTINE 
THIS   ROUTINE    WILL   REAO   HP    AOORESSED   PORT    ANO    SEND    THE    OJRRKP 
D6CIMAL    VALUE  -UN    ASCII)-    TO   THE   NORTH    STAR    FOLLOWED    BY   2?" 


A******************* 


110034 
TIMOUT 
I  OAT 
AC1AS 


LO0KC1 

8100 

ACIAT 


'**••»****•*•  immmmmmmmmmmm  **«*••«. 
GET    A   DELAY    COUNT 
Give    N.STAR    HANCICAP    OF    1/10    S6C 
READ    HP    ANO    SENO    TO   N.STAR 


6200   CE   0034  RONSNO    LOX 

6203   80    6476  JSR 

6206    BO    E216  JSR 

E209   B6    8000  L00KC1    LOA    A 

HOC   46  ROR   A 

6200   46  ROR    A 

E20E   24    FS  BCC 

E210   86    00  LOA    A 

E212    B7    8001  STA   A 

E215   39  XRONSN   RTS 

••■•>••■  ■•■••■••II  ■■,.. »«,»».,, •»*.«» ..... 


******mii»i»)m«nmi 


"XZm*   ISn^Sr.'fAlf    0F    H*«    C0Ny"T    12    BIT   BINARY   TO   BCD   TO 
ASCII,    ANO    SENO    IT    SERIALLY   TO   NORTH    STAR 


6216   B6   4000    tOAT 
E219   43 
E21A    97    00 
E21C    B6    4002 
E21F    43 
E220   84    OF 


LOA  A 
COM  A 
STA   A 

LOA  A 
COM  A 
ANO   A 


00RA2 

SUBLO 
00RB2 

ftOF 


READ   LOWER   BYTE   OF  MULTIPROGRAMMER 

l'S    COMPLEMENT    OF    LOW    BYTE 

LOW    BYTE    OF    SUBTRAHEND 

REAO   UPP6R    BYTE   OF   MULTIPROGRAMMER 

l'S   COMPLEMENT    OF    HIGH    BYT6 

MASK   OFF    UPP6R    4    BITS 
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MCE      7 


00319    E222  97   01 

STA    A 

SUSHI 

00320    E224  CE    0006 

WSEND 

LDX 

1LS03 

00321    E227   86    03 

IDA   A 

f*03 

00322    E229   C6    00 

LOA   B 

MOO 

00323    E22B    E7    00 

INC 

STA    B 

o.x 

00324    E220   OS 

I  NX 

00325    E22E    4A 

OEC   A 

00326    E22F   2A    FA 

a  pi. 

INC 

00327    E231    36    E3 

LOA    A 

«*EB 

00323    E233    97    02 

STA    A 

MINLO 

00329    E235    86    03 

LOA    A 

•  $03 

00330    E237   97    03 

STA    A 

MINHI 

00331    E239   BO    E300 

P0S1 

JSR 

SUB16 

00332    E23C   96    04 

LOA    A 

NEO 

00333    E23E    81    OC 

CMP   A 

«»00 

00334    =240   26    03 

BNE 

NEGRI 

00335    E242    7C    0006 

INC 

LS03 

00336    E245   20    F2 

BRA 

P0S1 

00337    E247   96    00 

NECR1 

LOA    A 

SUBLO 

00338    E249   8B    ES 

AOO   A 

»*E8 

00339    E248    97    00 

STA    A 

SUBLO 

00340    E240   96    01 

LOA    A 

SUSHI 

00341    E24F    89    03 

AOC    A 

IS03 

00342   E251   97    01 

STA    A 

SUSHI 

00343    E2S3    86    64 

LOA    A 

It64 

O0344    E255   97    02 

STA    A 

MINLO 

00345    E257   86    OC 

LOA    A 

HCO 

00346    E259  97    03 

STA   A 

MINHI 

00347    E25B    80    E300 

P0S2 

JSR 

SUB16 

00348    E25E   96    04 

LOA    A 

NEC 

0C349    E260   81    00 

CMP    A 

1*00 

00350    E262    26   05 

BNE 

NESR2 

003  51    E264   7C    0007 

INC 

LS02 

00352   E267   20  F2 

BRA 

P0S2 

00353    E269    86    64 

NECR2 

LOA    A 

**64 

003  54    E26S    9S    00 

AOO   A 

SUBLO 

00355    E260    97    00 

STA   A 

SUBLO 

0C3S6    E26F   96    01 

LOA   A 

SUSHI 

00357    E271    89    00 

AOC   A 

JSOO 

003  58    E273    97    01 

STA   A 

SUBHI 

00359    E275    86    OA 

LOA   A 

•  SOA 

00360    E277   97    02 

STA   A 

MINLO 

00361    E279    7F    0003 

CLR 

MINHI 

00362    E27C   BO    E300 

P0S3 

JSR 

SUB  16 

00363    E27F    96    04 

LOA    A 

NEG 

00364    E281    81    OS 

CMP   A 

MOO 

00365    E283    26    OS 

BNE 

NEGR3 

00366    E2S5    7C    0008 

INC 

LS01 

00367   E28S  20   F2 

BRA 

P0S3 

00368    E2SA    96    00 

NESR3 

LOA    A 

SUBLO 

00369    E28C    88    OA 

AOO    A 

fSOA 

00370    E2  8E   97    09 

STA   A 

LSD 

00371    E290  C6    04 

LOA    8 

#*04 

00372    E292    09 

ASCII 

3E.X 

00373    E293    86    30 

LOA    A 

1*30 

HIGH   BYTE   OF    SUBTRAHEND 


CLEAR  CURRENT   DIGIT 


LOAO    1000    INTO    MINUENO 

BRANCH    TO   16    BIT    SUBTRACTION    ROUTINE 

CHECK   TO    SEE    IF    SUBTRACTION    RESULT    IS   NEGATIVE 


RESULT    MAS   NEGATIVE i    ADO    1000    BACK 


LOAO    100    INTO   MINUENO 

BRANCH    TO   16    BIT    SUBTRACT    ROUTINE 


RESULT    HAS   POSITIVE,    GO    SUBTRACT    ANOTHER    100 
RESULT    HAS   NEGATIVE.    AOO    100   BACK 


LOAO    10    INTO   MINUENO 

BRANCH    TO  16    BIT    SUBTRACT    ROUTINE 
CHECK   TO    SEE    IF    THE    RESULT   WAS    NEGATIVE 

RESULT    WAS   POSITIVE,    SUBTRACT   ANOTHER    10 

RESULT    WAS    NEGATIVE,    AOO    10    SACK 
RESULT    IS  LSO 
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RAGE      « 


00374 

00375 

00376 

00377 

00378 

00379 

00380 

00381 

00382 

00383 

00384 

00385 

0C386 

00387 

00388 

00389 

00390 

00391 

00392 

00393 

00394 

00395 

0C396 

00397 

0  0398 

00399 

00400 

00401 

00402 

00403 

00404 

00403 

00406 

00407 

00408 

0C409 

00410 

00411 

00412 

00413 

00414 

0C41S 

0C416 

00417 

0C418 

00419 

00420 

00421 

00422 

00423 

00424 

00425 

00426 

00427 

00428 


E295  AB 
E297  A7 
E299  5A 
E29A  26 
E29C  CE 
E29F  C6 
E2A1  07 
E2A3  A6 
E2A5  81 
E2A7  27 
E2A9  F6 
E2AC  56 
E2A0  56 
E2AE  24 
E280  B7 
E2B3  08 
E2B4  7A 
E2B7  26 
E2B9  39 
E2BA  06 
E2BC  CI 
E2BE  26 
E2C0  86 
E2C2  20 
E2C4  7A 
E2C7  08 
E2CS  20 
E2CA  A6 
E2CC  20 
E300 


00 
00 

F6 

0006 
04 

OS 
00 
30 

11 

8000 


F9 

8001 

OOOB 
II 

08 
01 
04 
30 
E5 
OOOB 

09 

00 
OS 


NEXT! 


E300 
E302 
E303 
E305 
E306 
E309 
E30C 
E300 
E30E 
E310 
E312 
E315 
E317 
E318 
E31A 
E31C 
E31E 
E321 
E322 
E324 


AOO  A 
STA  A 
OEC   B 

BNE 
LOX 
LOA    8 

STA    B 
LOA    A 
CMP   A 
BEO 
L00KC4   LOA    B 
ROR    B 
ROR    B 
SCC 
STA   A 
I  NX 
DEC 
BNE 
RTS 
LOA    B 
CMP    B 
BNE 
LOA    A 
BRA 
OEC 
INX 
BRA 
LOA    A 
8RA 
ORG 


XI  BAT 

ZERO 


0,X 

O.X 

ASCII 

8LS03 

•  104 

KNT 

O.X 

«130 

ZERO 

ACIAS 


L00KC4 
AC  I  AT 

KNT 
PATCH 

KNT 
*S01 
NOT 
«S30 

LOOK C 4 
KNT 

NEXT  I 
O.X 

L00KC4 
SE300 


BCD    TO    ASCII 


GET   READY  TO    SEND    IT   BACK 


INITALIZE   COUNTER   KNT 
COMPARE   THE    CHARACTER   WITH 


REAO   AC  I A  STATUS   REG 


WRITE    CHARACTER    TO    ACIA 


NO   MORE    ZEROS   CAN    8E    SUPPRESSED 


ZEROS    STILL    BEING    SUPPRESSED 
NC   FURTHER    ZERO    SUPPRESSION 


*»*»»»*■******•■**••»*•**■»»», 


•  •***•*••*••***•  »*«*•»*•*  ■«****»  n.« 


16  8IT  SUBTRACTION  SUBROUTINE 


96  03 
36 

96  02 
36 

7F  0003 

7F  000  4 

43 

OC 

8B  01 

24  03 

7C  0005 

06  03 

53 

OB  05 

97  02 
D7  03 
7F  0005 
OC 

96    00 
9B    02 


SU816 


LOA 

A 

MINHI 

PSH 

A 

LOA 

A 

MINLO 

PSH 

A 

CLR 

CARRY 

CLR 

NEG 

COM 

A 

CLC 

AOO 

A 

•  SOI 

BCC 

LI 

INC 

CARRY 

LOA 

B 

MINHI 

COM 

a 

AOO 

B 

CARRY 

STA 

A 

MINLO 

STA 

B 

MINHI 

CLR 

CARRY 

CLC 

LOA 

A 

SUBLO 

ADD 

A 

MINLO 

PUT    HIGH    ANO   LOW   BYTE    OF    MINUENO    ON    THE    STACK 


TWOS'S    COMPLIMENT    OF    MINLO 


TWO'S    COMPLIMENT    CF    MINHI 


SUBTRACT    THE    LOW    BYTES 


100 


NICR01 


MOTOROLA   M68SAM   CROSS-ASSMBLER 


00429 

00430 

00431 

00432 

00433 

00434 

00435 

00436 

00437 

00438 

00439 

00440 

00441 

00442 

00443 

00444 

00445 

00446 

00447 

00443 

00449 

00450 

O0451 

00452 

0C453 

0C454 

00455 

00456 

00457 

0C45S 

00459 

00440 

00461 

00462 

004  63 

00464 

00465 

00466 

00467 

00468 

00469 

0C470 

00471 

00472 

0C473 

00474 

0C475 

00476 

0047T 

00478 

00479 

00480 

00481 

0C482 

00483 


E326  24  03 
E328  7C  000  5 
E32B  06  01 
E320  08  03 
E32F  08  05 
E331  97  OC 
E333  07  01 
E335  2A  03 
E337  73  0C04 
E33A  32 
E338  97 
E330  32 
E33E  97  03 
E340  39 
E400 


02 


BCC 
INC 
LOA  8 
ADO  a 
aoo  a 

STA  A 
STA  a 

a  pl 

COH 
PUL  A 
STA  A 
PUL  A 
STA  A 
XSU816  RTS 
ORG 


L2 


L3 


L2 

CARRY 
SUBH1 
MINHI 

CARRY 

SUBLO 

SUSHI 

L3 

NEC 

HINLO 

MINHI 

SE400 


SUBTRACT    THE    HIGH   BYTES 

STORE   THE   RESULT 

SET    NEG    TO   FF    IF   THE   RESULT    IS    NEGATIVE 

RETRIEVE   THE   ORIGINAL   MINUENO 


E400   7F    0016    TIMER 
E403   80    E18F 
E406    80    30 
E408  48 
E409   CE    E6EE 
E40C   DF    14 
E40E    9B    15 
E410   97    15 

* 
E412    80    E18F 
E415    81    30 
E417   27    04 
E419    86    20 
E418    97    16 
E410   CE    60F0   ARND 
E420   OF    OE 
E422    BO    E16B 

* 
E425    86    00 
E427   97    OF 
E429    86    05 
6423    BO    E476 

* 
E42E    86    01 
E430    BO    E476 

a 
EA33    86    09 
E435    BO    E476 

* 
E438   OE    14 
E43A    EE   00 
E43C    BO    E47E 

E43F  C6   05 


1  HS1"! J°   M0NIT0*    VYF   CARD   FOR    SPECIFIED   TIME      AND    REAO      IT         » 

•  BACK.    ANO    SENO    THE    RESULT    BACK   TC    THE    NORTH   STAR  • 

*******************  »•»*«*»**•»•«»  •*••»•»••»...«*»..,..,.. .»..„.! 

GET    THE    NEXT    CHAR    FRCH   NS    WHEN    SENT 
GET   THE    HEX    EQUIVALENT 
ANO    MULTIPLY    IT    8Y    2 
GET    INTO    X    AOORESS    TO    1ST   COUNT    MINUS    2 


JSR 
SUB  A 
ASL   A 
LOX 
STX 
AOO    A 


GETONE 
IS30 


dCAOO-2 

TEMP 

TEMP*1 


STA    A      TEMP+1 


JSR 
CMP   A 

BEQ 
LOA    A 


GETONE 

**30 

ARNO 

«S20 


STA    A      TYPE 
LOX  JS60FO 


STX 

JSR 

LOA  A 

STA  A 

LOA  A 
JSR 

LOA    A 
JSR 


VALL 
STOHP 

»soo 

VALH 

ttos 

SENO 

**01 
SENO 


LOA    A      »J09 
JSR  SENO 


LDX 

LOX 
JSR 


TEMP 

o.x 

TIMOUT 


6ET    THE    NEXT    CHAR    FRCH   NS    WHEN   SENT 
CHECK    IF    IT    EQUALS    ZERO 


GET    THE    CONTROL    WORO 

SENO    VALUE    IN   VALL    t.    VALH  TO   HP 

OD    -    13    -    SLOT    AOORESS    FOR    V/F   CARD 

SENO    A      05   TO      V/F   CARD 

SENO    A    01    TO    V/F   CARO 

SENO    A    09    TO    V/F    CARC 

GET    THE    POINTER    TO   THE   COUNT 
GET    THE    CORRESPONDING   COUNT    IN    X 


LOA   8     l$05  7060  COUNTER  TO   BE  REAO  5  TIMES 


101 


XICRC1 


MOTOROLA   M68SAN   CROSS-ASSMBLER 


PAGE  10 


00464 
00435 
00486 
00487 
00488 
O0489 
00  490 
0C49X 
00492 
00493 
00494 
00495 
00496 
0C497 
00498 
0C499 
00500 
0C501 
00502 
00503 
00504 
00505 
0CS06 
00507 
0C508 
00509 
00510 
00511 
00512 
00513 
00514 
00515 
00516 
00517 
00518 
00519 
00520 
00521 
00522 
00523 
00524 
00525 
00526 
00527 
00528 
00529 
0G530 
00531 
00532 
00533 
00534 
0C535 
0C536 
00537 
0C53S 


6441    07    14 

£443    86    00 
E445    BO    E476 

E448    86    OF 
E44A   80    E476 

E440   CE    ACFO 
E450    OF    OE 
E452    BO    E16B 

E455    CE    OODO 
E458   OF    OE 
E45A    90    E16  8 

E450    BO  E466 

E460    7A  0014 

E463    26  OE 
E465    39 


E466   SO    E200 
E469   CE    6CF0 
E46C    OF    OE 
E46E    BO    E16B 
E471    86   00 
E473    97    OF 
E475    39 


STA    S      TEMP 


LOA   A 
JSR 

LCA    A 
JSR 

LOX 
STX 
JSR 

LOX 

STX 

JSR 


I  SOD 
SENO 

i*OF 
SENO 

*»AOF0 
VALL 

STOHP 

I 400 00 

VALL 

STOHP 


SEND    A    13   TO    V/F   CARO 


SENO    (SL    ANO    SYE    FOR   REAOING 
SENO    VALUE    IN   VALL    C   VALH   TO    HP 
TO   TRIGGER       INPUT    FROM    V/F   CARO 
SEND    A    9   TO    THE    V/F    SLOT 


JSR  READ 

OEC  TEMP 

BNE  AGAIN 
XTIMER    RTS 

*  * 

*  ROUTINE  USED  BY  TIMER  TO  READ  A  BYTE  FROM  V/F  CARO  S  SENO  TO  NS  • 

*  « 

TIMEOUT,  REAO  HF  ANO  SENO  TO  NSTAR 


JSR 

RONSNO 

LOX 

f J60F0 

STX 

VALL 

JSR 

STOHP 

LOA 

A 

*>00 

STA 

A 

VALH 

SENO  VALUE  IN  VALL  t  VALH  TO  HP 


XREAO   RTS 


E476   9A    16 
E478    97    OE 
E47A   BO    E16B 
£470    39 


E47E   86    SO 
E480   4A 
E481    26    FO 
£483    09 
E484   26    F 8 
E486    39 


ROUTINE   USED    8T    TIMER    TO    ACORESS    THE    V/F    SLOT 


SENO        ORA    A  TYPE 

STA   A  VALL 

JSR  STOHP  SENO   A    5    TO    THE    V/F    SLOT 

XSENO      RTS 

•  *•**  "a  «»a»a»«*»*»**. *«»■*«**  »n<*l«)*M>tsit«M»tMttllti»t«*»  a: 

»  TIMER    SUBROUTINE    TO   TIHECUT   FOR    1/10.    U,    OR       10    SEC    DEPENDING      o 

•  ON    THE   CONTENTS    OF    REG   X  • 

•  o 
********************** **«*a*«w***«**a*«*s*****»a**««««a*as**a* ***«»»«« 
TIMOUT    LOA  A      #150 

BACK        OEC   A 

BNE  BACK 

DEX 

BNE  TIHOUT 

XTIMOU    RTS 
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00539 

00540 

00541 

0CS42 

00543 

0054* 

00545 

00546 

005*7 

00548 

0C549 

00550 

00551 

00552 

00553 

00554 

00555 

00556 

00557 

00558 

C0559 

0C560 

00561 

00562 

00563 

00564 

00565 

00566 

00567 

0  0568 

0C569 

00570 

00571 

00572 

00573 

00574 

00575 

0CS76 

00577 

00578 

0C579 

00580 

0CS81 

0CSS2 

00583 

00584 

00585 

00586 

00587 

00588 

0C589 

00590 

0C591 

00592 

00593 


•  RELAY      SUBROUTINE 

!  rIHi!.J2UTINE   l"LL   SSTABLISHIHRITE)    -OH-   *EAO    BACK   RELAY    POSITIONS 

•  TO    WRITE    :-    IT    ACCEPTS    NEW    POSITION    PROM    N.STAR    .CONVERTS    Tn 

I      %mF.l£?.SZal    IT   0UT   0N    SL°T    A      (Soon    OF    HP    ?NPUT    PORT 

•  TO      REA01-   FETCHES    CONTENTS    OF   SLOT    A   OF    HP    ANO    SENOS    IT    TO    NS 


ORG 


E500 

E500   BO    E18F    RELAY      JSR 

E503    81    57  CMP 

E505    26    20  BNE 

E507   80    E10B  JSR 

E50A   OE    OE  LOX 

E50C   OF    17  STX 

» 

E50S   CE    70F0  LOX 

E5U    OF    OE  STX 

E513    BO    E16B  JSR 

• 

E516   OE    17  lox 

E518   OF    OE  STX 

E51A    96    OF  L0A    A 

E51C    84    OF  ANO    A 

E51E    8A    10  ORA    A 

E520   97    OF  STA    A 

E522   80    E16B  JSR 

E525    20    13  BRA 

• 
S527   CE    AOFO    ROREL      LOX 

E52A    OF    06  STX 

E52C   BO    E16B  JSR 

• 

E52F  CE    0C1  0  LOX 

E532   OF    OE  STX 

E534    BO    E16B  JSR 

E537  BO    E200  JSR 

E53A    39  XREL  RTS 

E600  0RS 


"        "iimimi 

SE500 


'••••••»»•••«•»»•««»..«. 


GETONE 
<*57 

ROREL 
NSREAO 
VALL 
TVALL 

J870F0 

VALL 

STOHP 

TVALL 

VALL 

VALH 

(SOF 

M10 

VALH 

STOHP 

XREL 

ISAOFO 

VALL 

STOHP 

•SOOIO 
VALL 
STOHP 
RONSNO 

IE600 


GET    THE   NEXT    CHAR    FRCM   NS    WHEN    SENT 

IF   CHARACTER    IS    NOT    «W«    THEN    REAO   RELAY 

ELSE    GET   NEW   RELAY    POSITIONS      FROM        NS 
SAVE   THE    RELAY    POSITIONS    FOR    LATER    USE 
IF    NOTHING    SENT    THEN   OPEN    ALL    RELAYS 

SEND   OTE    SYE    ANC    THE    TO   HP 

SENO   VALUE    IN    VALL    t   VALH   TO    HP 

GET   THE    SAVEO   RELAY   POSITIONS 

ATTACH    SLOT    A   AOORESS   TO    OATA 

SENO   VALUE    IN   VALL    C   VALH   TO    HP 
OONE    -    SO   EXIT    FRCM    HERE 

SENO    ISL    AND    SYE    TO    TRIGGER    HP    FOR    INPUT 

SENO   VALUE    IN   VALL    C   VALH   TO    HP 

REAO    SLCT    A    I0001-CONTROL) 

SENO   VALUE    IN    VALL    C    VALH    TO    HP 

REAO    RELAY    ANO    SENO      RESULT    TO    N.STAR 


E600   BO    E10B 
5603    96    OE 
E605    84 
E607   48 
E608   48 
E609   SA    _ 
E608    97    16 
6600    7F    0019 
E610    BO    E10B 


vco 


03 


BO 


_,..    .  V    C   0         SUBROUTINE  ««».»»•»... 

"*C«W  «•»   CHARACTERS    SUBROUTINE 


JSR 
LOA 
ANO 
■ISL 
ASL 
ORA 
STA 
CLR 
GcTFRQ    JSR 


NSREAO 

VALL 

»J03 


IF    RANGE   NOT    BETWEEN    0    AND 
THEN    RETAIN    ONLY    THE    LEAST 


BITS 


I1B0 
TYPE 
FFLAG 
NSREAO 


IT    WITH    SLOT    AOOR 


POSITION    RANGE    ANO    COMBINE 
SAVE    IT    FOR    LATER    USE 

r^°ri2°iCAT,NC    N0    CH*NGE    W    "EOUENCY 
GET    THE    MAGNITUDE   OF    FREO    IN    SANARY 
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MOTOROLA    M68SAM   CROSS-ASSHBLER 


PAGE  12 


00594    6613    96 

10 

LOA 

A 

CP 

OC595    E615    26 

03 

8NE 

AROUN 

00596    E6I7   7C 

0019 

INC 

FFLAG 

0C597    E61A   96 

OE         AROUN 

LOA 

A 

VALL 

00598    E61C    97 

17 

STA 

A 

TVALL 

00599. E61E   SO 

E10S    GETAMP 

JSR 

NSREAO 

00600    E621    96 

10 

LOA 

A 

CP 

00601    E623   27 

IF 

BEQ 

NOCHNG 

00602    E625    06 

OE 

LOA 

B 

VALL 

00603    6627    37 

PSH 

3 

00604    E623   CE 

60F0 

LOX 

IS60F0 

00605    E62B    OF 

OE 

STX 

VALL 

00606    E620   BO 

E16B 

JSR 

STOHP 

00607 

* 

00608    E630   33 

PUL 

B 

00609    E631    02 

NOP 

C061O    E632   07 

OE 

STA 

3 

VALL 

00611    £634   96 

16 

LOA 

A 

TYPE 

0C612    E636    8A 

01 

ORA 

A 

•  to  I 

00613    EMS    97 

OF 

STA 

A 

VALH 

00614    E63A    30 

E16B 

JSR 

STOHP 

0C61S 

* 

00616    E630    96 

16 

LOA 

A 

TYPE 

0C617    E63F    97 

OF 

STA 

A 

VALH 

00618    E641    SO 

616  3 

JSR 

STOHP 

0061 9 

* 

0C620    E644    96 

19         NOCHNG 

LOA 

A 

FFLAG 

00621    E646    26 

10 

BNE 

XVCO 

0  0622    E64S   CE 

6  of  a 

LOX 

•  t60F0 

00623    E64B    OF 

OE 

STx 

VALL 

00624    E640    SO 

E16  3 

JSR 

STOHP 

0C62S 

* 

00626    E650    96 

17 

LOA 

A 

TVALL 

0C627    E652    43 

COM 

A 

0062B    E653    97 

CE 

STA 

A 

VALL 

00629    E655    96 

16 

LOA 

A 

TYPE 

00630    E657    9A 

02 

ORA 

A 

<»02 

00631    E659   97 

OF 

STA 

A 

VALH 

00632    E65B    BO 

616  3 

JSR 

STOHP 

0C633 

* 

00634    E65E    96 

16 

LOA 

A 

TYPE 

00635    E660   97 

OF 

STA 

A 

VALH 

00636    E662    3D 

E16B 

JSR 

STOHP 

00637    E665    39 

XVCQ 

RTS 

00638    E6F0 

ORG 

IE6F0 

00639    E6F0    0238           CAOO 

F08 

10288 

0C640    E6F2    1958 

FDB 

11958 

00641    E6F4    F070 

FC8 

tF070 

00642    E7FE 

ORG 

tE7FE 

00643    E7FE    EO 

FC8 

tEOi tOO 

E7FF    00 

0C644 

END 

SYMBOL   TABLE 

IF  NOTHING  SENT  THEN  ASSUME 
NO  CHANGE  IN  FREQUENCY  OF  VCO 

SAVE  IT  FCR  NOW  AND 

GET  THE  AMPLITUOE  FROM  NORTH  STAR 

IF  NOTHING  SENT  THEN 

DO  NOT  CHANGE  THE  AMPLITUOE  OF  VCO  OUTPUT 

SAVE  ONLY  LAST  8  BITS  OF  AMPLITUOE 

ON  STACK 

SEND      SYE   ANO   DTE    TO    HP 

SEND   VALUE    IN   VALL    C    VALH   TO    HP 

GET    AMPLITUOE    IN   VALL 
MAY    REQUIRE    A   COM    B    HERE 

GET    SLOT    ADORESS    WITH   RANGE 
MAKE    WR(AMP)     LOW 

SEND   VALUE    IN    VALL    S    VALH    TO    HP 


HOLD    AMPLITUDE    DATA 

SEND   VALUE    IN    VALL    C    VALH   TO    HP 


IF    NO   CHANGE    IN    FRO    THEN    QUIT 

SENO   VALUE    IN    VALL    C    VALH    TO    HP 
GET    THE    GREQUENCY    MAGNITUOE 


GET    SLOT    ADDRESS    A80   RANGES 
MAKE    HRIFRQI    LOU 


SENO   VALUE    IN   VALL    C    VALH   TO    HP 


HOLD    FREQUENCY    DATA 

SENO    VALUE    IN   VALL    C    VALH    TO    HP 
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BQB1EICBIIQM  HSCE  lb  IHE  E2ISIING  EBSQB  IB  "5  §823  INISSEBBI 
The  following  modification*  war*  mad*  in  tha  existing 
assembly  program  in  the  MC  6S0S  interface  EPROM.  The 
modifications  Mara  mad*  in  the  ADC  assembly  subroutine  so  that 
the  ADC  could  be  triggered  by  an  external  signal  and  samplas  be 
taken.  The  memory  address  at  which  the  modifications  ware  mad* 
ar*  clearly  indicated.  Tha  original  listing  of  the  assembly  is 
also  shown.  Tha  naw  EPROM  contents  are  also  shown. 
MEMORY  LOCATION  El  AS 

JMP      MODPRGMST      !  JUMP  TO  THE  MODIFIED  PROGRAM  ADDRESS 
MEMORY  LOCATION  E7  la 

ORG      MODPRSMST       !  START  OF  THE  MODIFIED  PROSRAM  AT  E7  lu 
LDX  asfcmaa 

!SET  INTERRUPT  MASK 

!  SET  BIT8  TO  ENABLE  INTERRUPT 

!  CLEAR  BIT7  OF  CONTROL  REG. 
I  CLEAR  INTERRUPT  MASK 
!  UAIT  FOR  INTERRUPT 
'  SET  INTERRUPT  MASK 
t  GO  BACK  TO  THE  ORIGINAL  PRGM. 
1  CONTINUE  ORIGINAL  PRGM.  AT 
!  RETURN  ADDRESS  El  AB 
As   soon   as  an  interrupt  occurs  tha  microprocessor   fetches 
tha  interrupt  service  vector  and  this  is  initialized  at   location 
E7  F8.   Tha  contents  of  this  location  are  E5  3a  and  tha  interrupt 


SEI 

LDA  A 

CRB2 

ORA  A 

•*ei 

STA  A 

CRB2 

LDA  A 

DDRB2 

CLI 

uo: 

SEI 

JMP 

RTNADDR 

RTNADDR  IS  El 

AB 
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t«rvic*  routine  starts  at  this  location  ia.  ES  Se. 
INTERRUPT  SERVICE  ROUTINE  STPRTS  PT  ES  58 


INTROUTIN 

NOP 

LDP.  A      DDRBS 

LDfl  P      CRBS 

AND  R       »«FE 

STB  A       CRBS 

RTI 


STPRT  OF   INTERRUPT  ROUTINE (ES  5B> 

NO  OPERATION 

CLEAR  BIT7  OF  CONTROL  REG. 

LOAD  ACC.  WITH  [CONTROL  RES.  B. 3 

CLEAR  BITB  TO  DISABLE  INTERRUPTS 


!  RETURN  FROM  THE  INTERRUPT  ROUTINE 
Thus  thasa  modifications  mada  in  tha  existing  ADC  subroutine 
enabled   the   use  of  a  trigger   reference   signal   to   initiate 
sampling.  The  aodified  EPROM  contents  are  shown  on  the  next  page. 


CBfilSCSli  RE  Z--  !20D'''r:E3 


333 

00  5F  00  17  00  03  00  OS 

8E  88  SO  85  88  B7  28  01 

310 

97  20  03  36  FF  87  20  00 

B7  20  02  86  04  B7  28  33 

£23 

35  3C  37  20  01  85  00  37 

40  31  37  40  03  37  48  00 

2123 

B7  40  02  SS  04  B7  40  03 

86  3C  B7  40  01  36  83  37 

944 

SO  00  BS  01  B7  80  08  B6 

40  82  64  30  44  44  44  3B 

333 

00  37  12  96  EO  89  00  97 

11  OE  11  EE  GO  36  93  B7 

360 

A3  01  88  00  B7  flo  00  FF 

A3  84  86  33  37  08  01  3D 

373 

21  3F  81  26  26  OS  BD  si 

83  20  F4  81  3F  25  83  BD 

333 

E2  00  20  EB  81  23  26  33 

BD  E4  80  28  E2  81  23  26 

393 

3!  BD  E3  00  20  D9  81  21 

26  OS  BD  E6  08  23  08  31 

WW 

41  26  CC  3D  £1  93  20  C7 

3E3 

F~  W  FF  FF  FF  FF  FF  FF 

FF  FF  FF  FF  *F   FF   F*   *F 

3C3 

FF  ri-    FF  FF  FF  FF  FF  FF 

FF  FF  FF  FF  FF  IF  FF  FF 

303 

FF  =F  FF  FF  FF  FF  FF  FF 

FF  F=  "  FF  FF   F=  FF  FF 

320 

FF  FF  FF   FF  FF  FF  FF  FF 

CF3 

~  FF  FF  FF  FF   FF  FF  FF 

FF  FF  FF  FF  FF  FF  FF  FF 

133 

BD  El  03  36  10  27  03  BD 

El  S3  39  86  08  37  18  37 

1  13 

OE  97  OF  37  OC  97  OE  37 

OS  4C  37  30  3D  El  8F  3: 

:23 

24  27  08  80  30  36  7C  00 

13  20  Fl  36  10  27  33  22 
OB  36  10  31  OB  27  2B  20 

139 

D6  OB  2S  CD  97  OE  7C  08 

-  14a 

EE  SF  OC  48  19  97  OC  D7 

3D  48  S3  48  S3  OC  9B  3C 

:so 

D9  3D  7S  00  03  26  E3  9B 

OE  D9  OF  97  3£  D7  8F  7C 

i£3 

00  OB  96  00  97  83  91  18 

25  CS  33  36  OE  37  20  00 

173 

96  8F  37  20  02  02  82  36 

48  81  84  F7  37  40  81  B6 

133 

40  02  84  40  27  F9  B6  48 

01  80  88  B7  48  01  33  B6 

133 

80  00  47  24  FO  B6  80  01 

33  FF  FF  CE  BO  FO  DF  OE 

ISO 

BD  El  SB  CE  08  SO  DF  OE 

7E  67  10  BD  El  6B  36  40 

133 

00  A7  00  88  36  48  02  07 

08  08  SC  54  00  26  EC  CE 

:C3 

60  88  06  88  43  37  88  08 

06  38  43  34  OF  97  01  88 

108 

DF  IS  CS  80  3F  BD  E4  7E 

3D  E2  24  3D  E2  89  DE  IS 

1E3 

8C  64  00  26  D3  33  FF  CE 

68  08  B6  40  03  20  F3  B6 

1F3 

40  02  7S  El  OB  FF  FF  FF 

FF  FF  FF  FF   FF  FF  FF  FF 

220 

CE  38  34  3D  E4  7E  3D  E2 

16  B6  38  08  46  46  24  F3 

213 

3c  23  B7  80  81  33  B6  40 

08  43  97  88  B6  48  02  43 

239 

84  OF  37  01  CE  00  06  86 

83  CS  38  E7  33  38  40  20 

220 

FO  86  E3  37  32  36  03  37 

33  3D  23  30  35  84  81  88 

243 

26  OS  7C  03  06  28  F2  =6 

30  3B  £3  97  00  96  31  83 

£38 

03  37  01  36  64  37  82  36 

33  97  03  3D  33  33  95  34 

2S3 

81  80  25  33  7C  08  87  23 

F2  36  64  93  38  97  88  96 

270 

01  83  00  37  81  86  80  37 

32  7F  80  03  BO  S3  03  3S 

230 

84  81  03  26  OS  7C  08  88 

23  F2  36  83  33  00  37  03 

290 

CS  04  39  86  38  SB  08  07 

38  SO  25  FS  OE  03  OS  CS 

11  r6  30  30  36  35  24  pg 
11  33  DS  3B  CI  31  25  ?4 
23  03  06  00  2^   03  fs  <rr 

MS 

133 

34  D7  3B  06  88  81  33  27 
B7  80  01  08  70  08  OB  36 

2C0 

85  33  23  23  70  80  SB  38 

220 

"F  FF  FF  FF  FF  FF  FF  FF 

FF  FF  FF  FF  FF  FF    Fr    « 

223 

FF  FF  FF  FF  FF  FF  FF  FF 

FF  FF  FF  FF  —  FF  F=  FF 

2F0 

FF  FF  FF  FF  FF  FF  FF  FF 

FF  FF  FF   FF  FF  FF  FF  FF 

300 

36  33  36  36  32  36  7F  08 

03  7F  88  04  43  OC  33  01 

310 

24  83  7C  00  OS  06  83  S3 

DB  OS  37  02  D7  03  7F  83 

320 

OS  OC  36  88  3B  02  24  03 

7C  88  33  D5  31  DB  03  DB 

330 
340 
333 

3S  37  08  D7  01  20  83  73 

33  F^  -s  FF  FF  FF  FF  FF 

30  04  32  37  32  32  37  03 

2S0 
373 
383 
39* 

FF  FF  FF  FF  FF  FF  FF  F* 
FF  F-  FF  FF  FJ»  FF  FF  F- 

FF  FF  FS   r?  rr  cp  sr  re 

II  II  II  II  'F  '"    "    -rF 

F=  FF   FF  F-    FF   FF   FF   FF 

303 

350 

--  ~   ::  FF   ~   il   7-   ~ 

ssjSSSssS 

188 


3C3 

=  = 

FF 

" 

-F 

F- 

rr 

F= 

rr 

FF 

-r 

rr 

=  F 

-- 

rr 

rr 

rr 

2u3 

FF 

=  = 

FF 

FF 

FF 

F- 

F  = 

-- 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

-F 

323 

r* 

rr 

FF 

FF 

?F 

FF 

FF 

FF 

FF 

" 

-tr 

-F 

FF 

FF 

Fr 

F- 

3.-0 

F~ 

rr 

FF 

FF 

FF 

FF 

FF 

rr 

FF 

FF 

FF 

rr 

rr 

FF 

rr 

-F 

.03 

7- 

33 

16 

3D 

21 

3F 

80 

30 

4a 

C2 

25 

22 

DF 

14 

93 

13 

-13 

37 

13 

BD 

21 

3F 

31 

30 

27 

34 

35 

23 

37 

16 

C2 

60 

F3 

in 

DF 

SE 

3D 

21 

53 

36 

DO 

37 

OF 

35 

03 

BD 

24 

75 

86 

31 

433 

BD 

E4 

75 

36 

39 

3D 

24 

76 

DE 

14 

EE 

30 

3D 

24 

72 

C5 

444 

33 

D7 

14 

B6 

3D 

3D 

E4 

76 

36 

3F 

3D 

24 

76 

C2 

SO 

FO 

433 

DF 

3E 

3D 

El 

53 

22 

00 

DO 

DF 

32 

3D 

21 

53 

BD 

24 

66 

463 

7Q 

00 

14 

25 

D2 

33 

BD 

ro 

00 

C2 

60 

FO 

DF 

M 

BD 

El 

473 

S3 

36 

DO 

37 

OF 

23 

30 

16 

97 

OE 

BD 

21 

53 

29 

36 

ID 

4aa 

4fl 

26 

FD 

03 

25 

F8 

39 

FF 

FF 

FF 

FF 

FF 

rr 

FF 

FF 

FF 

438 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

FP 

FF 

FF 

FF 

403 

FF 

FF 

FF 

FF 

FF 

-F 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

430 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

4C0 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

400 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

cr 

FF 

FF 

4Ce 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

-F 

F= 

FF 

FF 

FF 

FF 

4FB 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

FF 

F= 

FF 

FF 

FF 

FF 

FF 

FF 

saa 

BO 

2: 

3F 

31 

37 

25 

20 

BD 

El 

03 

DE 

OE 

DF 

17 

C2 

70 

sia 

F3 

DF 

N 

3D 

21 

53 

DE 

17 

DF 

M 

36 

OF 

34 

3F 

30 

10 

ssa 

97 

OF 

3D 

El 

63 

20 

13 

C£ 

SO 

FO 

DF 

3E 

3D 

21 

63 

C2 

53a 

00 

10 

DF 

BE 

3D 

21 

63 

3D 

22 

00 

39 
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APPENDIX  4 


TABLE  SH0WIN6  THE  READIN6S  OBTAINED  FOR  THE  ADOPTIVE  METHOD 


FREQUENCY 

CALCULATED  PHASE 

ACTUAL  PHASE 

ERROR 

F(Hi.  ) 

P(Deg) 

(Deg) 

(Deg) 

200 

60. 842437 

60.78165 

0.739213 

300 

50. 93527S 

50.013108 

0.922164 

400 

41.2S0897 

41.807545 

0. 526648 

588 

34.611122 

35. 58437 

0.  973248 

SM 

30. 582739 

30. 807028 

0. 224289 

700 

26. 862562 

27. 072203 

0.209641 

SM 

23. 599043 

24. 095689 

0.  496646 

900 

21.049831 

21.679956 

0.630125 

1000 

19.271924 

19.686781 

0.414857 

£000 

10.005898 

10. 142929 

0. 137031 

The  actual  phase  was  obtained  from  the  RC  circuit  used.  The 
values  of  R  and  C  were  found  using  a  digital  bridge. (Genrad  RLC 
Digibridge).  R=101. 15K   and   C=4. 394  nF. 
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APPENDIX  5 
Plots  obtained  for  the  FFT  method  with  actual  data  samples. 


1.1 


.E 


.A 


.3 


1. 


AMPLITUDE    X 


- — - 


FREQ.    HINDOH-0   TO   4.56300E+B3 
F  REQ .     I NTER VRL-7 . 1 3  2  8 1 E+0 1 
MRX   VflLllE=9.93624E-0] 


A. 


3.0 


p 


-i- 


-v 


—r 


0  a  IS  24  32  4<2  AS  5S  G4 

COEFFICIENT 

I 1 1 1  w» *iTJ= 1 1  I  I 

2  5.71E+02  1.14E+03  1.71E+Q3  2.28E+03  2.85E+03  3.«2E-H23  3.3SE+B3  4.5SE+B3 


FREQUENCY  CHxl 


111 


IBB 

150 

120 

SB 

SB 

3B 

I 

-3B 

-SB 

-30 

-iaB 

-150 
-IBB 


PHRSE  OF  X 


FREQUENCY  DOMAIN  DRTH  (PHASE) 

fREO.  WINDOH=0  Tip  4  .56300HN-B3 
f  REQ .  I  NTER VflL-7 1  1 3  2  ell  EtQJ 


e 

B   5 


16 


24  32 

COEFFICIENT 


46 


43 


36 


64 


71E+02     ]  .  ME+03     1.71E+03    2.2BE+03    2.B5E+03    3.H2E+B3    3.3SE+03    -1.3SE  +  B3 

FREQUENCY    CHz] 
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1.3 


.3 


0.0 


AMPLITUDE  Y 


FREQ.  WINDOW-0  TO  4.58300E-*B3 
F  REQ .  I NTERVRL-7 . 1 3  2  8 1 E+0 1 
MRX  VALUE-8.85345E-01 


^v 


-A- 


^- 


1 1 

1                     8 

i                       1 

. , 

18 
I 

34 

1 

1 

32 
COEFFICIENT 

40 

48 
1 

se 

i 

E4 

0  5-.71E+02  1.14E+Q3  1.7JE+03  2.28E+03  2.8SE+Q3  3.42E+B3  3.33EtB3  4.SBE+03 

FREQUENCY  CHzJ 
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FREQUENCY  DOMAIN  DHTH  (PHASE) 

FREQ.  WINDOW=e  TO  4  .S6509E-JJ03 
FREQ.  INTERVAL-? . 1328 !E*21 


.71E+02  1.14E+03  1.71E+B3  2.29E+03  2.83E+03  3.42E+03  3.39E+03  1.5EE+03 

FREQUENCY  CHzl 
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Printout    of   the   results   obtained    for   X 


FREQUENCY    DOMAIN    DATA 


FREQUENCY    UINDOU-0    TO    4. 34300E*03[ Hi ] 
FREQUENCY    INTERVAL"   7. 13281E->0UHt  J 


DC    TERM 

MAX    FREQ. 

I  • 


16 
17 
19 
I? 

SO 


23 
it 

27 
28 
2? 
30 


32 
33 

34 
35 
36 
37 


42 
43 


50 
SI 

32 
33 
34 

33 

36 
37 


FREQUENCYIHlI 
O.OOOOOE-f-OQ 
4.36500E-03 
7.13281E4-01 
1 .42636E-02 
2.13984E»02 
2  9S312E»02 
3.36641E-02 
4.27969E-*02 
4 .  99297E*02 

3.70623E-02 
6.41933E»02 
7.  I3281E-02 
7.84609E+02 
9.33937E»02 
».27244E*02 
9.98594E4-02 
1 . 04992E-.03 

I.14I23E*03 
1 .21238E*03 
I .28391E»03 
t .  35323E*03 
1.42434E*03 
1 . 497S9E»03 
1 . 36»22E«-03 
I  .6403IE-03 

I.7U87E-03 
1 . 78320E-03 
1  . 93433E-3  3 
t  .  923S6E4-03 
1.99719E»03 
2.06S32E*03 
2.13984E»03 
2.211t7E-HI3 

2.29230E-0  3 
2.33383E+03 
2.423I6E-*-03 
2.49648E+03 
2.36781E*03 
2.63914E*0  3 
2 .  71 047E-I-03 
2.781B0E+03 

2.83312E-03 
2.92443E+03 
2.99378E*03 

3. 3671 :E-0  3 
3.13S44E>03 
3.20977E-H3 

3.28109E-03 
3.332  42E.0  3 

3.42373E-H3 

3.49308E*03 
3. 36641E+03 
3 . 63773E»03 
3.7J9O6E.03 
3.79039E+03 
3.83172E*03 
3.92393E«03 

3.99437E*03 
4.06370E»03 
4.t3703E*03 
4.20834E*03 
4.27949E+03 
4.33102E-03 
4.42234E-03 
4.49347E»03 


REAL 
J.40134E-02 

-I.76363E-03 
1    22233E-03 

-1  99268E-04 
8.24246E-04 
8.33S04E-04 
3.80314E-04 
1.37764E-03 
1.S6304E-03 

-4.86491E-01 
-3.81833E-03 
-4.73131E-04 
-9.93944E-04 
-1 . 00074E-03 
-6.69683E-05 
-1 .28368E-03 
-1.11979E-03 

2 . 72667E-03 

2.03029E-03 
-1.17371E-04 
I . 39866E- 0  4 
6.49400E-04 
-2.38147E-04 
7.86S87E-04 
I . 43348E-03 

9.20344E-04 
-2.34123E-03 
-6.2922BE-04 
-3 . 97233E-04 
-I.I0891E-03 
-2.33238E-04 
-9.11379E-04 
-2.33126E-03 

4.21973E-04 
2.09620E-03 
6.83043E-04 
-1  .73746E-03 
B.91834E-04 
1.23044E-04 
3.98217E-04 
2.24207E-O3 

-3.34938E-02 
-1.74942E-03 
-1.04247E-03 
-3.46944E-03 
-9.46894E-04 
-4.398I7E-04 
-2.06128E-U4 
-2.41999E-03 

-1.64833E-03 
7.30390E-04 
8.92867E-04 

-3.08087E-04 
4.51612E-04 
2.90333E-04 

-3.7O1O3E-04 
1 .643S3E-03 

1 . 08919E-02 

2.017S3E-03 
-1  .02150E-03 

2.72912E-04 
-2.71833E-04 
-4.95025E-04 

3.43002E-04 
-1.13086E-03 


MAGNITUDE    PHASEIDECJ 


1 .87971E-03 
6.49393E-04 
2.80697E-05 
7.07762E-04 
-1  .  06828E-04 
-1 .  70660E-04 
7.391B1E-04 

3.39146E-01 
7.48294E-0S 

-3.71229E-04 
7.23897E-04 

-2.31644E-04 
4.64290E-03 
2.34433E-04 

-1.99492E-03 

-1.3S389E-S3 

1  .  06643E-O3 
I .20100E-03 

-I.42109E-03 
7.67064E-04 
3.99868E-04 

-1 . 0343SE-04 
2.02627E-03 

-2.20326E-02 
1.66U3E-04 

-7.83747E-04 
4.04247E-04 

-1.374I4E-04 

-3.10299E-04 
4.77342E-04 

-9.72983E-04 

3.00000E-04 
-3.84B44E-04 

7.26142E-04 
-3.32I70E-04 
-I.33993E-04 

3.41S47E-04 
-3. 3366 0E- 04 

1 -20244E-04 

-4 . 32263E-03 
1.90004E-03 

-2 . 0S737E-04 
5.16044E-04 
6.26218E-04 

-1.70736E-04 
1 . 07738E-03 
8.77726E-04 

-3.46327E-03 
-2.24680E-03 
-8.33841E-03 
-3.16747E-04 
-9. I3781E-04 
8.68241E-03 
-9. 14404E-04 
-1  .59044E-03 

-1.37I74E-02 
2.39334E-03 
4.S61B0E-04 
2.S9780E-04 
9.99892E-04 
6.24499E-03 
7.33391E-04 
2.14277E-03 


2.24219E-03 
4.79280E-04 
8.24724E-04 
1 . 09499E-03 
3.93034E-04 
1 . 38484E-03 
2.00618E-03 

8.83343E-01 
3.81904E-03 
6.0296  0E-04 
1.22962E-03 
I.02720E-03 
8.14889E-03 
1.30863E-03 
2.29946E-03 

3.04430E-03 
2.31106E-03 
1 .20674E-03 
1.39391 E-04 
1  . 00304E-03 
6.S3063E-04 
7.934S4E-04 
2.49347E-03 

2.20718E-02 
2.34667E-03 
1  00664E-03 
S.43182E-04 
1 . 1 1739E-03 
4.0033aE-a4 
1 .133S3E-03 
2.S2M6E-03 

4.54201E-04 
2.27609E-03 
9.94923E-04 
3.32424E-04 
9.  03031  E-04 
3.82221E-04 
9.  23491  E-04 
2.24329E-03 

3 . 37736E-02 
2.58299E-03 
1 . 06337E-03 
3.17209E-04 
1  . 13323E-03 
4.71801E-04 
1 . 09692E-03 
2.57413E-03 


3 . 83362E-03 
2.38778E-03 
8.96771E-04 
4.  418475-04 
9.30493E-04 
3.03037E-04 
9.84464E-04 
2.28873E-03 

1.91223E-02 
2.39362E-03 
1  . 13129E-03 
3.98042E-04 
1  . 03413E-03 
4.98974E-04 
9.29481E-04 
2.42287E-03 


36.94 

107.04 

1.93 

40.27 
-13.49 

-4.17 

21.42 

140.83 
178.33 

-142.00 
143.93 

-164.97 
143.27 
148.79 

-119.28 

-24.41 

27.48 
93.39 
-5.84 
49.73 
113.29 
-7.49 
34.33 

-37.41 

174.26 
-128.49 

134.36 
-172.94 
-129.22 

1 43 . 38 
-137.33 

49.84 
-22.93 

44.73 
-92.99 

-9.80 

71.22 

-44.32 

3.07 

-172.43 
132.44 

-169.68 
93.83 
146.32 

-138.78 
100.83 
140.04 


-113.43 

-71.48 

-3.33 

-134.21 

-60.97 

16.63 

-112. 04 

-44, 02 

-33.28 
89.32 

134.33 
46.72 

103.21 

172.79 
54.22 

117.82 
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Printout  of  the  results  ootainea  for  Y 

FREQUENCY    OOfldlN    DATA 


FREQUENCY    UINDOU-0    TO    4. 54500E»03[Hl! 
FREQUENCY    INTERVAL-    7.  13281EtOHHll 


COEFF. 

FREQUENCYIMll 

REAL 

MM 

MAGNITUDE 

PHASE! DEC  1 

DC  TERN 

0.0O0O0E*OO 

4.13313E-02 

NAX  FREQ. 

4.36330E-03 

-3.23438E-03 

1 

7.13281E-01 

1.44399E-03 

7.4q334E-04 

1.32126E-03 

23.99 

2 

1 . 42434E*02 

4.38783E-04 

1.03810E-03 

1.24443E-03 

38.09 

3 

2.13984E*02 

7.1S224E-05 

4.28124E-04 

6.32220E-04 

8  3.48 

4 

2.33312S-02 

I.49137E-04 

-3 .  44704E-04 

3.94021E-04 

-47.74 

S 

3. 36641 E« 0  2 

S.77433E-04 

-1.27403E-03 

1.39886E-03 

-43.41 

6 

4.27969E*02 

1 .23S34E-03 

-9.23392E-04 

1.33S31E-03 

-34.34 

7 

4.9»297E»02 

2.34984E-03 

3.314SSE-03 

2.37008E-03 

.71 

9 

3.70623E-02 

-3.02S29E-01 

8.42874E-01 

B.93624E-01 

109.74 

9 

6.419S3E»02 

6.80731E-04 

2 . 73284E-03 

2.81433E-03 

76.01 

10 

7.1 3231 E»02 

-7.42373E-04 

2.03421E-03 

2.16S30E-03 

110. 03 

11 

7.34609E-02 

-9.44134E-04 

1.27049E-03 

1.392BBE-03 

1 26 . 62 

12 

8.33937E»02 

-1 .8B722E-04 

4.74429E-04 

3.103B7E-04 

1 1 1  .  49 

13 

9.27244E»02 

3.90431E-04 

3.41098E-04 

6.81994E-04 

30.01 

14 

9 . 98594E*02 

9 .  94977E-0-4 

9.87334E-04 

1 . 40329E-03 

44.73 

13 

1 . 06992E»03 

O.92318E-04 

1.77414E-03 

1.90632E-03 

48.70 

It 

1.14123E*03 

-3.0002BE-03 

-3.12894E-03 

3.89836E-03 

-147.96 

17 

1 .21238E-03 

-9.44978E-04 

1.84333E-03 

2.09033E-03 

116.94 

18 

1.29391E-03 

-I.27492E-03 

9.6032SE-04 

1 . 3962SE-03 

143.  01 

1? 

1 . 33322E-03 

-8.74743E-04 

2.72213E-04 

9.18031E-04 

142.73 

20 

I -42636E+03 

-1.1I337E-04 

1.74192E-04 

2.08422E-04 

122.29 

21 

1 . 49789E*03 

3.89231E-04 

7.03330E-04 

8.03417E-04 

41.11 

22 

1 . 36922E-03 

1 .44240E-04 

1.S4439E-03 

1.3S330E-03 

33.93 

23 

I.44035E->03 

-4.40228E-04 

1.93474E-03 

2.0441SE-03 

1 08 . 82 

24 

1.71187E»03 

3.94011E-03 

3.72S49E-03 

7.03026E-03 

32.03 

23 

1.7B320E-03 

-1 .S8104E-03 

3.92221E-04 

1.47208E-03 

142.32 

26 

1.8S433E-03 

-1 . 47615E-03 

-1.79043E-04 

1  .  48497E-03 

-173.08 

27 

1  . 92586E-03 

-4.41444E-04 

-4.33421E-04 

7.73349E-04 

-143.83 

2S 

1.94719E.03 

-8.13307E-03 

9.00M3E-03 

1 .21329E-04 

132. 11 

29 

2.068S2E-0  3 

-2.07432E-04 

8.13439E-04 

8.41422E-04 

104.27 

30 

2.13984E»03 

-8.771 18E-04 

1.2I319E-03 

1 .31331E-03 

123.42 

31 

2.21117E«03 

-1.42919E-03 

8.98379E-04 

1 .8334SE-03 

131.40 

32 

2.232S0E-03 

-1.17188E-03 

1 . 1 I42SE-02 

1 . 11244E-02 

94.03 

33 

2.23383E-03 

-1 .69337E-03 

-7.92281E-04 

1 . 871 34E-03 

-154.93 

34 

2.42316E-03 

-9.86436E-04 

-1 .I0Q80E-03 

1.41333E-03 

-129.84 

33 

2.49648E+03 

-I.71417E-04 

-4.82582E-04 

7.03826E-04 

-104.11 

36 

2.567S1E-03 

-3.73342E-03 

4.80402E-05 

7.4833BE-03 

140. 01 

37 

2.63914E*03 

-6.94632E-04 

3.S2488E-04 

8.87681E-04 

141 .49 

3S 

2.71047E*03 

-1 . 49144E-03 

3.19423E-04 

1.52S30E-03 

147.90 

39 

2.781S0E*03 

-1.B39S2E-03 

-3.14091E-04 

1.91035E-03 

-144.33 

40 

2.33312E-03 

2.61S39E-03 

-1.2B2S3E-02 

1 .30893E-02 

-78 .  47 

41 

2.9244SE*03 

-4.33980E-04 

-1 .44877E-03 

1 .76645E-03 

-111 . 03 

42 

2 . 94S7BE»03 

2.10233E-04 

-1 . 32474E-03 

1 .34330E-03 

-31 . 00 

43 

3.06711E-03 

3.81743E-04 

-5.3B089E-04 

4.37740E-04 

-54 .  44 

44 

3. :38*4E-03 

-3.97700E-03 

3.22S14E-03 

7.94090E-03 

13B.82 

-  -  -  43 

3.20977E+03 

-9.54843E-04 

-4.13353E-03 

8.37890E-a4 

-177.24 

44 

3.28109E»03 

-1.24383E-03 

-7.44279E-04 

1 .43122E-03 

-149. 15 

47 

3.3S242£«03 

-9.73723E-04 

-1.39941E-03 

1.97371E-03 

-121.39 

48 

49 

3.42373E->03 
3.49308E-03 

3 . 93778E-03 
7.81389E-04 

-8.84769E-03 
-1.43332E-03 

1  . 06333E-02 
1 .810B8E-03 

-34.13 

-44.43 

-35.02 

-4.41 

143.59 

-133.19 

-103.03 

-75.99 

30 

3.34641E*03 

1.11439E-03 

-7.32330E-04 

1  . 3433IE-03 

31 

32 
33 
34 
33 

3 . 63773E»03 

4.74851E-04 

-3.20037E-03 

6.76952E-04 

3.70906E-03 
3.79039E-03 
3.8S172E-03 

-4.69471E-03 
-3.78839E-04 
-3.49291E-04 

1.97127E-03 
-3.73020E-04 
-1.3733SE-03 

6.97SB9E-0S 
8.13920E-04 
1 .42214E-03 

3.92205E+03 

4.34710E-04 

-1.73343E-03 

1 . 789S4E-03 

36 
37 

3.99437E»03 

-3.49340E-02 

1 . 14931E-02 

3.874I4E-02 

142.43 
-16.74 

4. 04370E-03 

1.44736E-03 

-3.013B3E-04 

1 .74112E-03 

39 

4. I3703E*03 

1.33333E-03 

3.08470E-04 

1 . 37030E-03 

13.01 

43.  19 

171 .46 

-92.39 

39 

4.20B34E-03 

3.23430E-04 

4.93000E-04 

7.20304E-04 

60 

4.27769E+03 

-8.34732E-03 

1.23437E-03 

B . 461 12E-03 

61 

4.33102E-03 

-3.47331E-05 

-7.69708E-04 

7  70491E-04 

62 

4.42234E*03 

7.31244E-04 

-1.196S7E-03 

1 . 41286E-03 

-57.89 
-29.44 

62 

4.49347E+03 

1.3449SE-03 

-B.83121E-04 

1  . 79696E-03 
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Correlation  (CXY)  results 


FREQUENCY  DOMAIN  DATA 


FREQUENCY  UINOOU-0  TO  6. 40000E»01 IHzI 
FREQUENCY  INTERVAL-  1  , OOOOUE'OO [Hi  1 


'  COEFF . 

FREQUENCY; Hi! 

REAL 

MM 

MAGNITUDE 

RHASEtDEGl 

DC  TERM 

0.90000E-00 

3.21394E-01 

nax  FREQ. 

6.40000£*01 

-2.18730E-U 

1 

1 . 0OO00E+00 

3.43436E-U 

3.63294E-U 

4.62483E-11 

38.37 

2 

2.00OO0E-O0 

-8.84103E-12 

3.30774E-11 

3.41794E-11 

104. 18 

3 

3.00000E*00 

-1  .80304E-U 

1.88733E-11 

2.61032E-11 

133.49 

4 

4.00000E+00 

-1.33038E-1I 

1.23183E-U 

1 . 827B3E-1 1 

137.43 

3 

3.50000E-00 

-6.S2333E-12 

9.84300E-12 

i.iattzs-n 

123.34 

4 

6.00000E-00 

-1 .22730E-11 

-1.97721E-I2 

1 .24313E-11 

-170.83 

7 

7. 0O0O0E*OO 

9.89933E-12 

3.48738E-12 

1.1318SE-U 

29.00 

3 

3. Q0000E*00 

-1 . 96406E-1 1 

9.32094E-12 

2.18267E-1! 

134. 14 

9 

9. 000O0E+00 

-9.60420E-12 

-1.23714E-11 

1  .36630E-1I 

-127.83 

10 

1 .OOOOOE-01 

-1.93432E-I1 

-1.03196E-11 

2.13192E-11 

-130.43 

11 

1.10000E*01 

-1.90448E-12 

-1.66749E-11 

1.47833E-11 

-96 . 32 

12 

1 .20000E-01 

1.26353E-13 

2.30401E-11 

2.90604E-1I 

89.74 

13 

1  .30000E*01 

-4.4I009E-U 

-3.96644E-11 

4.01484E-U 

-140.01 

14 

1 .  40oaoe*oi 

3.69434E-1I 

-3.33326E-11 

6.66984E-11 

-34.37 

13 

1 .30000E-01 

3.13103E-U 

8.78624E-13 

3.13227E-1I 

1.41 

11 

1  . 40000E-01 

4.79688E-U 

-1.87300E-U 

7.03073E-11 

-13.42 

17 

1.70000E-01 

3.10930E-11 

2.30238E-11 

3.84904E-U 

36.32 

ia 

1 ,80000E»01 

-3.03864E-13 

2.80  092E-U 

2.80109E-U 

90.42 

l? 

1 . 90000E-01 

-2.23839E-12 

1  .43991E-11 

1.47728E-11 

98.79 

20 

2.00000E*01 

-1.03219E-11 

1.36381E-11 

1 .88483E-11 

123.93 

21 

2.  '.3000E-01 

-1.77134E-11 

-1.00770E-U 

2.03772E-U 

-130.36 

22 

2.20000E»0! 

1 -48419E-11 

-1.73246E-I2 

1  .69308E-U 

-3.37 

23 

2.30000E.01 

-1.42040C-U 

3.02S34E-U 

3.34317E-U 

113.13 

24 

2.  40000E-01 

-1.37220E-11 

-1.38373E-I1 

2.23137E-11 

-134.77 

23 

2.30000E-01 

-4.2463BE-14 

3.21803E-12 

3.21842E-12 

90.49 

26 

2.60000E»01 

-6.34926E-I2 

2.21849E-14 

6.34930E-12 

177.80 

27 

2.7ooooe»oi 

-9.78391E-12 

7.90984E-12 

1.23814E-U 

141.03 

28 

2.80000E+01 

-9.40631E-I2 

-1.S8281E-U 

1 .43431E-11 

-130.98 

2? 

2.9ooooe+oi 

2.91247E-U 

-1.04332E-11 

3.19438E-11 

-19.74 

30 

3.QO0OOE*')! 

2.S9814E-12 

1.8397IE-U 

1.83777E-11 

91.94 

31 

3. 10000E-01 

-1.08310E-U 

1.29749E-U 

1 .49030E-11 

129.83 

32 

3.23000E-01 

-2.63423E-U 

1.36230E-12 

2.  66384E-1 I 

176.43 

33 

3.30000E»01 

-I.39132E-I1 

-1  .31612E-U 

2.1780SE-H 

-136.39 

34 

3. 4OOO0E-O1 

-1 .40990E-12 

-1.42334E-U 

1 .63150E-1I 

-93.64 

33 

3.30000E-01 

3.36480E-12 

-1.20273E-I1 

1.32323E-11 

-43.17 

36 

3.60000E»01 

1.8013IE-11 

2.83336E-12 

1 .82377E-11 

9.  00 

37 

3.70000E-01 

-4.03243E-I2 

1.73009E-U 

1  . 77446E-1 1 

103.12 

38 

3.30000E-01 

-1.32042E-11 

-3.420712-13 

1.32142E-U 

-177.34 

3T 

3.90QO0E+Q1 

-1.33241E-U 

-1.04938E-11 

1.71194E-U 

-142.20 

40 

4.00000E*01 

1  . 02134E-11 

-3.20374E-11 

2.42892E-I1 

-63. 13 

41 

4. :aooaE+oi 

6.21642E-12 

2.37291E-11 

2.43299E-11 

73.32 

42 

4,;ooooe*oi 

-3.39422E-11 

-6.11333E-I2 

3.64384E-U 

-170.33 

43 

4.30000E+01 

9.38641E-12 

-6 . 48287E-1 1 

6.74847E-I1 

-92.00 

44 

4.  40030E  +  01 

3.44S77E-11 

-4.88726E-U 

7.46726E-U 

-40.88 

43 

4.30000E-01 

8.13342E-U 

4.49011E-U 

9.3097BE-U 

29.84 

44 

4.40000E-01 

-9.03427E-12 

4.34026E-U 

4.40264E-11 

97.98 

47 

4.73000E-01 

8.73789E-12 

-3.36321E-I2" 

9.36466E-12 

-21.03 

48 

4.30003E-01 

4.98097E-01 

4.3S779E-02 

3. 000OOE-01 

5.  00 

49 

4.93000E+01 

-6.5494SE-U 

-2.27697E-11 

6.93397E-I1 

-140 .83 

30 

3. OQOOOE'01 

-2.12238E-11 

-1.42044E-11 

2.53413E-11 

-144.  21 

31 

3. 13000E-31 

-3.33S31E-11 

-2.99243E-12 

3.54780E-11 

-173.20 

32 

3.20000E-01 

-4.43890E-U 

-4.31932E-11 

6. 19358E-11 

-133.78 

33 

3. i0003E*01 

I.94428E-1! 

-6.90197E-11 

7.17U3E-U 

-74.23 

34 

3. 400QOE-01 

3.88924E-11 

-1 .04234E-I2 

3.89043E-1! 

-1  .34 

33 

s.sooooe-oi 

4.33378E-12 

-1.63883E-U 

1.72023E-1! 

-74.44 

36 

3.40000E-01 

2.82742E-11 

-1 .1*7731-11 

3.03914E-11 

-22.  44 

37 

5.70000E-01 

2.02272E-U 

-7.99883E-13 

2.  02430E-U 

-2.26 

38 

3.30003E-0  1 

2.B8099E-11 

1.18144E-U 

3.U3S3E-U 

22.  30 

39 

3.90000E-01 

1.12313E-I1 

2.39436E-11 

2.443S4E-U 

64.83 

60 

4.00000£*01 

1 .97191E-12 

2.42373E-11 

2.43373E-1I 

83.  33 

61 

4.10000E-01 

-1.23739E-11 

2.33479E-11 

2.43194E-U 

118.31 

62 

6.20000E*Ol 

3.28107E-12 

7.39344E-12 

8. 09a6aE-12 

66.  07 

63 

4.3oaooe*oi 

-1.331921-11 

3.60833E-11 

3.84600E-11 

1 1 0 . 23 
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APPENDIX  7 
Program  listing  and  results  obtained  for  the  Numerical  Method. 

tO  !  THIS  PROGRAM  DETERMINES  THE  PHASE  BY  SIMULATION  OF  TUO  SIGNALS 

12     !  2iKi.,SNS.I(l<>  BY  FITTING  A  THIRD  DEGREE  POLYNOMIAL  AT  THE  ZERO  CROSSINGS 

30     !  OF  THE  SIGNALS  AND  USES  THE  NEUTON'S  METHOD  TO  FIND  THE  ROOT  AND  THEN 

40     !  CALCULATES  THE  PHASE...  SHOUS  ALSO  THE   ERROR  BETWEEN  THE  ACTUAL  PHASE 

50     !  AND  THE  CALCULATES  PHASE. .. SHOUS  ALSO  THE  EFFECT  OF  VARYING  THE  RATIO 

40     I  OF  SAMPLING  FREOUENCY  TO  THE  FREQUENCY  OF  THE  SIGNALS .... SHOUS  THAT 

70     !  THE  ERROR  IS  MINIMUM  AT  A  RATIO  OF  6. 

80     !  GENERATE  Y1(K>  AND  Y2<K>  THE  TUO  SIGNALS  OF  INTEREST 

90     !  Y2  LEADS  Yl 

100  DIM  YK200)  ,Y2<200>  ,X(200> , I<30> , J<30> 

110   I  SAMPLING  FREQ.  IS  Ft  AND  THE  FREOUENCY  OF  THE  SIGNALS  IS  F 

120   !  BY  VARYING  Fl  UE  CAN  STUDY  THE  EFFECTS  OF  THE  RATIO  Fl/F. 

130  F«100 

140  INPUT  Fl 

ISO   PRINT  "SAMP. FREQ. IS  ",F1 

160    PRINT  -FRE8.  OF  SIGNAL- ",F 

170  FOR  P»0  TO  90  STEP  10 

180    I  HERE  P  IS  THE  PHASE  IN  DEGREES 

190  F0"2»(22/7)»F 

200  T-l/Fl 

210    P0"P»<22/7)/180 

220    FOR   K"l    TO    199 

230    Y1(K)»SIN(F0»K»T) 

240    Y2(K)»SIN<F0»K»T*P0) 

230   NEXT   K 

2J0   Ll-1 

270    L2"t 

2S0   P3-0 

290      !    DETERMINATION   OF  ZERO   CROSSINGS 

300  FOR   K«l    TO    197 

310  IF  SGN<Y1CK)>-SCN<Y1<K+1>>  —  2  THEN  GOTO  330 

320  GOTO  3S0 

330  HLD-K 

340  Ll'Ll-M 

33 »  NEXT  K 

360  FOR  K-l  TO  197 

370  IF  SGN<Y2<K))-SCN(Y2<K-M>>»-2  THEN  GOTO  390 

380  GOTO  410 

390  J(L2>»K 

40  0  L2-L2-M 

410  NEXT  K 

420    !  FIND  THE  PHASE  AT  THE  FIRST  10  ZERO  CROSSINGS  AND  TAKE  THE  AVERAGE 

430    i    FIRST  SET  UP  THE  REBD.  VALUES  FOR  THE  DIFFERENCE  TABLE 

440  FOR  N"l  TO  10 

430  K-ICN) 

460  X(K)«Y1(K) 

470  XCIC-l)-Yl(K-l) 

480  X<K>1)»Y1<K*1> 

4»0  X(K*a>"YKK*2) 

300  I  USE  SUBROUTINE  FOR  FITTING  THE  POLYNOMIAL 

310  GOSUB  630 

520    1  HERE  Tl  DENOTES  THE  ZERO  CROSSING  OF  FIRST  SIGNAL 


330 

T1-X0 

340 

!  REPEAT 

FOR 

SECOND 

SIGNAL 

S30 

K«J(N) 

560 

X<K)»Y2<K> 

370 

X<K-1)-Y2CK 

-1  > 

380 

X(K*1)»Y2CK 

♦  1  > 

190 

X(K*2)-Y2(K 

♦  2) 

600 

GOSUB  650 
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610    I  T2  REPRESENTS  ZERO  CROSSING  OF  SECOND  SIGNAL 

620  T2»X0 

630  GOTO  990 

640    !  SUBROUTINE  FOR  FITTING  THE  POLYNOMIAL  USING  THE  DIFFERENCE  TABLE 

630  REN  FORM  DIFFERENCE  TABLE 

660  IF  ABS<X(K)X.001  THEN  XCK)»0 

670  IF  A8S(X(K+1)X.001  THEN  X<K»1)"0 

680  AO-X(K-l) 

670  A1"X(K> 

700  A2-X<K-M> 

710  A3-X(K*2) 

720  B0-A1-A0 

730  B1-A2-A1 

740  B2-A3-A2 

730  C0-<Bl-B0>/2 

760  Cl»<B2-Bl)/2 

770  DO"<Cl-C0>/3 

780  M1»F0»T 

7?0  U2"U1»U1 

800  U3"U2»U1 

810  E0-A0-B0»<K-1  )<tUl+C0«K»<K-l  >»U2-D0«K»<K-1  )»(K*1  )»U3 

820  E1"B0-C0»(2»K-1 )«M1*D0»<3»K»K-1 )»«2 

830  £2»C0-D0»<3»K>»U1 

840  E3-D0 

830  REM  THE  ABOVE  ARE  THE  COEFFS.  OF  POLYNOMIAL. 

360  REM  NEUTONS  METHOD 

870  DEF  FNP(V>"E0+E1»V+E2*V»V+E3»V*V»V 

880  DEF  FNa(V)»ei*E2»V*2i-E3»V»V»3 

890  X0"«l»K»Ul/2 

900  X0-X0-FNP<X0)/FN9(X0) 

910  IF  ABS(FNP(X0))-tE-4<0  THEN  930 

920  M1«N1+1 

930  IF  Ml>10  THEN  770 

740  GOTO  900 

930  H1»0 

960  RETURN 

970    !  OBTAIN  THE  DIFFERENCE  BETWEEN  THE  ZERO  CROSSINGS  OF  THE  TWO  SIGNALS 

980    !  AND  CALCULATE  THE  PHASE 

990  P1-T1-T2 

1000   P2«P1/U1»360»F/F1 

1010   P3-P3+P2 

1020    NEXT    N 

1030  PRINT  'AVERAGE  OF  THE  PHASE  VALUES  IS-,P3/(N-1) 

1040  PRINT  -ACTUAL  PHASE  IS", P, "ERROR  IS  -,P-P3/(N-1) 

1030   PRINT  •  • 

1060   P3-0 

1070   I  REPEAT  FOR  NEXT  VALUE  OF  PHASE 

1080  NEXT  P 

1090  END 
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SAMP.FREQ.IS  600 

FREQ.  OF  SIGNAL=      100 
AVERAGE  OF  THE  PHASE' VALUES  IS 
ACTUAL  PHASE  IS       0 

AVERAGE  OF  THE  PHASE  VALUES  IS 
ACTUAL  PHASE  IS       10 

AVERAGE  OF  THE  PHASE  VALUES  IS 
ACTUAL  PHASE  IS       20 

AVERAGE  OF  THE  PHASE  VALUES  IS 
ACTUAL  PHASE  IS       30 

AVERAGE  OF  THE  PHASE  VALUES  IS 
ACTUAL  PHASE  IS       40 

AVERAGE  OF  THE  PHASE  VALUES  IS 
ACTUAL  PHASE  IS       50 

AVERAGE  OF  THE  PHASE  VALUES  IS 
ACTUAL  PHASE  IS       60 

AVERAGE  OF  THE  PHASE  VALUES  IS 
ACTUAL  PHASE  IS       70 

AVERAGE  OF  THE  PHASE  VALUES  IS 
ACTUAL  PHASE  IS       80 


AVERAGE  OF  THE  PHASE 
ACTUAL  PHASE  IS 


VALUES  IS 
90 


ERROR  IS 

9.77897919153 
ERROR  IS 

19.3094805375 
ERROR  IS 

28.7222999216 
ERROR  IS 

38.1637038994 
ERROR  IS 

47.7834279363 
ERROR  IS 

59.9999991558 
ERROR  IS 

69.7789819638 
ERROR  IS 

79.3094892334 
ERROR  IS 

88.7222917675 
ERROR  IS 


.2210208085 
.6905194625 
1 .2777000784 
1 . 8362961  006 
2.2165720637 
. 0000008442 
.2210180362 
.6905107666 
1.2777032325 
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APPENDIX  8 
Listing  of  the  program  used  for  the  Adaptive  Predictor. 

0  REM  THIS  PROGRAM  IS  AN  ADAPTIVE  PREDICTOR  WHICH  ACTS  AS  A 

30  REM  SIGNAL  ENHANCER. 

4(9  REM  THE  INPUT  SIGNAL  IS  STORED  IN  THE  ARRAY  X)  WHICH 

50  REM  CONTAINS  500  SAMPLES  OF  THE  INPUT  SIGNAL  X 

BO  REM  THE  ARRAY  Xl(*>  IS  USED  TO  GENERATE  THE  N  PAST  VALUES 

70  REM  OF  THE  INPUT  AND  THIS  ACTS  AS  THE  INPUT  TO  THE  FILTER 

80  REM  THE  ARRAY  G(*>  IS  USED  TO  STORE  THE  OUTPUT  OF  THE 

90  REM  FILTER  WHICH  IS  AN  ENHANCED  VERSION  OF  THE  INPUT 

100  REM  SIGNAL  -  GETS  RID  OF  THE  NOISE  IN  INPUT  SIGNAL 

110  REM  SINCE  THE  ADAPTIVE  FILTER  TAKES  A  FEW  ITERATIONS  TO 

1£0  REM  ADAPT  ITSELF  .  THE  OUTPUT  G<*>  MUST  BE  TAKEN  AFTER  A 

130  REM  FEW  ITERATIONS  FOR  SUBSEQUENT  PROCESSING. 

140  REM********************************************************** 

150  REM  FORM  THE  DIFFERENT  ARRAYS 

160  REM********************************************************** 

170  DIM  X(500),  X1(16),G<500),  W(16> 

1 80  REM***************************************************** 

190  REM  INITIALIZE  FILTER  COEFFICIENTS 

£00  REM***************************************************** 

£10  FOR  1=1  TO  16 

£20  W(I)=0 

£30  NEXT  I 

£40  REM***************************************************** 

£50  REM  DEFINE  THE  NUMBER  OF  WEIGHTS  OF  THE  FILTER  TO  BE 

£60  REM  USED  (N) 

£70  REM***************************************************** 

£60  N=10 

£90  rem***************************************************** 

300  REM  DEFINE  THE  CONVERGENCE  PARAMETER  VI 

310  REM***************************************************** 

3£0  V1=.00£ 

330  REM***************************************************** 

340  REM  START  PROCESSING 

350  REM***************************************************** 

360  FOR  K=N  TO  499 

370  rem***************************************************** 

380  REM  COMPUTE  INPUT  SIGNAL  ENERGY  AND  USE  IT  TO  UPDATE 
390   REM  THE  CONVERGENCE  PARAMETER 

400  REM***************************************************** 

410  VS=V£*(1-V1)+V1*X  (K)*X  (K) 

4£0  IF  V2<.001  THEN  450 

430  V=V1/V2 

440  GOTO  470 

MKJ  V=0 


iai 


460  REM***************************************************** 

470  REM  GENERATE  XI  THE  INPUT  TO  THE  FILTER 

480  REM***************************************************** 

498  FOR  J=l  TO  N 

508  M=J->1 

510  XKJ)=X  (K-M) 

5S0  NEXT  J 

530  G0=0 

548  REM***************************************************** 

550  REM  FIND  THE  OUTPUT  OF  THE  FILTER  G0 

560  REM***************************************************** 

570  FOR  L-  1  TO  N 

580  G1=W(L)*X1 (L) 

590  G0=G0+G1 

600  NEXT  L 

610  REM***************************************************** 

6£8  REM  COMP.  ERROR 

630  REM***************************************************** 

640  E=X(K+1)-G0 

650  REM***************************************************** 

660  REM  STORE  THE  OUTPUT  OF  THE  FILTER  IN  G<*> 

670  REM***************************************************** 

680G<K+1)=G0 

690  FOR  1=1   TO  N 

700  REM***************************************************** 

710  REM  UPDATE  COEFFS.  USING  THE  ERROR  E 

7£0  REM***************************************************** 

730  U(I)=W(I)+V*E*X1(I) 

740  NEXT  I 

750  NEXT  K 

768  REM***************************************************** 

770  REM  AFTER  THE  REQUIRED  NUMBER  OF  ITERATIONS  G(*) 

780  REM  CONTAINS  THE  ENHANCED  OUTPUT ALLOW  ABOUT  100 

790  REM  ITERATIONS  FOR  THE  ADAPTATION  OF  THE  FILTER. 

800  REM***************************************************** 

818  RETURN 
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APPENDIX  9 
Pr-onr-am  listing  for  the  PLOT  package 


NAME   PRO  TYPE  RIO     REM  THIS  PROGRAM  SHOUS  HOU  THE  Plot  SUBROUTINE  CAN  BE  USED 

20  REM  TO  PLOT  THE  REQUIRED  DATA 

30  REH  HERE    RANDOM  IS  TO  BE  PLOTTED 

40  OPTION  BASE  t 

50  DIM  Ho<100>,V*UDQ> 

60  RANDOMIZE 

70  INPUT  Xpoints 

80  FOR  1-1  TO  1 00 

90  Ho(I)»I 

100  V»U>"SIN<2»PI«10«I> 

110  NEXT  I 

120  EOSUB  Plot 

130  END 

140  Ploti   REM    THIS  IS  A  SUBROUTINE  TO  PLOT  DATA  IN  THE  ARRAYS 

130  REM  Ho  AND  Vs  .  THE  PROCEDURE  IS  TO  FILL  THE  ARRAYS  AND  THEN 

160  REM  CALL  THE  SUBROUTINE  Plot 

170  REM  THE  DATA  MAY  BE  ON  FLOPPY  DISK  OR  OTHER  STORASE  DEVICE 

ISO  REH  THIS  DATA  HAS  TO  BE  ACCESSED  AND  USED  TO  FILL  THE  ABOVE 

1?0  REM  MENTIONED  ARRAYS 

200  REM  Find  A4I  *nd  «ln  valuss  of  ths  data  In  ths  arrays 

210  Ho«ln-Vswin*lE99 

220  REH  Npoints  IS  THE  NUMBER  OF  POINTS  TO  BE  PLOTTED 

230  FOR  1-1  TO  Npoints 

240  IF  HodKHonln  THEN  Honin-Ho(I) 

230  IF  VKIXVs-nln  THEN  v*nin-v»i  I  > 

260  NEXT  I 

270  Hon«»"Vs«ax«-lE99 

280  FOR  1-1  TO  Npolnts 

270  IF  Ho<I)>Honax  THEN  Honjx-Ho(I) 

300  IF  V»CI)>V6naz  THEN  V»naK-Vs(I} 

310  NEXT  I 

320  PRINT  Honin,Honai,Vsnln,Vsn4i 

330  REM  TO  find  Units  Hoi-  ,Hol  ,Vst ,  Vsp 

340  GOTO  620 

330  !  ACTUALLY  A  SUBROUTINE 

360  REM  THIS  SUBROUTINE  USES  THE  VALUE  IN  Valus  AND  RETURNS  A 

370  REM  NEU  VALUE  DEPENDING  WHICH  OF  THE  LIMITS  Hor ,Hol , Vst ,V»t> 

380  REM  IS  TO  BE  DETERMINED 

•370  IF  ABS(INT<Val»s)>-ABS<Valos>«0  THEN  Return 

400  IF  ABS<Valu»>>t  THEN  Thsrs 

410  IF  yaluo-1  THEN  Rsturn 

420  GOTO  430 

430  Dscinal-1 

440  Dplaca«l»Dsclnal 

430  DplacsM/Dplacs 

460  IF  ABS<Valu»)>Dplacs  THEN  Hsrs 

470  Dsclnal-DscinalalQ 

480  GOTO  440 

490  Thsrsi   Dplac*"l 

S00  Hsrsi   Mark-Dplacs 

310  I 

320  IF  Lln-1  THEN  Eotu<   • 

330  IF  Lln-0  THEN  Getnln 

S40  Eotnn ;   Hsiad-Valo*  DIV  Jilititl 
330  IF  Moxod-0  THEN  Hsiad— 2 

560  Valos-MaiadsDplacs 

570  SOTO  Rsturn 

580  Gstnini    Mlnad-Valus  DIV  Do  loco- 1 

390  IF  Hinad-0  THEN  Mlnad- . 9 
600  ValoS'HlnadsDplacs 


611    Rbtarn I         RETURN 

620  !    Find    the   Unit* 

630  !     Hor 

640  Valv«>Haff«l 

690  LiM"l 

660  GOSUB   3S0 

670  Hor-v*iua 

680  I    Hal 

690  yalua-Homn 

700  LiM-0 

711  COSUB    3S0 

720  Hal-Ualub 

730  I    V*T 

740  V41U6-V6H4I 

730  Lix-1 

760  COSUB  3S0 

770  U«1aVllu« 

780  !  U»b 

790  Valua-vamn 

800  Lin-0 

'  810  GOSUB  330 

820  Ubb-Valtt* 

330  REM  THE  TIC  SPACINGS  FOR  DRAUINC  THE  AXES  ARE  THEN  FOUND 

840  !  FIND  XTlcTTlc  ETC 

8S0  HoTic-(Hor-Hol)/20 

860  Hn4j"S 

870  UtTic-<y»T-v»o>/20 

880  Wp)aJ"3 

890  PRINT  Hal,Hai-,U.b,Ua>T,Hatlc,V«Tle 

900  REM  THESE  SUBROUTINES  DETERMINE  THE  DIFFERENT  PARAMETERS  TO  BE 

910  REM  USED  IN  FINDING  INTERSECTIONS  AND  THE  SPACINGS  NEEDED  FOR 

920  REM  LABELING  THE  AXES 

930  IF  <H»1<0>  AND  (Har<0>  AND  IIWIO)  AND  (U»t<0>>  THEN  GOSUB  On. 

940  IF  1*11101  AND  <Hel->0>  AND  ( < Uao < 0 >  AND  (Willi)  THEN  COSUB  T«a 

930  IF  (Hal)O)  AND  (Har)O)  AND  ((Va>b(0)  AND  (V»T(0))  THEN  GOSUB  Three 

960  IF  (HaKO)  AND  CHor(O)  AND  ((U*b<0)  AND  (V«T>0))  THEN  GOSUB  Faur 

970  IF  (HaKO)  AND  (Hei->0)  AND  ((Va}b(0)  AND  <V9T)0))  THEN  COSUB  Fiv» 

980  IF  (Hal>0)  AND  (Har>0>  AND  ((Vbb(O)  AND  (V»f)0))  THEN  GOSUB  Six 

990  IF  (HaKO)  AND  (Har(O)  AND  ((Ubb)O)  AND  (U»T)0))  THEN  GOSUB  S«w«n 

1000  IF  (HaKO)  AND  (Har)Q)  AND  ((U»b>0>  AND  (U«T)0))  THEN  GOSUB  EignT 

1010  IF  (Hal)O)  AND  (Har)O)  AND  ((U»b>0)  AND  (U.t>0))  THEN  GOSUB  Nln« 

1020  PRINT  Xlnt.Tint 

1030  FOR  I  "5  TO  20  STEP  3 

1040  PRINT  Hal+Hetical 

1030  PRINT  Vbb*V«Tic»I 

1060  NEXT  I 

1070  REM  START  PLOTTING  OPERATIONS 

1080  PLOTTER  IS  7,3,-9872A* 

1090  GRAPHICS 

1100  LIMIT  30,300,30,200 

1111  FRAME 

1120  LOCATE  23,143,23,90 

1130  FRAME 

1140  SCALE  Hal,Har,V«b,U«T 

1130  AXES  Hatic,U»Tlc,Xint,Yll)t,Hn4j,UmJ,3 

1160  FRAME 

1170  REM  ACTUAL  PLOTTING  OF  DATA 

1180  FOR  1-1  TO  Npalnt* 

1190  PLOT  H.III.HC) 

1200  NEXT  I 


1210      SOTO    1540 

1220   On.:Xint-Hor 

1230      Ytnt-M.T 

1240      RETURN 

12S0    Two:    Xjfit-0 

1260       TllfVlt 

1270      RETURN 

1280    Thr»»:    XlnT-Hol 

1270      UnT-V»t 

1300      RETURN 

1310    Fbbp;     Xlnt-Hor 

1320      Tinr-0 

1330      RETURN 

1340   Fiv«>    Int.rini.o 

1330      RETURN 

1360   Six :    Xint"Hol 

1370      Tini.l 

1380    RETURN 

1370   3>«>ni    xlnT-Hor 

1400       Tlnr-VW 

1410      RETURN 

1420    Eight!    Xtnt»0 

1430      Ylnt-VM 

1440      RETURN 

1450    Mint:     XinT-Hoi 

1460      Ylnt-U«b 

1470      RETURN 

MJJ  |S    0FHTHEU"^N%0°J   KtSS    T„E'ftXEESEBMINE   «    ""™   ~ " 

1500  IF   HoD-0   THEN  Sgnx— 1 

1510  IF   HiK.O    THEN    Sgnx-1 

1320  IF   V»<-0   THEN   Sgny-1 

1330  IF    Urt>-0    THEN    Sgny  —  1 

1360  '  MOVE   Hol,Yrt 

1370  FOR    I-t    TO    20    STEP    5 

1380  Xd-Hol»I»H«Uc 

r3?0  Ya-YinTfl.3PV«Uc»Sgny 

1600     houe  xa.Ya 

1610  LABEL  USINC    "ND.DDE'lXa 

1620  NEXT  I 

1630  FOR  I«o  TO  20  STEP  5 

1640  Xa»XlnT»3»HoUc«Sgm 

1650  ra-u«»*i,u.,lc 

1660  HOVE    Xd,To 

1670  LABEL  USINC    -MD-DDE-lYd 

1680  NEXT  I  ' 

\Hl  22VE  Hol*5*M«'i=.Tiin.4.«.rlc.Sgi), 

1700  PRINT  "XLABEL- 

1710  INPUT  B* 

1720  LABEL  B» 

1740   DO*  L*eEL  ™E  Y  *"Si  "0T«TE  K0Jlm   <*   'EN  BY  70  DECREES 

1730   LDIR  70 

Im   PSISTXiT,LA.E5L-H•,lc*s,", •i'"*5"v"ic 

1780  INPUT  At 
1771  LABEL  A« 
1800   LDIR  1 


RAD 

512  JNJTIflL!2E  pEN  DIRECTION 
REM  END  OF  PLOTTING  OPERATION 
PEN  0 
RETURN 
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ABSTRACT 

The  applications  of  signal  processing  techniques  in 
microprocessor  or  computer  based  instrumentation  were  studied  and 
a  specific  application  in  the  measurement  of  phase  was 
considered.  The  analyses  were  carried  out  in  the  discrete  time 
(OT)  domain  using  "samples"  of  the  original  signals.  Algorithms 
based  on  Adaptive  filtering,  FFT,  Correlation  and  Numerical 
methods  have  been  developed.  These  methods  were  verified  by  using 
computer  simulation  and     implemented  on  an  automated  test  system. 

The  algorithm  based  on  Adaptive  filtering  was  developed 
using  the  Widrow's  LMS  (least  mean  square)  method  in  the 
cancellation  mode.  Samples  of  the  original  signals  form  the 
primary  and  reference  inputs  to  the  adaptive  filter  and  after  a 
sufficient  number  of  processing  steps  the  phase  difference  is 
obtained  from  the  coefficients  of  the  filter.  The  FFT  method  uses 
the  information  from  the  phase  spectra  of  the  two  signals  in  the 
frequency  domain  and  determines  the  phase  difference.  The 
Correlation  approach  uses  the  amplitude  spectra  of  each  signal 
and  also  that  of  the  cross  correlated  sequence  obtained  by 
multiplying  the  two  signals  in  the  DT  domain  and  then  calculates 
the  phase  difference.  The  Numerical  method  involves  curve  fitting 
and  Newton's  method  to  obtain  an  estimate  of  the  zero  crossings 
of  the  signals  from  which  the  phase  difference  is  obtained. 

The  software  development  and  system  implementation  was  done 
on  a  North  Star  Horizon  II,  Z80  based  microcomputer  (NS)  in 
conjunction  with  an  HC6802  interface  and  a  Hewlett-Packard 
Multiprogrammer   (HPM);    the  whole  forming  an  automated   test 


system.  The  HC  6802  interface  serves  as  a  communication  link 
between  the  NS  and  the  HPM.  Data  samples  of  the  signals  are 
obtained  by  using  the  High  Speed  Analog  to  Digital  Converter  card 
available  on  the  HPM  and  these  samples  are  then  transmitted  to 
the  NS  for  subsequent  processing  by  the  algorithms  discussed 
earlier.  Computer  simulation  was  done  on  a  Hewlett-Packard  9845B 
Desktop  minicomputer. 

The  design  aspects  of  a  suitable  system  for  implementing 
these  methods  is  also  discussed.  The  role  of  the  microprocessor 
as  the  control/communicating  device  is  emphasized.  Other 
applications  in  instrumentation  are  considered.  A  general  purpose 
plot  package  was  developed  for  the  HP  9845B. 


